ON BLOW-UP "TWISTORS" 
FOR THE NAVIER-STOKES EQUATIONS IN M^: 
A VIEW FROM REACTION-DIFFUSION THEORY 



V.A. GALAKTIONOV 

Abstract. Formation of blow-up singularities for the Navier-Stokes equations (NSEs) 

Uf + (u • V)u = -Vp + Au, divu = in x R+, 

with bounded data Uq is discussed. Using natural links with blow-up theory for nonlinear 
reaction-diffusion PDEs, some possibilities to construct special self-similar and other 
related solutions that are characterized by blow-up swirl with the angular speed near the 
blow-up time (this represents simplest w-limits of rescaled orbits as periodic ones) 

(f{t) crln(r-t) 93(i)-^->oo as t ^ T- {(J^Q). 

This is done in cylindrical polar coordinates {r, ip, z} in , using the restriction of the 
NSEs to the linear subspace W2 = Span{l,2:}. Similarly, blow-up twistors with axis 
precessions in the spherical geometry {r, 9, ip} are introduced. 

It is shown that other blow-up patterns (a "screwing in tornado" ) may correspond to a 
slow "centre-stable manifold-like drift" about Slezkin-Landau singular or other equilibria 
of the NSEs. Some approaches to blow-up singularities can be applied to 3D Euler's 
equations and to well-posed Burnett equations in 7D (i.e., the NSEs with A 1-^ — A^). 
Though most of blow-up scenarios were not justified even at a qualitative level, the author 
hopes that the proposed approaches to families of blow-up and other patterns, including 
those with blow-up swirl, will give some extra insight into the micro-scale "turbulent" 
structure of the NSEs. 

The discussion of possible types of blow-up patterns for the NSEs is going in conjunc- 
tion with some other classic nonlinear PDEs of mathematical physics. 



1. Introduction: Navier-Stokes equations and blow-up 

1.1. Two faces of the open problem via blow-up formulation: first discussion 
around self-similar blow-up. It is well- understood (see e.g., [148i Ch. 5]) that the 
fundamental open problem of fluid mechanic^ and PDE theory on global existence or 
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1.1^ 



nonexistence of bounded smooth solutions of the Navier-Stoke^ equations (the NSEs), 

ut + (u ■ V)u = -Vp + Au, 
divu = in x {u={u,v,w)), 

with arbitrary bounded divergence-free LF'-data uq, is two-fold: 

From a standard evolution PDE and blow-up point of view, there exist two possibility: 
the positive answer, i.e., existence of a global bounded solution (then necessarily it is 
unique by the existing local semigroup of bounded smooth solutions) is equivalent to 

fl.2) 



3 global bounded solution = ^ finite-time blow-up at any T > 0. 



On the other hand, the negative answer, i.e., nonexistence in general of a global 
bounded solution, is equivalent to the following: 



1.3) ^ global bounded solution = 3 a finite-time blow-up pattern, 



which corresponds to bounded L^-initial data. In both cases, we mean blow-up in L°°(M^) 
of bounded smooth solutions at the first blow-up time t = T. 

It seems obvious that both scenarios (11.21) and (11.31) assume a detailed study of possible 
blow-up behaviour of solutions a.s t —>■ T", and this cannot be avoided and represents a 
clear alternative face of this fundamental open problem of PDE theory. Of course, (11.21) 
would be achieved without any blow-up study provided that a new technique (e.g., a new 
conservation law and/or a monotonicity formulae) could be invented for the NSEs. How- 
ever, the long history of this Millennium Open Problem suggests that this is expected to 
be extremely difficult. Anyway, this can happen since the Navier-Stokes equations can 
indeed inherit from their universal physical nature (a system comprising Newton's Second 
Law, the continuity, and a basic viscosity) some extra still hidden and unknown continu- 
ous, discrete, or other symmetries/monotonicity/symplectic/dissipative, etc. features. 

It is definite that, during a few last years, the direction of the attacking this open 
problem was clearly partially changed and a seriously increasing number of papers along 
the evolution blow-up scenarios (11.21) and (11.31) were published. In particular, a rather 
complete negative answer was achieved supporting somehow (II. 2p in the following way: 



1.4) for (11.11) . a standard self-similar blow-up as t ^ T is impossible. 



In a most general evolution setting, this is due to the work by Hou and Li [ 103] , based on 
the crucial nonexistence result in Necas-Ruzicka-Sverak [162] proved via the Maximum 
Principle (the MP). The ban (11.41) was proven in |lU3j by solid semigroup theory that is 
adequate to classic evolution approaches to nonlinear PDEs with blow-up (we present a 
spectral discussion of (11.41) in Section [2]). In fact, it was proved that any convergence in 
L^(M^), with p > 3, as t — T~ to a rescaled self-similar profile means that the solution 
remains bounded at t = T (so it is a removable singularity). This approach does not 
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exhaust all the possibilities of, say, almost self-similar or other (Type II) blow-up, but, 
nevertheless, is a new convincing fact on possible singularities in the fluid model fll.ip . 

The negative result (11.41) completes a long remarkable history of the study of similarity 
blow-up singularities for the NSEs that was initiated by J. Leray in 1933-34 [136[ I137j . 
who actually posed a deeper problem on both backward and forward phenomena: 

Leray's blow-up scenario: self-similar blow-up as t ^ T~ (t < T) 
and similarity collapse of singularities as t — > (t > T); 

see his precise statements and a discussion on these principal issues in Section [2l 

1.2. Similar open regularity problems for = 4 and for the well-posed Bur- 
nett equations. It is important to recognize that, unlike its great sounding and well- 
established reputation even among non-experts, the NSEs Millennium problem is just a 
"remnant" of a general fundamental difficulty of PDE theory in the twenty first century, 
which seems cannot principally be understood by classic modern techniques, to say noth- 
ing about a rigorous proof. Therefore, most general new concepts covering a wider PDE 
area are in great demand. 

To justify that, one does not need to address essentially other classes of nonlinear 
PDEs and/or systems, and can just slightly generalize the classic NSEs. First, let us 
mention that a similar and not less difficult (but indeed more exotic) open problem on ex- 
istence/nonexistence of global classical solutions persists for the Navier-Stokes equations 
in dimension = 4; see Scheffer (1978) [188j and recent developments in |48j . 

Second, as a next neighbouring example of a similar nature associated with applications, 
consider the well-posed Burnett equations 

(1.6) ut + (u • V)u = - Vj9 - A^u, div u = in x M+, 

where, in comparison with the standard model (11. II) . the Laplacian squared with the 
correct sign — on the right-hand side is posed. Note that here one deals with the linear 
solenoidal bi-harmonic flow induced by the operator + A^, which is parabolic but 
assumes no any order-preserving and other properties of the semigroup e~^ * (which are 
somehow naturally partially inherited from the positive Gaussian for — A in the NSEs 
(11. ip ). Surely, both problems, the NSEs (11.11) in and the Burnett equations (11.60 in M7 
look rather ridiculous and seem cannot have a real application. But for PDE theory these 
can be considered as some key and canonical representatives exhibiting the necessary 
principal difficulties (or course, there are other examples like that, e.g., supercritical 
nonlinear Schrodinger or Ginzburg-Landau equations to be used and discussed as well). 

The more complicated model (11.60 appears on the basis of Grad's method in Chapman- 
Enskog expansions for hydrodynamics. In particular, equations (11.60 were studied in [71] 
as a particular case of Kuramoto-Sivashinsky-tjpe PDEs (see on a derivation therein also), 
where it was shown that, for the smooth orbits, the following embedding holds: 

(1.7) Lf(M^) =^ L~(R^) {N>Q). 
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(1.5) 



This is an analogy for (11.11) . where the idea of a similar transition L^(R'^) =^ L°°(M^) 
goes back to Leray (1934); see a survey in the next section. Fixing 

(1.8) = 7 (or any > 7), 

we arrive at a similar open problem on existence/nonexistence of global smooth solutions. 
For A^ < 6 the proof is easier. In fact, this is the analogy of A^ = 2 for (II. ip . where 
existence- uniqueness in the Cauchy problem is due to Leray (1933) |135j (extended by 
Hopf in 1951 [lOT]) and Ladyzhenskaya (1958) for the IBVPs [I28l[n9]. Obviously, with 



such a huge growth of the order of PDEs involved and the dimension A^ of the Euclidean 
space for the Laplacian and the convective term (to say nothing about the non order- 
preserving properties of the higher-order semigroup e~^ *), a detailed and convincing 
analysis of (II. 2p and (11.30 . together with a complete description of blow-up patterns, 
seems entirely illusive and non- achievable. 

1.3. On universality of the open =^ L°° problem in PDE theory. As we have 
just mentioned above, roughly speaking, the Millennium Prize Problem, posed specially 
for the NSEs, is, in a loose sense, "non-unique", since similar open regularity problems 
(or not that lighter significance) occur for many evolution PDEs of various types. We list 
a few of them, where the difficult open mathematical aspects of global existence and/or 
blow-up are associated with the following factors: 

(i) supercritical Sobolev parameter range of the principal operator (hence, standard or 
very enhanced embedding-interpolation techniques fails), and, in fact, as a corollary, 

(ii) multi-dimensional space x G M^, with A^ > 3, at least (this leaves a lot of room 
for constructing various blow-up patterns via self-similarity, angular swirl, axis pre- 
cessions, linearization, matching, etc.). 

We now list those PDEs, where we give a few recent basic references to feel the subject. 

(I) Supercritical defocusing nonlinear Schrddinger equation (NLSE) (see |150[ 1215] ) 

(1.9) -iwt = Am- with p > ps(2) = ^ (iV > 3); 

(II) 2mth-order supercritical semilinear heat equation with absorption (m = 1 is covered 
by the MP; see [Sn] and [33], where the result in § 4 for p > ps(2m) applies to small 
solutions only): 

(1.10) = -(-A)™u - with p>ps(2m) = f^ (A^ > 2m, m > 2); 

(III) The semilinear supercritical wave equations (see |106| 1221] , as most recent guides) 

(1.11) utt = l^u - \u\'P-^u, with p>Ps(2) = f^ (A^>3). 

One can add to those "supercritical" PDEs some others of a different structure such as 
the Kuramoto-Sivashinsky equations for / = 1, 2, ... [74] 

(1.12) u, = -{-Aru + {-Ayu + lZik)dkD..iHn, |d| = l, p>po = l + ^. 

Here, po is not the Sobolev critical exponent, though precisely for p > po, L'^ L°° by 
blow-up scaling, [74, § 5]. On the other hand, a more exotic applied models exhibit similar 
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fundamental difficulties such as the following nonlinear dispersion equation (see [691 [75] 
for references and some details) 

(1.13) Ut = -D.,,[{-Ar^u]-D,,{\u\P-'u), with p>ps(2m) = f±|^. 

In view of the conservation properties for the models (11.121) and (ll.lSp . these, though 
being local, can be more adequate to the nonlocal NSEs (II. ip . than the others above. 
In most of the cases, the operator on the right-hand sides satisfying for u G C^(M^) 

(1.14) A(n) = -(-A)™M- =^ (A(u),m) = -/ - / < 0, 

is indeed coercive and monotone in the metric of L^(M^), which always helps for global 
existence-uniqueness of sufficiently smooth solutions of these evolution PDEs. For the 
NLS (II. 9p . this gives a stronger conservation laws than for the focusing equation with 
the "source-like" term = Evidently, replacing A in (11.91) and (II. lip by — (— A)"^, 

m > 2 moves the supercritical range to that in (ll.lOp . On the other hand, introducing 
quasilinear differential operators —{—A)'^\u\'^u with cr > moves the critical exponent to 
ps{2m, a) = (cr + l) j^^^^ . Similar supercritical PDEs can contain 2mth-order p-Laplacian 
operators, such as the one for m = 2, with cr > 0, 

(1.15) A(m) = -A(|An|'^An) - (A(m), n) = - / | Am|'^+2 _ j |^|p+i < 0. 

However, the lack of embedding-interpolation techniques to get L°°-bounds, which can 
be expressed as the lack of compact Sobolev embedding of the corresponding spaces for 
bounded domains Q C (this analogy is not straightforward and is used as a consistent 
illustration only) 

(1.16) H'^iQ) LP+\n) for p>ps{2m), 

actually presents the core of the problem: it is not clear how and when bounded solutions 
can attain in a finite blow-up time a "singular blow-up component" in L°°. For the 
operator in (I1.15p . a similar supercritical demand reads 

(1.17) W^^,{n) <^ L^+\n) for p>ps(4,a) = (^^15±||±^, iV > 2(a + 2). 

In the given supercritical Sobolev ranges, finite mass/energy blow-up patterns for (II. 9p - 
(11.130 are unknown, as well as global existence of arbitrary (non-small) solutions. 
It is curious that for the NSEs with the same absorption mechanism as above, 

(1.18) ut + (u- V)u = -Vj9 + Au- |u|P-%, divu = in ^ m+, 

by the same reasons and similar to (ll.lOp . the global existence of smooth solutions is 
guaranteed [27] in the subcritical Sobolev range only: for 

(1.19) p < 5 = ^^lljy^g (and p > | by another natural reason). 

We thus claim that, even for the PDEs with local nonlinearities (I1.9p -( |LTT1) (and similar 
higher-order others), the study of the admissible types of possible blow-up patterns can 
represent an important and constructive problem, with the results that can be key also for 
the non-local parabolic flows such as (11.10 . (II. 6p . (11.180 . etc. Moreover, it seems reasonable 
ffist to clarify the blow-up origins in some of looking similar and simpler (hopefully, yes, 
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since ( II .ip is both nonlocal and vector-valued unlike the others) local supercritical PDEs, 
and next to extend the approaches to the non-local NSEs fll.ll) : though, obviously, the 
former ones are not that attractive and, unfortunately, are not related to "millennium" 
issues (however, many PDE experts very well recognize how important these are for 
general PDE theory). 



1.4. Main synthetic goal of the paper and on the proposed style of research. 

For the author, who for almost thirty years dealt with blow-up singularities in various 
nonlinear PDEs and mainly in react ion- diffusion systems of different orders, the appear- 
ance of the paper |103j^ was a crucial sign. Actually, this announced that, if blow-up 
singularities in the model (II. ip are possible, these must be of interesting, complicated 
enough, non-similarity, and non-symmetric nature, though anyway a long story of var- 
ious unsuccessful attempts to reconstruct those blow-up singularities (if any) suggested 
that this is expected to be a difficult affair. This convinced us to look now carefully at 
such a complicated model at the Navier-Stokes equations (a system of four PDEs with 
{x, t) G X M+!) from the point of view of standard blow-up theory. Though the author, 
who was dealing with uncomparably simpler PDEs, which however generated a number 
of still open problems (further comments on this matter will be given later), naturally 
expected that a definite, to say nothing about a rigorous, singularity construction might 
not be convincingly done just in view of a general complexity of the model. 

Let us specify the actual main goal of this essay. Of course, clearly, some part of the 
author general motivation is associated with the Millennium (ll.2p - (ll.3p Problem. On 
the other hand, there exists the second half of the motivation, which the author honestly 
regards as not as less important (and even more valuable for general PDE theory). 

More clearly, in his opinion, performing a partial or most complete classification of 
singularities for (II. ip and other related nonlinear models of practical interest along the 
lines of blow-up theory becomes nowadays a fundamental mathematical direction. Overall, 
this is about a description of possible complicated micro configurations that an evolution 
system can create, or. 



;i.20) Goal 



for NSEs: describe all "turbulent" incompressible fluid 
configurations on possibly minimal micro-{x,t) -scales, 



regardless whether these are developed in finite-time blow-up or non blow-up manners. 
Actually, (ll.20p does not assume that the solutions must be of a finite kinetic energy; it 
is just necessary that the behaviour of solutions at infinity (as x — > oo) does not play 
any essential role for formation of patterns. Otherwise, this would mean posing various 
"boundary conditions" at infinity that can immensely increase the variety of admissible 
patterns. Describing some approaches to (11.200 is our actual main synthetic goal. 



'The first preprint of tliis paper, wliicli was available to the author, dated 24th August 2007. 
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Concerning the style of a structural characterization of our concepts and ideas, we are 
oriented to perform our research in a such unified manner that, at least, formally, 



1.21) Style : approaches cover both NSEs (11.11) and the Burnett equations (11.61) . 



This gives a room for possible future author's apology in the sense that 



1.22) if (11.21) were proved, the goal would be the 7D Burnett equations (I1.6P 



(or other complicated nonlinear evolution PDEs with similar principal difficulties of L°°- 
bounds). Of course, (11.61) will be studied less, but we will indeed seriously present basics 
of necessary related linear and nonlinear operator theory covering the case of the bi- 
harmonic operator — A^u or the 2mth-order one — (— A)™u for any m > 2. It is clear 
that any rigorous justifying the conclusions for (11.61) would be incredibly difficult or even 
will be never achieved, but anyway this cannot prevent us from performing the research: 
the necessity of a deep mathematical study of various nonlinear PDEs without any hope 
of strict formulations of many results is already clearly an inevitable feature of modern 
PDE theory associated with several types of higher- order singular or degenerate equation^. 

1.5. Layout of the paper. The main steps of our analysis are as follows: 

• Blow-up survey (Section\^. We begin with a necessary short survey devoted to classic 
and recent blow-up results for the Navier-Stokes and some for Euler equations. 

•First application of blow-up scaling: Hermitian structure of zero sets (Section\3\). This 
is devoted to a first application of blow-up scaling to the NSEs. Namely, we show that 
the local structure of multiple zeros of regular solutions can be governed by special vector 
solenoidal Hermite polynomials as eigenfunctions of the adjoint rescaled linear operator. 
We claim that this study is a natural and unavoidable step for further deeper discus- 
sion and classification of all types of micro-scale, single-point configurations that can be 
generated by the evolution system in finite time. These sets of singular patterns are as- 
sumed to include also possible L°°-blow-up patterns, which are essentially nonlinear and, 
for revealing of those, demands complicated matching procedures, where the Hermitian 
polynomial space-time structures will be key. 

• Third blow-up scaling: singularities in NSEs and EEs (Section^. This is another ver- 
sion of blow-up scaling showing that any blow-up in the NSEs must be supported by a 
bounded "NS-entropy" solution of the EEs, which is defined on larger space-time subsets. 

• Blow-up twistor mechanism (Section\B^. According to typical blow-up results, which 
are well-known for a wide audience of mathematicians working with nonlinear parabolic, 
hyperbolic, dispersion, Boussinesq, and other evolution PDEs, it is clear that, regardless 
a pretty strong negative result in (11.41) . the story of the open problem is far away not only 
from being solved but even reasonable understood. We intend to show some extra ways 
how the Navier-Stokes equations can create complicated blow-up patterns. The main 
difficulty is indeed to detect a suitable and adequate for fiuid vortex models mechanism 



^"The main goal of a mathematician is not proving a theorem, but an effective investigation of the 
problem...", A.N. Kolmogorov, 1980s . 



of blow-up swirl (rotation) leading to an essential reduction of the system (11 .ip . which 
can be at least discussed to clarify the way of construction of such blow-up patterns that 
have been called in Section [5] blow-up twistors. Namely, we show that, in cylindrical polar 
coordinates {r,(p,z}, NSEs (11.11) allow a consistent restriction to the 2D linear subspace 

(1.23) 1^2 = Span{l,z}, 

and moreover this subspace is partially invariant under the nonlinear operators involved 
(this is understood in the sense of mappings, [76, Ch. 3, 7]). Such blow-up patterns belong 
to W2 and lead to involved nonlinear systems, which are very difficult to study, but, 
anyway, this admits the natural similarity logarithmic travelling wave (TW) mechanism 
of blow-up vortex swirl about the z-axis, with the angular dependence 

(1.24) ^(t) = -aln(T-t) =^ (^(r) = ^ ^ 00 as t T- , where a ^ 0. 
In other words, we pose an extra logTW angular dependence in cylindrical coordinate^ 

(1.25) ip = fi — o"ln(T — t), where ^ G (0, 2n) is the rescaled angle, 

which inserts into the system a new "nonlinear eigenvalue" a G M. As usual, assuming 
such an extra evolution freedom in this swirling dynamical system extends the overall 
possibility to get suitable blow-up patterns, possibly even in the self-similar forn|^. In- 
serting the blow-up swirl (11.251) into (11.51) revives this Leray scenario, since now, in the 
rescaled variables, we observe not a stabilization to a point (already prohibited by the 
MP), but convergence to a periodic orbit. Roughly speaking, this falls into the scope of a 
much more difficult dynamics corresponding to the "Poincare-Bendixson Theorem" (ex- 
istence of blow-up), or to "Dulac's Negative Criterion" (nonexistence), as in classic ODE 
theory [T76J . but the current PDEs are infinite-dimensional and nonlocal. But this is not 
the end of the story even if a nonexistence Dulac-like result would have been proved: the 
rescaled orbits may converge to various quasi-periodic orbits with arbitrarily large number 
of fundamental frequencies (some details to be also discussed). 
Meantime, we state this our observation as: the NSEs 



1.26) (11.11) on W2 admit a (self-similar) mechanism of blow-up swirl at a point. 



We are not able to study somehow rigorously the evolution dynamical systems that occur 
and propose a few ideas how various blow-up twistors can occur via evolution close to cer- 
tain invariant manifolds associated with the similarity flow. Though the fluid model (11.11) 
naturally supports formation and evolution of vortices ("von Karman's streets"), it should 
be noted that, mathematically speaking, taking into account the rotational torsion-like 



^Recently, the blow-up ln(T — t) factors are more boldly appear for the NSEs; see e.g., § 1-2]- 
^Introducing "twistors" , we mean a specific self-similar angular blow-up mechanism (|1.25p and around 
with further perturbations and approximate similarity features to be described shortly. The fact that 
blow-up in the NSEs could be connected with a "tornado-type" structures is well-recognized. Q.v. Sinai 
|195|. p. 730]: "... A negative answer [to existence of strong solutions of the 3D NSEs] could be connected 
with solutions which develop singularities in finite time like a tornado- type solution where infinite vorticity 
appears at some particular points in time and space." (Underlying is author's.) 
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"spiral wave" mechanism (11.261) and other related axis and vertex precession phenomena 
includes into the nonstationary rescaled system extra velocity and other parameters, being 
nonlinear eigenvalues, that, as usual, improves the overall probability to get a necessary 
pattern by matching of various local blow-up and non-singular flows. 

We introduce and discuss these examples as a warning showing that a reasonable sim- 
ple treatment of the scenario (11. 3p cannot be expected. We postpone until Section [8] 
construction and analysis of a most involved blow-up pattern with the swirl for (II. ip in 
the spherical coordinates {r, 6, ip}, where this gets much more complicated and seems does 
not admit any lower- dimensional reductions. Overall, these lead to (ll.26p . 

Actually, our blow-up concepts are rather general and are not bounded by the frame- 
work of the classic model (II. ip . In particular, as a key ingredient, we show possible 
generating mechanisms of countable families of other blow-up patterns, which can ex- 
hibit different properties. In other words, we propose the following statement, which is 
well-understood in react ion- diffusion theory (see comments below): 



1.27) in general, "self-similarity ban" (11.41) does not prevent blow-up in NSEs (II. ip . 



For higher-order systems of PDEs with (II. 4p . the blow-up story is about to begin, 

and the present paper pretends to be just a first step along a "blow-up R-D direction". 

Of course, (11.11) is a dynamical system that admits the strong a priori control of the L^- 
norm of the solutions at the blow-up time t = T, to say nothing about the evident presence 
of the MP |162j in the stationary rescaled form of the equations, but, possibly, the system 
is complicated enough to get over such an obstacle in an evolution way. Therefore, we 
will present in Section im a first discussion of various ways how to overcame the ban (11.41) 
keeping blow-up similarity rescaled variables. We show that a thorough study of blow-up 
evolution on the quasi-stationary manifold of Slezkin-Landau's singular solutions for a 
submerged jet0 [ 196| 1132] (see also |200j ) could provide us with other types of patterns. 

Thus, the introduced blow-up twistor mechanism shows that the Navier-Stokes equa- 
tions in M^, which naturally support vorticity-type evolution (as a necessary feature of this 
basic fluid model), can develop blow-up rotational angular phenomena in finite time at a 
fixed stagnation point of the flow (Section [5]), which naturally leads to periodic cu-limits. 
It is not still clear if this twistor construction may lead to a truly localized blow-up swirl- 
like singularity. By the partial invariance of the twistor (it belongs to the 2D subspace 
(I1.23P ). it has a velocity field that is unbounded in the ^-direction. In our discussion, we 
show that the branching phenomena that could lead to 2;-localization of such a perturbed 
twistor, imply a very difficult evolution matching-like problem, which seems cannot be 
tackled rigorously still. In Section [HI in the spherical geometry, we involve extra precession 
axis and vertex mechanisms that make the problem more complicated but built a bridge 
to more realistic generalized quasi-periodic (or periodic) blow-up twistor behaviour. 

We do not rule out existence of other types of rotational singularities for (II. ip . especially 
in view of the "multiplicity curse" of blow-up scenarios; see (ll.29p below. The mathematics 



See bibliographical comments and references in Section |7| and on further extensions in Appendix C. 
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of such singularity structures promises to be unbelievably difficult since assumes sorting 
out an infinite (at least a countable) number of various possibilities. This "countability 
curse" for blow-up asymptotics will be explained later. In view of that, possibly, it is 
not an exaggeration to say that it would be much easier and much more pleasant to 
settle (11.21) . and hence would forget about this "awkward" multiple blow-up stuff. This 
is a warning again: in reality, blow-up theory is well-known to be very difficult even for 
looking very simple PDE models! However, for the NSEs, in author's opinion. 



1.28) both claims (II. 2p and (11.31) assume equally difficult proofs by blow-up scaling. 



In other words, to find the right answer, one needs to pass through a sequence of very 
difficult steps of analysis that in most issues are overlapping in both approaches. 

The following possible important feature of the blow-up twistors deserves mentioning: 
in view of their extreme rotational nature, they can create vorticity and sometimes velocity 
fields that tend to zero as t ^ T~ in the standard weak (integral) sense in L^^^. Therefore, 
though we observe that the vorticity gets infinite as t ^ T~, the local total mass (the 
integral over a small shrinking neighbourhood of the stagnation point) becomes negligible 
and just disappears at the blow-up time. 

Remark: On infinite family of patterns with regional and global blow-up. It 

turns out ^70j that there exists an infinite countable family of blow-up patterns on W2 in 
the cylindrical coordinates. As a compensation for a lack of proper mathematics here, in 
[To] , we managed to justify much simpler and rigorously conffimed construction of blow- 
up space jets that exhibit effective regional or global blow-up in the radial variable. This 
example underlines another important feature of blow-up for (11.11) (though these always 
are of infinite energy): 



1.29) there can exist an infinite countable family of blow-up patterns. 



In particular, there exists effective regional blow-up of the z-component of the vector field 
in the radial r-direction, i.e., the blow-up wave does not propagate and is of a standing 
type. This is a good sign, but of course the velocity field is unbounded in the z-direction, 
so that the patterns have infinite energy. A proper "bending" of the 2;-axis to create a 
kind of "blow-up ring" is rather suspicious, since also will exhibit infinite kinetic energy 
as t ^ T~ . A more complicated geometry is necessary for the next refined blow-up 
construction. Another important feature is 



1.30) these blow-up patterns converge to similarity solutions of Euler's equations. 



In the rescaled sense, the convergence turns out to be uniform in M^! We do not know 
whether the phenomenon (11.301) is expected to be generic for other hypothetical blow-up 
patterns of finite energy. 

Similarly, in [70], for the 3D Euler's equation^ (EEs), 

(1.31) + (u ■ V)u = -Vp, divu = in x M+, 



^Leonhard Paul Eulcr, 1707-1783. 
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the restriction to W2 is shown to admit patterns with single point blow-up in the r-variable. 
•Blow-up patterns with swirl and convergence to EEs fll.3ip (Section\B\j. This is about 



possible extensions of the resuhs in [70] to the non-radial geometry with the blow-up 
swirl included, which is a training for more complicated blow-up structures to come. 

• Blow-up about Slezkin-Landau and other equilibria (Section^. This can lead to blow- 
up patterns by a kind of "centre-stable manifold" analysis by using various singular or 
regular stationary solutions of the NSEs. Since some results are derived on the basis of the 
famous Slezkin-Landau singular steady states, we include some history of these solutions 
and put translations of Slezkin's rare notes of 1934 and 1954 in Appendices A and B at 
the end of the paper to be followed by a further discussion in Appendix C. 

As a rather general conclusion, we state the following "steady-TW exercise" for the 
rescaled NSEs, which should be solved before even talking about (11.21) or fll.3p : 



;i.32) 



to describe all families of equilibria in M'^ \ {0} or M^, 
crv^ + |z/-Vv + iv-F(v- V)v = - Vp + Av, div v = 0, 



where a G M is the nonlinear eigenvalue introduced as in (11.251) being the blow-up angular 
velocity. For cr = 0, (ll.32p is Leray's classic rescaled stationary problem with a number 
of nonexistence results, jl62[ l2U8"t I152j (though not everything is still known, [2U4]). 

For cr 7^ 0, suitable profiles in (ll.32p represent the simplest case of non self-similar blow- 
up with the omega-limit set consisting of a periodic orbit. Nevertheless, plugging into the 
system a single new real eigenvalue a can be non-sufficient for settling the question on 
possible non-trivial cu-limits of rescaled blow-up orbits, so we will need to introduce further 
mechanisms of axis precessions leading to quasi-periodic and other motions. 

In general, even for a = 0, the stationary (nonlocal elliptic) problem: 



1.33) 



describe all singular equilibria in \ {0}, 
(U ■ V)U = - VP + AU, div U = 0, 



remains still open (though some key steps have been already made, [2081 12041 11171 IllSj . 
which inspire some optimism that the situation is under a proper control; see below). 
As a clue exercise to such a complexity, in Section [TJ we demonstrate that complicated 
oscillatory sign changing singular equilibria occur for the simplest elliptic problem from 
R-D theory in the supercritical Sobolev range: 

(1.34) Am + \u\P"^u = in \ {0}, where p > ps = > 3. 

Non-radial singularity structures for (11.341) were not addressed in the literature and are 
unknown. These can generate extremely complicated blow-up patterns in the correspond- 
ing parabolic R-D flow; see (11.630 below. 

In fact, the proposed swirl behaviour (I1.25P naturally corresponds to the rotational 
invariance of some symmetric singular equilibria in (11.330 . In general, for possible more 
complicated singular stationary profiles U(x), new types of "rotations" are necessary to 
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introduce and understand. One of a possible scenario of blow-up in the NSEs is as follows: 



(1 35) orbit {u(t)} evolves as t — > T "close" to some singular stationary 

^ ' ^ manifolds fll.33p being "trapped" in their complicated structure. 

Then "swirling type" of the orbit behaviour close to x = as t — on shrinking compact 
subsets will naturally and entirely depend on the singular steady manifold involved and 
can be extremely complicated leading to any of multi-dimensional quasi-periodic or even 
chaotic attractors. 

• Twistors in spherical geometry ( Section We next perform a formal study of a blow- 
up twistor in the spherical coordinates that are necessary for creating a truly spatially 
localized pattern. We postpone our final conclusions inherited from the previous analysis 
until Section [9l As a by-product, we again naturally arrive at the problem fll.33p . which 
is difficult and open, regardless a good progress made recently on understanding of the 
scaling nature and uniqueness of famous Slezkin-Landau singular solutions of a submerged 
jet, |208l I204t 1155] : see Section [5] for details. Thus, we again discuss the possibility 
that finite kinetic energy blow-up twistors may be generated in a small shrinking as 
t — > T~ vicinity of every "proper" singular steady states from (11.331) . which possibly is 
not precisely of homogenuity — 1 in r = |x| and possessing torsion-precession mechanism 
of its swirl axis. Other, more regular "steady states" are also of importance. 

Thus, on the basis of the given fundamental model (11. ip . we will show some mathe- 
matical tools that are necessary to tackle general difficult problems of blow-up and non 
blow-up (then a local smooth solution becomes global and hence unique) for complicated 
higher-order PDEs. 

1.6. On some reaction-diffusion and parabolic analogies to be applied: first 
exercise on Type I and Type II blow-up patterns. Nowadays, blow-up PDE 
theory, as a self-contained subject, embraces a wide range of various nonlinear evolu- 
tion models; we list a few monographs from the 1980s and later periods up to 2007, 
m [m EHl ESI [US [M UMl HSZl M , where further extensions and references can be 
found. Most of them are mainly and specially devoted to blow-up behaviour in nonlin- 
ear partial differential equations of parabolic and hyperbolic types (those two PDE areas 
are well-established for blow-up since the 1950s) and contain key literature and various 
blow-up results achieved in the last fifty years. During this long time, a huge amount 
of beautiful and difficult conclusions were obtained resulting in deep and sometimes ex- 
hausting understanding and complete classification of blow-up patterns for some nonlinear 
parabolic and other models. 

Therefore, for those who are constantly working in these PDE areas, it is not surprising 
that the crucial two the so-called Millennium Problems such as 

(I) Global existence or nonexistence of smooth solutions of the NSEs in M^, and 

(II) The Poincare Conjecture (see further comments below), 

both, in their already existing (for (II)) and possible (for (I)) ways of the solution, heavily 
rely on the study of blow-up solutions of the nonlinear evolution equations involved. 
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Frank-Kamenetskii equation: self-similar (Type I) and Type II fast blow-up patterns. In 
discussing further consequences related to both problems ( 11. 2p and ( 11. 3p . we apply some 
general results of blow-up theory developed for various reaction-diffusion equations and 
systems, which turn out to be fruitful in application to blow-up for the more complicated 
PDEs such as fll.ll) . In particular, it is well-known that even for simple looking semilinear 
parabolic reaction-diffusion PDEs such as the classic Frank-Kamenetskii equation (1938) 
[60] developed in combustion theory of solid fuels (also called the solid fuel model; first 
blow-up results in related ODE models are due to Todes, 1933) 

(1.36) Mt = AM + e" in x R+, 

the property (11.291) holds. Namely, the first family of "linearized" blow-up patterns is con- 
structed by linearization techniques and further matching of centre and stable manifolds 
orbits, where the latter one is infinite-dimensional causing the eventual countability of 
the family; see explanations below. Moreover, it is crucial that this family of blow-up 
asymptotics can exhaust all possible types of blow-up behaviour that is available in the 
model. Such delicate results for = 1 or 2 and for other R-D PDEs (11.631) are due 
to Velazquez |213] . Hence, the family of blow-up patterns for (I1.36P from |213i I214j is 
evolutionary complete (a notion introduced in [HS], where further references can be found). 

For > 3, (11.361) possesses non-trivial self-similar blow-up (Type ]@) patterns, which 
are called nonlinear eigenf unctions, and, in addition, the set of linearized patterns is more 
involved and can include other countable families. The total family of blow-up patterns 
gets more complicated, so its evolution completeness is unknown representing a difficult 
open problem for (11.361) . 

It is curious that that this remains open even for the radial equation (11.361) in M.^, 

(1.37) Ut = Urr + ^Ur + e'' m ]R+ X M+ {u = u{r,t), r = \x\ > 0, Ur{0,t) = 0). 

To illustrate this fact and our future arguments in Section I5.12[ we briefly explain how a 
countable family of non self-similar patterns of Type II can occur. 

Similarity blow-up patterns. Performing in (I1.37P the standard self-similar blow-up scaling 
yields the rescaled equation, which are known to exist for > 3 [TTj, 

uir,t) = -\n{T-t)+viy,T), y = T = -\n{T-t) 

(1.38) 

=^ Vr = H(f ) = Anv - i^yvy + - 1. 

The self-similar rescaled profiles f{y) ^ are its good equilibria (we omit details since 
will be talking a lot about similarity blow-up later on): 

(1.39) H(/) =0 in R^, f{y) has at most logarithmic growth as y ^ oo. 

According to (ll.38p . each /(?/) defines a Type I self-similar blow-up (see below). 

Non-similarity blow-up patterns. Equation (11.391) in dimension > 3 admits a singular 
equilibrium of the form 

(1.40) v(2/) = in^ in M^\{0}. 



'The terms "Type I, II" were borrowed from Hamilton [M], where Type II is also called slow blow-up. 
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Then some blow-up patterns can evolve as r — > +00 {t —>■ T^) close to this "stationary 
manifold" . 

Inner Region I expansion: linearization. To see the applicability of this idea, we perform 
the linearization in (11.381) . 

(1.41) v{T) = V + Y{r) =^ Yr = ANY-lyYy+'-^Y + B{Y), 

where D is a quadratic perturbation as F ^ in a suitable metric. By classic theory of 
self-adjoint operators [15], the linear operator with the inverse square potential in (11.411) 

(1.42) H'(y)=A^-il/-V + ^/ 

is well posed in if^, (M^) (with the weight as in (I2.24p ). i.e., has a compact resolvent 
and a real discrete spectrum, provided that the inverse square potential is not too much 
singular at the origin y = 0. Namely, one needs that 2(A^ — 2) is less than the constant 
ch of embedding Hq{Bi) C L^^|_i(i?i) by the Hardy classic inequalit}0: 

(1.43) /^^g^dx<CH/5jVM|2dx =^ 2(iV-2)<CH= (^)', i.e., iV > 10. 

Thus, in integer dimensions > 11 (A^ = 10 has own peculiarities; see [73]), the 
operator (ll.42p has a discrete spectrum and radial eigenf unctions satisfying 

a(H'(V)) = {Afc, A; = 0,2,4,...}; rkiv) - hy" + ■ ■ ■ , y ^ +00, 

i^kiy) -Cky ^ + ... , y-^ 0, where 6 = ^\ - 

and bk and Ck are some normalization constants. The orthonormal set of eigenfunctions 
$* = {i'l} is then complete and closed in L^*(M^) in the radial setting. 

Hence, for > 11, equation (11.371) admits very special asymptotic patterns; see [I71I72] 
(global solutions and a survey) and |98li59] (blow-up solutions and a survey). In particular, 
concerning fast blow-up patterns of Type if^, these are obtained by matching of the 
linearized behaviour on the steady singular manifold for any A; = 2, 4, ... such that Xk < 0: 

(1.45) Vk{y,t) = V{y)-Cke'''^rk{y) + --Fk{y,T) = -2\ny-e^>=^y'' as y ^ 0, 

with a bounded flow at the origin r = 0; see [HHl ESI HZ] for first results in this direction, 
and further references above. The correct L°°-rate of blow-up behaviour is obtained by 
calculating the absolute maximum of Fk{y, r): for r ^ 1, 

(1.46) {Fk)'y = at ~ e— =^ supy Fk{y, r) = Fk{yk, r) ~ a^r, = > 1. 

Inner Region II: quasi-stationary regular flow. Thus, in Region II close enough to r = 0, 
one needs to solve the original equation (I1.37P with the condition as in (ll.46p . which 



"'^''it's idea goes back to the 1920s, [95], and in this form was aheady used by Leray in 1934 [137) : see 
Section [7?3l for extra details. 

^"'^In blow-up R-D theory, Type II assumes faster non-self-similar growth. 
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suggests the following scaling therein: 

m(0, t) = — ctfc ln(T — t) ^> u{r,t) = —akln{T — t) + w{^, s), 
(1-47) ^here ^ = (^z^ and , = -L_ (r - ^ +00. 

Substituting this into the F-K equation, after elementary manipulations yields the fol- 
lowing perturbed problem for w: 

(1.48) Ws = A^w + 6"" - ^ 7(1 + 1 * > 1- 

Recalling that according to fll.461) . we have to have that w{0,t) = 0, general stability 
theory for such blow-up singularity problems [79] suggests, that since (11.481) is a per- 
turbed gradient flow, there is the stabilization to the unique bounded stationary solution: 
uniformly on compact subsets in as s ^ +00, 

(1.49) w(e,s) ^ Vr(0, where A^VT + e'^ = 0, 1^(0) = (W^(0 ~ -2 In^, ^ > !)• 
Eventually, in this Region II, the current blow-up pattern behaves as: 

(1.50) u{r,t)^ak\HT-t)\ + W{j^-^) as t ^ T". 

Thus, (11.381) gives the asymptotic structure of such blow-up patterns {uk{r,t), k > 2} 
with the following fast blow-up rate: 

(1.51) Mfc(0,t)~ (l + ^)|ln(T-t)| as t ^ T~ (A^ ~ -|, A;>l). 

These Type II non self-similar blow-up patterns can have arbitrarily fast growth for /c ^ 1 
than the standard Type I similarity divergence associated with the ODE 

(1.52) M' = e" =^ u{t) = \\n{T -t)\ as t ^ . 

For = 10, the origin ?/ = is in the limit-point case for the operator (ll.42p . so that 
its deficiency indices are (2, 2) jl60] . and there still exists its proper self-adjoint extension 
with a discrete spectrum and compact resolvent, so that (ll.44p holds; see details of an 
application in [73] . This also allows to construct blow-up patterns with a slightly different 
blow-up rates. 

On the other hand, for G [3, 9], instead of (ll.44p . the origin is oscillatory: as y — 0, 



(1.53) V^*(i/) ~y-^[CiCos(61ny) + C2sin(61ny)], where h = ^Kio jv) ^ 

This yields that any eigenfunction of any self-adjoint extension of (11.420 (existing by 
classic theory |160j ) will make the function (ll.4ip to be sign-changing and oscillatory near 
the matching origin y = 0. Indeed, for nonnegative solutions, this prohibits any matching 
at all, so that Type II blow-up patterns are then nonexistent, which is rigorously proved 
in Matano-Merle |146j . Therefore, it is plausible that for radial nonnegative solutions in 
dimension 3 < < 9, the only blow-up patterns are exhausted by the family of linearized 
ones (as in Section [5.120 and the nonlinear self-similar solutions (ll.39p . though there is no 
proof of such a completeness. However, in the nonradial geometry for solutions of changing 
sign such a classification blow-up conclusion is far away from being even formally justified. 
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In view of (11 ■4p . one of our tricks is to check weather Type II non-similarity blow-up 
patterns can be constructed for the NSEs, though this should be much harder of course. 

Quasilinear reaction- diffusion extension: countable sets of blow-up patterns. For a quasi- 
linear extension of fll.361) in ID, with the so-called p-Laplacian operator (here p \—>- a), 

(1.54) ut = i\u,\''u,)^ + e^ (a>0), 

a countable set of blow-up patterns was constructed in [21]. As an intrinsic feature of 
essentially quasilinear problems, it was shown that, depending on cr > 0, first patterns 
represent nonlinear eigenfunctions, i.e., are self-similar, while the rest are constructed by 
the linearization techniques and matching as for (ll.36p . For all a > a^o = 0.60..., all 
blow-up patterns are nonlinear, i.e., then (ll.54p admits countable families of self-similar 
solutions. Moreover, it was also shown that, in view of the strong degeneracy of the 
p-Laplacian operator in (ll.54p . there exist other blow-up similarity solutions, which can 
compose uncountable families. It is still an open problem to prove evolution completeness 
of these blow-up families. We reflect this discussion as follows: in the family (ll.29p . 



1.55) there may exist infinitely many nonlinear or "linearized" blow-up patterns. 



Higher-order parabolic extensions. Some types of blow-up singularities are already known 
for higher-order react ion- diffusion equations such as 

(1.56) ut = -i-A)"'u+\u\P-^u in x M+ (m > 2, p > 1) 

(1.57) and Ut = -A'^u - A{\u\p-'^u) in x M+ {p > 1); 

see [22l \56\ [55| [M] and references therein concerning other models with blow-up. The 
latter limit unstable Cahn-Hilliard equation (I1.57P has the divergence form and hence the 
flow preserves the total mass of L^-solutions (in this sense this mimics a similar feature of 
(11. ip ). Equation (11.570 is shown to obey Leray's blow-up scenario (II. 5p [70] . Note that 
the questions on a complete description of nonlinear and linearized blow-up patterns for 
(11.561) or (I1.57P and evolution completeness of the whole countable family remain open 
even in the one-dimensional case = 1 and m = 2 (the fourth-order diffusion only). 

1.7. Second discussion around the NSEs. Thus, the negative conclusion (ll.4p rules 
out a standard self-similar way of blow-up in the Navier-Stokes equations. Anyway, 
regardless known deep ideas on interaction of "vortices tubes" and others (mainly coming 
from fluid dynamics and less mathematically developed), it seems we cannot still imagine 
how complicated other non self-similar individual "linearized" blow-up patterns can be 
(for instance, these can be partially invariant, i.e., "partially nonlinear", or non-invariant 
at all), to say nothing about possible interaction of a number of different related blow- 
up "vortex tubes". This is our aim to give, on the basis of a reaction-diffusion blow-up 
experience, a possible new insight into the world of imaginary singularities for (II. ip . 

In view of a more complicated nature of the Navier-Stokes equations that are a system 
of four equations in the four- dimensional space-time continuunj^ R'^ x R_|_, the slogan 



'I.e., a "dynamical system 4 x 4 = (4 dependent x 4 independent variables)". 

16 



(11.291) should contain more surprises. At least, one can expect a formally infinite number 
of possible types of blow-up to be checked out, and, to achieve the positive answer (11. 2p . 
all of them should be ruled out by the assumption of finite energy and/or others. It is not 
an exaggeration to assume that (II. ip may admit a countable set of different countable (or 
even uncountable) scenarios of formation of single point blow-up patterns that admit local 
construction (on shrinking subset to the blow-up point) by some kind of combination of 
linearized and nonlinear eigenvalue techniques. Then the crucial fact on the finiteness of 
the energy of the globally extended patterns can be checked only after matching (or non- 
matching) of these local singularity structures with surrounding bundles of less blowing up 
or even bounded "tails" . Such matching procedures (which are responsible for existence 
or nonexistence of the patterns) are also expected to be extremely difficult. 

Vice versa, to achieve the negative claim (II. 3p . in order to construct a suitable blow-up 
pattern, which prohibits global evolution of bounded smooth finite kinetic energy solu- 
tions, one then needs to sort out, possibly again, a huge (infinite and even uncountable) 
number of a priori unknown blow-up solution structures. 

Note that, for many other models including various nonlinear higher-order PDEs, for 
which there is no any hope to control existence and/or nonexistence of blow-up singu- 
larities by some conservation/monotonicity/order-preserving/etc, mechanisms (that are 
nonexistent for sure), the careful study of (ll.2p and (II. 3p is unavoidable, is very difficult 
indeed, and, in many cases, in a full generality, is not doable at all. The latter means that 
this problem is not analytically solvable, i.e., constructive conditions that guarantee (11.20 
(or (ll.3p ) cannot be derived. Nevertheless, this is not a manifestation of any kind of a 
"PDE agnosticism" , since the most stable and generic asymptotic structures can and must 
be studied and understood by any, rigorous or not, mathematical means. The illusive is 
a full description of all the possible singularity patterns and their evolution completeness 
(see a proper setting for this below) to eventually guarantee (11.20 or (11.30 . 

1.8. On a related blow-up parabolic area: Ricci flows and the Poincare Conjec- 
ture. Another amazing geometry and PDE area, which is not less famous nowadays, is 
the study of blow-up structure of Ricci flows of metrics in M^. Following the logic of evo- 
lution blow-up PDE analysis, it can be characterized that Perel'man's new monotonicity 
formulae and principals of his blow-up surgery at finite-time singularities of the Ricci fiow 
made it possible to guarantee to get a proper global extension beyond all blow-ups and 
with a kind of necessary "symmetrization" at the finite extinction time. In particular, 
this allowed to prove the Poincare Conjectur^ see [21j for details, history, references, 
and recent development. Proposed for this kind of analysis by Richard Hamilton in 1982 
the Ricci flow for a family of Riemannian metrics g{t) = {gij{t)}: 

(1.58) g{x,t) : gt = — 2Ric(7, or, for components, {gij)t = —'2Rij, 

where Ricg = {Rij} is the Ricci curvature of g, represents a system of parabolic PDEs 
for the components, which overall obeys the Maximum Principle (the MP) and other 

^■^"^4 closed, smooth, and simply connected 3-manifold is homeomorphic to S"^. This still remains the 
only solved Millennium Problem among other seven. 
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related classic properties of parabolic flows. The scalar curvature equation for the scalar 
curvature R = {g'^^ Rij){x,t) is a semilinear parabolic PDE with a quadratic nonlinearity, 

(1.59) Rt = AR + 2\Ricg\'^, 

where A < is the corresponding Laplacian. Moreover, in dimension 2, the scalar 
curvature R{x, t) then takes the form of a standard semilinear quadratic R-D equation 
[SI p. 227] 

(1.60) Rt = AR + R^ {N = 2). 

This shows a (obviously, well-known) link between Ricci flows singularities with standard 
and well-developed blow-up R-D theory. Similar to (11.361) . blow-up for this equation 
has a long history and embraces hundreds of publications in the 1980-90s in almost all 
leading world journals on nonlinear PDEs and a few monographs mentioned above. For 
future use, let us note that a key idea of its generic non-trivial blow-up behaviour with 
logarithmic (the appearance of ln(T — t)-term is due to a centre subspace behaviour as 
explained in Section [5.121 so we avoid further comments), 

(1.61) ~ V|ln(r-t)| as t^T-, 

deformation of similarity blow-up structures dates back to Hocking, Stuartson, and Stuart 
in 1972, [100] . It was proved rigorously twenty years later that the actual structurally 
stable blow-up behaviour for (11.601) as t — > T~ close to the blow-up point x = is given 
by 

(1.62) R{x,t) = ^j4-^^=^=^){l+o{l)), where /*(?/) = tt^, 

and c* > is a constant depending on the dimension only (c* = ^ for = l); see 
|213j and \IS7\ p. 312] for references and results. A full account of results on blow-up 
for semilinear equations such as (ll.60p can be found in the most recent monograph by 
Quittner and Suplet [183j . Formula (ll.62p is universal for dimensions A^ < 6, i.e., in the 
subcritical Sobolev range for the elliptic operator in (11.601) . with the general term, 

(1.63) Ut = Au+\u\P-\ with p = 2 < psoboiev = ^ {N>2). 

As an accompanying key feature for further use here, the structurally stable (generic) 
blow-up pattern (11.621) perfectly serves as a powerful confirmation of the slogan (ll.27p . 
Indeed, the similarity scalings for the equation (ll.60p are standard and simple, 

(1.64) i?(a;,t) = _L_^, y=^v 

so do not contain any of logarithmic factors as in (ll.6ip and (11.620 . which are created in 
the blow-up evolution as t — >■ T~ for almost arbitrary suitable solutions with flL^-data. 

It is indeed also worth mentioning that since the function f^:{y) in (11.620 is radial, this 
expresses the phenomenon of strong symmetrization of blow-up structures near blow-up 
time that is true for almost all admissible non-constant solutions of (ll.60p (classification of 
singularities is also a key feature of Perel'man's proof, which the crucial blow-up surgery 
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is based upon, [3T], Ch. 7]). We will constantly use the generic blow-up behaviour (11.621) 
as a source for further speculations concerning the model (II. ip . 

It is remarkable that one of the key ideas of Perel'man's proof was surgery followed by 
a necessary classification of possible blow-up singularities, which can occur for the Ricci 
flow ffL58|) . Recall our main target for the NSEs. 

1.9. Reaction-diffusion: a new Type II blow-up patterns for p = ps- The par- 
abolic equation (11.631) in the critical case p = Ps for = 3 reveals another type of 
construction of Type II non self-similar blow-up patterns. Namely, in the critical Sobolev 
case, (11.631) admits the Loewner-Nirenberg stationary solution (1974) [142] 

This is indeed a very special and remarkable explicit solution, so that, similar to the sce- 
nario in Section [L6l blow-up can occur about the ID manifold of such rescaled equilibria. 
Namely, following [59j, we explain how this can happen. First, as in (11.381) . we perform 
the standard self-similar scaling in (ll.63p . with p = 5 for = 3, 

u{x,t) = {T-t)-^^v{y,T), y=^^, T = -\n{T-t), so 

(1.66) ^ 

Vr = H(t') = Av ~ \ y ■ Vt" — ^;3Y V + |f 

Second, let us assume that f (y, r) behaves for r ^ 1 being close to the stationary manifold 
composed of equilibria (11.650 . i.e., for some unknown function ifij) — > -|-oo as r — > +oo, 

(1.67) v{y, t) ^ v5(r)MLN (v^^ {T)y) 

p — 1 

on the corresponding shrinking compact subsets in the new variable C, = ip^~{T)y. It then 
follows that, in the case = 3, on the solutions (11.650 in terms of the original rescaled 
variable y (see computations in [591, p. 2963]; our notations have been slightly changed) 

(1.68) \y(y^r)\P-^v{y,r)-^^6{y) as r ^ +oo 

in the sense of distributions. Therefore, equation (11.660 takes asymptotically the form 

(1.69) v^ = U{v)=Av-ly-\/v-\v + ^5{y) + ... for r > 1. 
Thus, we get an eigenvalue problem for Hermite's classic operator fTEl p. 48]: 

(1.70) B* = Av-lyVv, where a(B*) = {A^ = = 0, 1, 2, ...}, 

defined in ^^.(M^), p*{y) = e-l*'!"/^, with the domain H^,{M.^), to be used lateS It 
follows from (I1.69P that one can try the following regular parts of such Type II blow-up 
patterns: balancing two terms in ( 11.690 yields 

(1.71) Vf,{y,T)=e^^^--^>Hf,{y) + ... =^ V5^(r)~e^^ for r > 1, 

^^The same operator occurs for the rescaled equation (|1.60p . where the factor ^/\]Il(J^~^t)\ — ^/r is 
due to the eigenspace behaviour with A2 = — 1 and the second Hermite polynomial H2{y) = (1 — ^). 
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1.65) 



N{N-2) 
N{N-2)+\x\^ 



where Hpiy) are Hermite poljTiomials as the eigenfunctions of (11.701) . Together with 
the scahng in fll.66p . this yields a countable family of complicated blow-up structures. 
To reveal the actual space-time and changing sign structures of such Type II patterns, 
special matching procedures apply. In j59j , this analysis has been performed in the radial 
geometry, though (and this is also key for us to recognize) still no rigorous justification 
of the existence of such blow-up scenarios is available. 

1.10. Further comments on the "RD-sense" of this essay. As was already an- 
nounced, the present paper is a general view from blow-up reaction- diffusion theory to 
possible "singularities" (including, e.g., multiple zeros) of the 3D NSEs. In other words, 
the main intention of the author is to involve some RD-experience into the NSEs study. 
As is easily detectable by experts in the area, this special issue dictated, in particular, a 
quite special kind of survey and references involved in the list, which by no means reflects 
the actual history and modern trends of NSEs and Euler's equations theory (to say noth- 
ing about R-D theory as well). Nevertheless, the author hopes that a "detached observer 
view" can deliver some new useful accents to the area. 

The subjects of the R-D and the NSEs do not coincide at all and seems their overlapping 
is suspicious in many places. Since the NSEs subject is huge and extremely difficult, the 
author apologizes for any inconvenience caused by his approximate constructions and/or 
speculations that, from classic fluid dynamic PDE views, can be classifled as rather non- 
consistent or well known somehow. It seems understandable that, in view of not that 
optimistic final conclusions (cf. (19.41) ). which were predicted in advance, the author, with 
a completely different background, did not have a very strong motivation to make more 
clear some of the matching constructions when this looked being possible (not often, 
unfortunately). The "semi- mathematical" (and sometimes logically incomplete) language 
of presenting the speculations reflects the obvious truth: each elementary step in this long 
story of singularities/blow- up for the NSEs could take years to flx (and never known how 
many). So, the author took the risk to show quite a discontinuous route (a "jumping 
frog" -style) to the end of this story, of course, with a full understanding how risky this 
way could be in mathematics, when many steps and concepts of his speculations can be 
attacked by justifled incinerating critics from some experts from the fleld. Anyway, the 
author hopes that the attentive Reader will flnd some ideas and concepts of blow-up from 
the R-D systems useful, even if these were not presented on a sufficiently rigorous basis 
and costs; and even will find satisfactory a rigorous ("almost", i.e., if fixed in a reasonable 
finite time) construction of a countable family of blow-up space jets for (11.11) . 

Indeed, the combination of all the possible tools of singularity analysis, leading to a 
success for the NSEs (11.11) . will serve as a solid and reliable basis for further development 
of modern theory of higher-order PDEs. In this context, blow-up scaling methods become 
more and more penetrating into the core of PDE theory being natural and unavoidable 
tools of evolution analysis of higher-order nonlinear equations in all its three "hipostases" : 

(i) existence, (ii) uniqueness, and (iii) asymptotic behaviour. 
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2. Some facts of singularity history for the Navier-Stokes equations: 
Leray's blow-up solutions, singular points, spectra, and patterns 

We present here a short survey on global solvability, singularity formation, and other 
classic facts concerning the Navier-Stokes equations (11. II) . and stress a special attention to 
some blow-up issues, which will be used in what follows. A perfect and detailed overview 
of main mathematical results concerning the NSEs is available in Taylor |206[ Ch. 17], 
which includes several aspects to be quoted below without proper referencing. 

2.1. Leray— Hopf ("turbulent") solutions of finite kinetic energy. It is a classic 
matter that the energy L^-norm is natural for (11. ip . After multiplication by u in the metric 
(■, ■) of L^(]R^) and integration by parts, the convective and pressure terms vanish on 
smooth enough functions u(x, t) with sufficiently fast (say, exponential) decay at infinity, 

(2.1) ((u-V)u,u)=0 and - ( Vp, u) = (p, V ■ u) = 0. 

Therefore, on such smooth solutions, we have the instantaneous rate of dissipation of 
kinetic energy given by 

(2.2) liMtm = -\\Dn{t)\\l =^ ||u(t)|h + 2/pu(.)||id. = ||uo|b, t > 0. 



Actually, the estimate in ( 12. 2p with the inequality sign " < " is the energy inequality 
for Leray-Hopf weak solutions of (11. ip (in 1933, Leray also called such weak solutions 
"turbulent" |1371 p. 231, 241] and compared these with regular ones; see also Lions |14m 
Ch. 1, § 6] for a discussion); q.v. e.g., |182] and references therein. (12. 2p is also a crucial 
identity for general turbulence theory. E.g., the famous Kolmogorov-Obukhov power "K- 
41" law (1941) [119[ I165j for the energy spectrum of turbulent fiuctuations for wave 
numbers k from the so-called inertial range, 

(2.3) E{k) = Celk-i, where e = {\Dm\^) , 

uses this rate of dissipation of kinetic energy for high Reynolds numbers Re ((■) is an 
invariant measure of calculating expected values). 

Note another, weaker "conservation laws" , since the convective terms in the NSEs are 
in divergent form (® is the tensor product of vectors in R^): 

(2.4) (u • V)u = div (u (g) u) = {uxi)^^ + {vu)^^ + {wu)^.^. 

Hence, integrating over with the necessary decay at infinity, one can get that the total 
"masses" of the velocity components are preserved: 

^/M(x,t)dx, ^/t;(x,t)dx, ^ / w(x, t) dx = 0, i.e., 

(2.5) 

/ u{t) = Juo, J v{t) = Jvq, J w{t) = Jwo for t> 0. 

Of course, as estimates, (12.50 are weaker than (12. 2p . and the mass "semi-norms" by taking 
the absolute values of the integrals there are difficult to apply for solutions of changing 
sign. Though, as sharp conservation properties, these can be valuable in construction (or 
prohibiting) special sensitive blow-up patterns. Moreover, for some parabolic problems 
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with the mass conservation (and no L^-control), estimates such as fl2.5p can play a key 
role for global extension of blow-up solutions beyond blow-up, for t > T. A discussion 
on the application of Leray's blow-up scenario (11. 5p to the limit Cahn-Hilliard equation 
(11.571) is presented in [70] . 

Nevertheless, the L^-bounds as in (12. 2p are not sufficient to control the L°°-non-blowing 
up property of solutions, and this is the origin of extensive mathematical research in the 
last seventy years or so. Recall that global weak solutions of (II. ip satisfying 

u e L°°(M+; L2(M^)) n L2(M+; H\R^)) 

were already constructed by Leray jl36[ I137j (1933), and Hopf |101j (1951). We recom- 
mend recent papers P, [161 ISll [371 US] as a guide to various results of local and global (for 
small and other data) theory of (II. ip . including analyticity results in both spatial and 
temporal variables; see [IHl 1224] for a modern overview. 

2.2. Blow-up self-similar singularities with finite energy are nonexistent: on 
Leray's scenario of backward and forward blow-up self-similar singularities. It 

seems that the original idea that the classic fundamental problem of the unique solvabil- 
ity of (II. ip in M^, i.e., existence of a global smooth bounded L^-solution, is associated 
with existence or nonexistence of certain blow-up singularities as t T~, goes back to 
Th. von Karman; see [T09l[TT0l[TTT]R and survey [13]. 

As usual, similarity solutions, as a next manner to further specify the behaviour, can 
be attributed to an invariant group of scaling transformations: if {u{x,t), p{x,t)} is a 
solution of the NSEs ([HI]), then 

(2.6) {ux{x,t) = Au(Ax, A^t), px{x,t) = \'^p(\x, A^t)} is also a solution for any A G M. 
Setting here A = formally yields two types of self-similar solutions of (II. ip . the 



blow-up and the global ones. Let us now focus on key historical aspects of such a study. 

Namely, in 1933, J. Leray [1361 1137] proposed a mathematical question to look for 
blow-up in (II. ip driven by the self-similar solution^ 

(2.7) u(x,t) = ^U(y), pix,t) = ^^Piy), where 1/ = 7=. 

This Leray's statement is well known and was stressed to in many papers [T371 p. 2250 



"... la solution des equations de Navier dont il s'agit est: 

(3.12) Ui{x,t) = [2a{T -t)]-^W,[{2a{T -t))-^ x] {t < T) 

(Xx designe le point de coordonnees Axi , Xx2, Aa;3 . ) " 



^^The author apologies for not being able to trace out von Karman's original work (or a lecture?), 
where the idea of singularity of the velocity field ^ {T — t)~" (a = |, |, or |?) appeared first. 

"'^^Such a blow-up backward continuation variable y for ID linear parabolic equations was already 
systematically used by Strum in 1836 [203j : on his backward-forward continuation analysis, see [66l P- 4]. 

^'^Here and later, boxing and underlying are author's. 
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However, at the end of the same paper in Acta mathematica, Leray returned once 
more to this similarity blow-up problem and now his question is also truly remarkable 
|137l p. 245] (here (3.11) is the system (12. 9p below for the similarity profiles in (12.71) ): 

"Remarque: Si le systeme (3.11) possede une solution non nulle Ui{x) 
cette solution permet de construire un exemple tres simple de solution 
turbulente c'est le vecteur Ui{x,t) egal a 

(2.8) [2a{T -t)]-^ U,[[2a{T -t)]-^ x] pour t < T et a o pour t > T; 

11 existe une seule epoque d' irregularite : T." 

Therefore, Leray in (12.81) posed both principal questions on existence of self-similar 
solutions of blow-up backward type for t < T and of the standard forward type for t > T, 
which are naturally supposed to "coincide" (in which sense? - in general, a difficult 
question of extended semigroup theory) at the unique singularity point t = T. Thus, this 
is a principal setting not only for a similarity way of formation of a blow-up singularity 
as t ^ T~, but also for self-similar continuation of the solution for t > T, i.e., beyond 
blow-up time, when it becomes again regular and bounded. Even for simple parabolic 
reaction-diffusion equations such as (11.361) and (ll.60p . though such "peaking" blow-up self- 
similar solutions are known, a theory of such an incomplete blow-up is far away from being 
well-understood and contains a number of open problems; see references and results in 
[77] and in survey [78] (as usual, further recent progress already achieved in this direction 
can be then traced out via MathSciNet). For the Cahn-Hilliard equation (11.571) . this is 
done in [70] in the lines of Leray's scenario (11.51) . 

Thus, substituting (12.71) into (II. ip yields for U and P the following "stationary" system 
(3.11) in jm\ p. 225]: 

(2.9) iU + |(?/- V)U + (U- V)U = -VP + AU, divU = in R^. 

During last twelve years, a number of enhanced negative answers concerning existence 
of such non-trivial similarity patterns (12. 7p . (12. 9p were obtained. The key ingredient [162] 
of such nonexistence proofs is the Maximum Principl^: 

(2.10) n = i|Up + i?/-U + P satisfies - AH + (U + § y) ■ VH = -|curl Up < 0; 

see further details in [311 1152[ 1208] . and the advanced and negative nonstationary PDE 
answer in [103] . Let us note an existence result in [38] for N = 4. 

Regardless the nonexistence of the similarity blow-up (12.70 and especially in connection 
with, the following Leray conclusion deserves the attention [1371 P- 224]. For the function 

(2.11) V{t) =snp^ \u{x,t)\ for t < T : 

The fact that ^ |up +p satisfies the MP for the stationary NSEs is well-known; see e.g., [90] and 
more references in |162) . 
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"Un premier caractere des irregularites : Si une solution des 
equations de Navier devient irreguliere a I'epoque T, alors V{t) 
augmente indefiniment quand t tend vers T ; et plus precisement: 



(3.9) V{t) > A ' 



T-t 



Thus, Leray was the first who derived this estimate from below for a solution to have 
blow-up at t = T. This estimate led Leray to pose the question on blow-up similarity 
solutions (12. 7p (not just the fact that this would be a simple dimensional way for the 
NSEs to develop a finite-time singularity as t — > T~ with an analogous extension for 
t > T). Indeed, otherwise, it was shown that the integral in the estimate on p. 223 for 

the successive approximations {u'['\x,ty^ 



t 



\/W(t) <^{t), where (^(t) > A' J dt' + V{0) 

will never exhibit a necessary divergence (more precisely, here a Dini-Osgood-type integral 
condition is supposed to occur). In other words, Leray's estimate (3.9) above then suggests 
that a fast-type blow-up, which is not self-similar (i.e., of Type II), is the only possible. 

We devote a notable part of the present paper to discussing such type of fast blow-up 
that is unbounded in the rescaled variables introduced in fl2.7p and fl2.18p below. 

Thus, nonexistence (11. 4p of Leray's similarity solutions (12.71) and other related local 
types of self-similar blow-up is a definite step towards better understanding of the sin- 
gularity nature for the Navier-Stokes equations. This does not settle the problem of 
singularity formation, since there might be other ways for (II. ip to create singularities as 
t T~ rather than the purely self-similar scenario (12.70 : see the monograph |148j for 
details. Other concepts of such a multiplicity are discussed below. 

Concerning blow-up of infinite energy solutions, consider the strain field {q.v. |161j ): 

(2.12) u={-C{t)x + u{x,y,t),-C{t)y + v{x,y,t),2C{t)z), where ({t) = \\uj{;t)\U 
and the vorticity is structurally associated with famous Oseen's vortex (1910) |17Qj . 

(2.13) a;(r,t) = 2(5^e"T^ {r^ = x"^ + y"^) , so that curl u = (0, 0, w;(r, t)). 

The additional velocity field associated with an analogous vorticity in the cylindrical 
coordinates {r, ip, z} is pL59j (here F > is a constant) 

(2.14) u= (0,t;(r,t),0), where y(r,t) = 2^ (l-e-^) (u; = (0, 0, u;(r, t))) . 



"'^'^This Leray's integral approach can be characterized as a forerunner type for a number of later famous 
estimates via "control of vorticity" including that by Beale, Kato, and Majda [10] for Euler's equations 
(on recent extensions, see |223| ). 
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2.3. On blow-up in Euler's equations. A similar nonexistence of self-similai@ and 
some other types of blow-up has been obtained for the Euler eguations fll.311) in 
(derived by L. Euler in the middle of the eighteenth century [54j^i: see the most recent 
survey by Constantin ^42j for a modern mathematical activity exposition around); see 
[Mf H5l 1102] , and references therein. Local existence of smooth solutions for fll.Sip has 
been known from 1920s; see Lichtenstein |139j (and [8] for a full history); see also later 
results of the 1970s by Kato [113] . Existence of global classical solutions is still opeij^- 
Pelz |175j and Gibbon ^5] contain interesting surveys on various mathematical, symmetry, 
and more physical ideas concerning possible blow-up scenarios including numerical aspects 
(for the latter, see Kerr |116] for completeness concerning recent discussions on numerical 
blow-up issues). 

Infinite energy solutions of fll.311) do blow-up, for instance, according to the following 
separable solution in the cylindrical coordinates [57]: 

(2.15) u=(M(r,t),0,2;7(r,t)), where ^{r,t) = u{r,t) = . 

Blow-up of more general solutions of this form u{x, y, z) = {u{x, y, t),v{x, y, t), z'y{x, y, t)) 
was studied in detail in [H], where, in particular, non-power blow-up rate was observed, 

~ {T-t)\ ln{T-t)\ t ^ T . 

Notice also that, in a bounded (interior or exterior) domain in rescaled variables, Euler's 
equations fll.3ip admit non-trivial similarity solutions [97] (cf. (12.71) for a = /5 = |) 

(2.16) u(x, t) = jjTz:tj^ U(?/), y = jj^tji^ in the range /3 E [|, l] and a + (3 = 1. 

In the original spatial x-variables, this self-similar blow-up is supported by boundary 
conditions for a domain that shrinks into a point as t 1^. On the other hand, condi- 
tions (looking still rather non-constructive) of pointwise (in L°°) blow-up for fll.3ip were 
introduced in [35]. 

We believe that some ideas of our blow-up swirl analysis can be also applied to Euler's 
equations fll.3ip . and this is reflected in [TOJ. Moreover, in [70j, we perform the asymptotic 
construction of blow-up patterns of the NSEs fll.ip converging as t ^ to similarity 
solutions of the EEs (1 i.j'ip . However, in general, fll.3ip is a special subject that, in 
several circumstances, demands different approaches for constructing families of blow-up 
patterns. 



It seems that first strongly physically motivated arguments in favor of blow-up phenomena for Euler's 
equations were due to Onsager [169] (1949). 

^^In published form the incompressible equations appeared only in 1761, while a preliminary version 
was presented to Berlin Academy in 1752; see the full history in [40] , 

"The blow-up problem for the Euler equations is a major open problem of PDE theory, of far greater 
physical importance than the blow-up problem for the Navier-Stokes equation, which of course is known 
to nonspecialists because it is a Clay Millennium Problem," 42, p. 607]. 
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2.4. Singular blow-up set has zero measure. There exists another classic direction 
of the singularity theory for the Navier-Stokes equations that was originated by Leray 
himself |137j (see also details in 153]) and in Caffarelli, Kohn, and Nirenberg [26]. It 
was shown that the one- dimensional Hausdorff measure of the singular (blow-up) points 
in a time-space cylinder is equal to zero and these are contained in a space-time set of 
the Hausdorff dimension < i. We refer to [931 I164[ I192j for further development and 
references. In particular, at a given moment Iq, the admitted number of singular points 
can be finite; see [1631 1191] and presented references. 

Incidentally, among other results including Leray's one in [137j . a refined criterion is 
obtained in [192] . saying that, if T = 1 is the first singular (blow-up) moment for a 
solution u{x,t) of (11.11) . then 

1 

(2.17) lim_ / / \u{x,t)fdxdt = +oo. 

t R3 

The condition (12.171) is consistent with Leray-Prodi-Serrin-Ladyzhenskaya regularity L^''^ 
criteria and other more recent researches; see key references, history, details, and results 
concerning this huge existence-regularity-blow-up business around the NSEs in [531 l63l 
11341 11481 11821 11921 12021 1219J : [57J represents a modern panorama of such studies, which 
also commented that Ohyama's result in 1960 [166] was obtained before Serrin's one in 
1962 [194] . (I2.17P is also associated with Kato's class of unique mild solutions (in M^), 
|115] : see details and key references in [29l [63t |2T6] . 

2.5. Leray rescaled variables. We perform the nonstationary scaling as in (12. 7p . T = 1, 

(2.18) u(x,t) = ^7^u(t/,r), y = ^£=, ^ = _ in(l - t) ^ +oo as t ^ 1". 

This yields the rescaled equations for u = (u^ , u"^ , ii^)'^ and P, 

(2.19) = Au - |(y • V)u - I u - (u ■ V)u - VP, div u = in x M+. 
In particular, after scaling, (12.171) takes the form 

(2.20) ^lini^ e-r"e-^(/R3 \u{y, s)\' dy') ds = +oo, 

so that, if t = 1 is singular, then the solution of the rescaled equations (12.191) must diverge 
(blow-up) as r ^ +oo in L^(M^). 

As a standard next step, we exclude the pressure from the equations (12.191) . 

= H(u) = (B* - i/)u-P(u- V)u in x M+, 
^'^'^^^ where Pv = v - VA-i(V • v) 



is the Leray-Hopf projector of (L^(]R^))'^ onto the subspace {w G {L"^)^ '■ divw = 0} of 
solenoidal vector fieldo- Another representation is Pv = (t>i — Ria, f 2 — -R20", "^3 — -^30")"^, 



^■^This emphasizes the unpleasant fact that the NSEs are a nonlocal parabohc problem, so that a 
somehow full use of order-preserving properties of the semigroup is illusive; though some "remnants" of 
the Maximum Principle (cf. (|2.10p ) for such flows may remain and actually appear from time to time. 
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where Rj are the Riesz transforms, with symbols and a = RiVi + R2V2 + R^v^,- 

We then first apply P to the original velocity equation in (11.11) and next use the blow-up 
rescaling f l2.18p . Using the fundamental solution of A in M^, > 3 (ctat is the surface 
area of the unit ball Bi C M^) 

(2.22) = where = 

the operator in ( 12.2ip is written in the form of Leray's formulation [1481 p. 32] 
H(u) ^ (B* - i /)u - (u ■ V)u + / ^ tr(Vu(^, T)f d^, 



(2.23) 

where tr(Vu(z,r))2 = Y.^.^.^ ff^^i and = ^ > 0. 



2.6. Hermitian spectral theory of the blowing up rescaled operator: point 
spectrum and solenoidal polynomial eigenfunctions. In fact, the rescaled equations 
fl2.21l) (or the original ones (I2.19P ) are truly remarkable. Writing first linear operators on 
the right-hand side of (12.211) in a divergent form, 

(2.24) B*u = (B* - i /)u = Au - i (y ■ V)u - i u = ^ V • (p*Vu) - | u, 

where the weight is p*{y) = e ^ , we observe that the actual rescaled evolution is now 
restricted to the weighted L^-space L^*(M^), with the exponential weight p*{y)- Here, 
B* = B* — I / is a shifted adjoint Hermite operator with the point spectrum [I5i p. 48] 

(2.25) a(B*) = {Afc = -| - |, = |/3| = 0, 1, 2, ...} (/5 is a multiindex), 

where each has the multiplicity for = 3, or the binomial number C^^^_]^. 

The corresponding complete and closed set of eigenfunctions $* = {ip^iy)} is composed 
from separable Hermite polynomials. Similar spectral and eigenfunction properties can 
be also attributed to blow-up problems for 2mth-order parabolic PDEs such as (I1.56P : see 
e.g., [51], where a correspondence between the "blow-up" space L^, and for the global 
evolution, with p* = ^, is more clearly explained. The bi-ortho normality holds: 

(2.26) (^^,^'7) = 5/37 for any /5, 7. 
Note another important for us property of Hermite polynomials: 

(2.27) ^i'p, any derivative D^ip*^ is also an eigenfunction with k = \P\ — > 0. 
Recall that [I5] 

(2.28) polynomial set $* is complete and closed in L^, (M^). 

We need to consider eigenfunction expansions in the solenoidal restriction 

(2.29) ^p.(K^) = ^J*(K^)^ n {divv = 0}. 

Indeed, among the polynomials $* = {ip^^} there are many that well-suit the solenoidal 
fields. Namely, introducing the eigenspaces 

r, = Span{^;, \P\ = k}, k>l, 
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in view of fl2.27p div plays a role of a "shift operator" in the sense that 



(2.30) 



div : $f 



k-l- 



We next define the corresponding solenoidal eigenspaces as follows (see also Section [2^ : 

(2.31) SI = {v* = {vl,v^,v*^) : divv* = 0, v* e $^}, where dim = k{k + 2); 

see [m [82] and further references therein. Actually, [81] deals with global asymptotics 
as t ^ +00, where the adjoint operator B in fl2.49p occurs. Since B is self-adjoint in 
L^(M^), several results from [821 Append. A] are applied to B* (cf. Section [2^8!) . For a 
full collection, see [20] for further asymptotic expansions and extensions of these ideas. 

In particular, those solenoidal Hermite polynomial eigenf unctions of B* can be chosen 
as follows [82, p. 2166-69] (the choice is obviously not unique; normalization constants 
are omitted): 



Ai — — i 



'11 



" " 




2/3 




-2/2 


-2/3 


> V*2 = 







2/1 


. ^2 . 












(dim^i* = 3); 



(2.32) 





4 - 2/1 - yl 




ym 




2/12/3 


A2 = -1 : = 


ym 




4 - 2/? - 2/3 




2/22/3 




2/12/3 




2/22/3 




4 - 2/? - 2/2 
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-2/12/3 
2/12/2 



'25 



2/22/3 



-2/22/1 





-2/22/3 




2/12/2 




2/2 - 2/1 


V26 = 


2/22/3 




2/i - 2/? 




-2/12/2 




2/? -2/i 




-2/22/3 




2/12/3 



etc. 



(dimcSI = 8 

We need the following final conclusion. By (12.281) . the set of vectors is complete 



and closed iio L'^p* 



(2.33) 



Vv G Llj 



i^Y,so that 

I)3N3 



*3 L 



> 1. 



It then follows from (12:261) -( 1230D that 



(2.34) polynomial set $* = $ fl {divv = 0} is complete and closed in L^*(]R^). 

In what follows, we always assume that we deal with "solenoidal" asymptotics involving 
eigenf unctions as in (I2.3ip . 

For Burnett equations (11.61) . as we have promised to go with in parallel, the blow-up 
rescaling and elements of linear solenoidal spectral theory are found in Section 13. 2[ 



2^Note a standard result of functional analysis: polynomials are complete in any weighted L^-space 
with an exponentially decaying weight; see the analyticity argument in Kolmogorov-Fomin 1120} p. 431]. 
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2.7. On a countable set of quasi-periodic singularities: first formalities. Thus, 
according to the criterion fl2.20l) . the moment T = 1 is not a singular (and hence regular) 
point, if the corresponding locally smooth solution of f l2.2ip does not blow-up as r — +00 
in a suitable functional setting. Thus, the problem of global existence and uniqueness of 
a smooth solutions of the NSEs in reduces to nonexistence of blow-up in infinite time 
for the rescaled system fl2.19p or fl2.2ip . In such a framework, this problem falls into the 
scope of blow- up/non-blow- up theory for nonlinear evolution PDEs. 

Let us first discuss a simple corollary that follows from the above spectral properties of 
fl2.24p . Since by assumption T = 1 is the first blow-up point of u(x, t), we study solutions 
of that are globally defined in X i.e., do not blow-up in finite r. Moreover, 

the scaling (12.180 implies that we are looking for orbits with very sharp L^-divergence, 

(2.35) ||u(r)||^ = cie^^(l + 0(1)) as r ^ +00, where Ci = ||m(-, 1)||^ > 0. 
Consider the energy identity for smooth rescaled solutions 

(2.36) i^||u(r)||^dr=-Pu(r)||i + i||u(r)||^. 

Solving it together with (I2.35P yields the following control of the gradient Du[t): 

00 00 

(2.37) ||u(t)||2= (ci + 2/e"f||Du(s)||2ds)ei =^ / e"! pu(s) ||2 ds < 00. 

T 

In particular, for any e > 0, the following measure is always finite and satisfies: 

(2.38) meas{s>l: \\Du{s)\\l > eei] = o{l) as e^0+. 

Indeed, in comparison with (12.350 . this shows certain "degeneracy" for r ^ 1 of the L^- 
norm of the gradient Aii(r) relative to that of u(t). At least, this shows that the gradient 
cannot have a uniform exponential growth as in (I2.35P and should be much slower. In 
general, this does not imply essential consequences, since, according to (12.240 

the actual blow-up evolution of u(r) as r ^ +00 

(2 39) 

^' ' is restricted to the space L^, (M^), rather than L2(M^). 

Therefore, in particular, the exponential L^-divergence of the orbit (I2.35P on expanding 
spatial subsets as ?/ — 00 can be "invisible" in the metric L^* (M^) with the exponentially 



decaying weight p*{y) = e"'^'^/^. 

Thus, this is the class of global orbits of interest, and we are looking for the structure of 
these a;-limit sets, which, for such smooth orbits, are defined via local uniform convergence. 

As usual in dynamical system theory, one first discusses the case when the orbit {u(r)} 
approaches the simplest invariant manifold being a point. Thus, assume that 

(2.40) u(r) ^ u as r +00, 

where, for future use, we suppose convergence in L''(]R^), with g > 3 (note that a standard 
topology is expected to be that of L^*(M^) or the corresponding Sobolev one Hp*(R^)). 
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Then u is necessarily a self-similar profile, so that by global and local nonexistence results 

u = 0. 

Note that the spectrum fl2.25p of the linearized operator implies that, with a proper con- 
trol of the quadratically small convection term (this is easy in L^, with the exponentially 
decaying weight), the non- stationary small solutions u(r) must satisfy 

(2.41) |u(r)| ~ 0(e-^^) as r ^ +00, 

where Aq = — | is precisely the spectral gap in fl2.25p . More precisely, as a new application 
of spectral theory, we have the following: 

Proposition 2.1. Assume that, along a subsequence {rfc} —>■ +00, 

(2.42) u(?/,rfc)->0 uniformly in L°°nL^p,. 

Then t = T is not a blow-up time for u(x, t) (in other words, the singularity is removable) . 

Proof. Consider the sequence of solutions {uk{y, s) = u(y, Tk + s)} with vanishing initial 
data in L°° according to (12.421) . Using and well-developed spectral properties of the 
linearized operator B* in (12.241) defined in L^*(R^), with generalized Hermite polynomials 
as a complete and closed set of eigenfunctions. Therefore, according to classic asymptotic 
parabolic theory (see e.g., [145] ). we conclude that for any sufficiently large k, 

(2.43) Ufc(?/,s) ~ O(e-i), s>l =^ u(y, r) ~ O(e-i), r > 1. 
Overall, taking into account Leray's scaling (12.181) . this yields (T = 1): 

(2.44) u(a;,t) ~ (T-t)-50(e-5) = 0(1) as t^T', 

so that u(x, t) is uniformly bounded at t = T. □ 

Thus, equation (12.211) does not allow stabilization to an equilibrium, since this corre- 
sponds to the no-blow-up case. Further, as usual in textbooks on dynamical systems, the 
next candidate for being the corresponding invariant manifold is a periodic orbit of finite 
period > (S- In |178l p. 1218], this conjecture was connected with the study of the 
complex Ginzburg-Landau equation 

(2.45) iMt+ (1 -i£)Au + (1 + i5)|np'^n = / in M^x(0,T), 

where e > 0, 6 > and a > j^. It was pointed out in [178] that the model (12.451) exhibits 
the same scaling, similar energy control, and local semigroup theorj@; see also [149[ § 4] 
for some related estimates of such a periodic behaviour of an unknown structure. It was 
pointed out in [220] (see also [141] ) that the CGLE such as (12.451) can be a good PDE 

^^The role of r-periodic motion for rescaled blow-up solutions of the Euler equations (|1.3ip was pointed 
out in j l80] (however, the nonexistence conclusion on p. 218 therein looks rather suspicious). 

^''However, unlike (jl.ip . it was shown in [178] and in other papers cited therein that (|2.45p admits 
a lot of (countable families of?) blow-up similarity profiles, and this makes it more analogous to the 
Cahn-Hilliard model (I1.57P : see [56], where a countable blow-up family was detected even for = 1, 
and [70j . where Leray's scenario was shown to apply. 
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system modelling regularity and other questions regarding the NSEs (11 .ip . In particular, 
global existence of weak solutions can be obtained analogously to Leray's proof. Inciden- 
tally, concerning the singular set of blow-up points for (12 .45^ of zero measure [220j, it is 
known that this set is restricted to a bounded domain, |14I] . 

Actually, this would mean proving for (I2.2ip the Poincare-Bendixson theorem saying, 
essentially, that if a rescaled orbit of (I2.2ip satisfying (I2.36p . (I2.37P does not stabilize 
to an equilibrium (and actually (I2.44p makes this impossible), then it convergence to a 
simple closed curvcJ, i.e.. 



(2.46) t^(uo) consists of a T,,-periodic orbit Fi. 



According to (I2.39p . the natural metric of convergence as r ^ +oo in (I2.46p is then 
assumed to be that of L^, , while the L^-divergence (I2.35P occurs on subsets with \y\ ^ 1, 
which does not affect the convergence. A more clear discussion of such blow-up periodic 
orbits will be postponed until Section [731 where a proper spectral theory is under scrutiny. 

Of course, (12.461) would be the best and very pleasant case. Indeed, for general DSs of 
such complexity, (I2.46P is a very difficult open problem. Hence, the first main principal 
difficulty is how to predict a possible structure of such blow-up "periodic orbits" . Actually, 
this is one of our main goals, and the blow-up angular swirl mechanism (11.251) in Section 
is a natural argument in support of the periodic motion (I2.46P . 

Continuing using the logic of standard dynamical system theory, we next would have to 
assume that a periodic ti;-limit set would have been ruled out (i.e., being nonexistent) for 
a given rescaled orbit {u(r)}. Then we should, e.g., conjecture a countable set of other 
possibilities of evolution geometrically related to tori in M""*"^: for n = 2, 3, 4, ... , 



(2.47) "-^(uq) consists of a quasi-periodic orbit r„ driven by n fund, frequencies. 



where n = 1 leads to (12.460 . To get a quasi-periodic motion for n = 2 on the invariant 
a;-limit set, we introduce in Section [8] the idea of precessions of the swirl axis, which leads 
to extremely difficult and open mathematics. For n > 3, such a clear visual geometric 
interpretation of the scenarios (I2.47P is not that easy or straightforward. A spectral 
background for bifurcation of such patterns is difficult and obscure (cf. Section [7^41) . 

Eventually, under the assumption of nonexistence of all of those types of cj-limits in 
(12.471) . passing to the limit n — oo would then have led to a kind of a strange attractor 
for the dynamical system ( 2.2ip . which can be an extremely complicated invariant man- 
ifold to be proved to exiscj. Recall that, for proving nonexistence of L°°-singularities 



^"^Is there any hope that a nonstationary version of the eUiptic differential inequahty (|2.f Op apphed to 
the nonlocal parabolic PDE (|2.2ip can rule out at least some of special quasi-periodic oscillations about 
the unique equilibrium 0? (seems, not sufficient, and no hope). 



^^The proof of existence of a robust strange attractor for the E. Lorenz dynamical system in R'^ [144] 
("3 X 1" , i.e., uncomparably easier), describing thermal fluid convection with some relations to the NSEs, 
proposed in 1963 was declared by Smale as one of the several challenging problems for the twenty-first 
century (1998), which eventually took nearly four decades to complete; see [2091 1210] and also [153] . 
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for (11.11) . all these infinite number of possibilities have to be ruled out (via a new en- 
ergy/monotonicity/spectral, etc. control). 



2.8. Global similarity solutions defined for all t > do exist. Solenoidal 
linearized patterns. In contrast to blow-up, as usual for typical parabolic reaction- 
diffusion-absorption problems, global similarity solutions of (11.11) without blow-up, i.e., 
(12. 7p with T—t^t>0 occur more frequently and correspond to the following scaling and 
rescaled equations (the invariant of the scaling group involved, y = x/\/t, is sometimes 
called the Bolzman suhstitutioi)Q\: 



-oo, 



u(x,t) = ^u(?/,r), y=fi, r = Int +00 as t - 
where = Au + i (y • V)u + ^ u - P (u ■ V)u in R-'^ x IR+. 
Here we obtain the linear operator 

(2.49) B = B - /, 

where B is adjoint to B* in the metric of the dual space L^(M^). The spectrum of B, 
which has a self-adjoint (Friedrichs) extension in L^, where p = [ISj, is 

(2.50) a(B) = {Afc = -| - 1, k= |/?| = 0, 1, 2, ...}. 

Non-trivial self-similar profiles, i.e., stationary solutions of (I2.48p . do exist and describe 
asymptotics as t — *• +oo of various sufficiently small solutions; see [301 1291 HH mi [82], |89l 
[92[ I177j . In [177] . a simple criterion of asymptotic similarity form for solutions was 
obtained. Note that Slezkin-Landau singular stationary solutions (17. ip are self-similar. 

In addition, by (12.501) . is asymptotically stable, and this makes it possible to construct 
fast decaying solutions on each ID stable manifolds with the asymptotic behaviouS 

(2.51) U/3(x,t) ~ t^'=~W/3(^) + ... as t oo, where V/3 = v|,F G 5^ 
are solenoidal eigenf unctions of B defined in Section 12. 6[ Namely, taking 

(2.52) V = {v,,V2.v^Y e Sk, v,e<^k = {i^p = ^D^F{y), \P\ = k}, 
where F stands for the rescaled Gaussian (see (15.601) with = 3), we have that 

(2.53) div V = {v,)y.^ + {V2)y, + {v;)y^ = div (v*F) = (div v*)F -ly.^r*F. 

This establishes a one-to-one correspondence between solenoidal eigenfunction classes SI 
in (I2.3ip for B* and Sk in (12.511) for B; see (I2.32p for the first eigenfunctions V/3 = v^F. 
Therefore, dim Sk = k{k+2), etc.; see details and rather involved proofs of the asymptotics 
(12.511) for /c = 1 and 2 in |8lj- We will deal with patterns such as (12.511) later on. 



^'^Similarity solutions were used by Weierstrass around 1870, and by Bolzman around 1890; this rescaled 
variable y in parabolic PDEs was widely used by Sturm in 1836 ^203j (and possibly even before?) 

""^We present here only the first term of expansion; as usual in dynamical system theory, other terms in 
the case of "resonance" can contain In i- factors {q.v. 4 for a typical PDE application); this phenomenon 
was shown to exist for the NSEs in ^82, p. 236]. 
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It follows from (12.321) that there exists the corresponding eigenf unctions of B with 
Ai = — I given by (I2.50p . Hence, by scaling (12.481) . for instance, there exists the asymptotic 

pattern for t ^ 1 (the rate 0(t~^) in is thus sharp for uq G I/^(M5j_)) 

(2.54) ui(x,t) ~ ^vi(^), where vi(y) = ^ (2/3, 2/3, 0)^ e"^. 

In the half space x G M^n{a;3 > 0} with no slip boundary condition u|j:3=o = 0, the decay 
rate of solutions is also of interest (clearly, some of the polynomials in (12.321) are good for 
that); see [39] for recent developments. 

Note again that the calculus of solenoidal eigenfunction classes look like being specially 
designed to suit the divergence-free flows not only for both types of scalings, blow-up 
(12.181) and the global one (I2.48p . but also for the Burnett equations (II. 6p . where classes 
Sk can be defined for any necessary 2mth-order linear operators [51]. The adjoint ones 
then are also composed from solenoidal generalized Hermite polynomials only. 

3. First application of Hermitian spectral theory: Sturmian local 
structure of zero sets of bounded solutions and unique continuation 

3.1. Nodal sets for the Stokes problem and NSEs. Here we perform a first step to- 
wards the classification problem ( ]1.20p . Namely, we assume that at the point (x, t) = (0, 1) 
the solution u(y, t) is uniformly bounded and is such that the eigenfunction expansion of 
the corresponding rescaled function satisfying (I2.40p . 

(3.1) u(y,r) = E(/3)C/3(r)v;(?/), where c^v^ = {civ;„ C2V*p^, CsV*^^^ G ^^.(RS)^ 

converges in L^, (R'^), and moreover, uniformly on compact subsets. These convergence 
questions of polynomial series are standard; see [511 [SZ], where further references and 
details are given. Then the expansion coefficients satisfy the following dynamical system: 

2) = (^/^ ~ ^)^^ + ^(°.7) ^"7/3CaCT, for any |/?| > 0, 

\ where = -(P (v^ ■ V)v;, v^) for all a, 7. 

It is natural to assume that the quadratic sum on the right-hand side converges for the 
given smooth rescaled solution u(?/,r). Recall that, according to scaling (12.180 . we deal 
with bounded and uniformly exponentially small functions satisfying 

(3.3) \u{y,T)\ < Ce~^ in R^ x R+. 

The system (13. 2p is difficult for a general study, and, of course, it contains the answer to 
the existence/nonexistence problem, provided that (0, 1) is a singular point of the solution 
u{x,t). For regular points, it can provide us with a typical classification of nodal sets of 
solutions. This kind of study was first performed by Sturm in 1836 for linear ID parabolic 
equations |203] : see historical and other details in [67] Ch. 1]. 

Thus, following these lines, we clarify local zero sets of solutions of the NSEs at regular 
points. Assume that 

(3.4) u(0,l) = O, 
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which can be always achieved by constant shifting u{x,t) u{x,t) — u(0, 1). In this 
connection, recall that the first eigenf unctions with X/s = 

(3.5) v*(|/)~(l,l,l)^ (1,1, 0)^ (1,0,1)^,..., 

are the only ones that have empty nodal sets of some of its components. Then, bearing 
in mind the blow-up scaling term {1 — t)~^ = in (12.181) . we have to assume that 

(3.6) Co(r) = or Co(r) ^0 as r — > +oo exponentially faster than e~5. 

Polynomial nodal sets for the Stokes problem. A first clue to a correct understanding of 
the DS ( 13. 2p is given by the Stokes problem, i.e., without the nonlinear convection term, 

(3.7) Ut = - Vp + Au, div u = 0. 
Then (13.21) becomes linear diagonal and is easily solved: 

(3.8) c^=(A^-i)c/3 =^ cpir) = cp{0)e-^^^ for any |/?| > 0. 

Therefore, according to (13. ip . all possible multiple zero asymptotics for the Stokes problem 
(its local "micro-scale turbulence" ) is described by finite solenoidal Hermite polynomials, 
and the zero sets of rescaled velocity components also asymptotically, as r — > +oo (i.e., 
t 1~) obey the nodal Hermite structures. 

NSEs . Consider the full nonlinear dynamical system (13. 2p . which on integration is 

(3.9) Cp{T) = Cf3{0)e 2 -e 2 J 2^(^ ^) rfa7/3(CaC^)(s)e 2 ds. 



It follows that the nonlinear quadratic terms in (13. 9p . under certain assumptions, can 
affect the rate of decay of solutions near the multiple zero. As usual in calculus, this 
indeterminacy can be tackled by L 'Hospital rule. 

Since we are mainly interested in the study of nodal structures of solutions by using the 
eigenfunction expansion (13. ip . we naturally need to assume that it is possible to choose 
the leading decaying term (or a linear combination of terms) in this sum as r — > +00. 
Then obviously these leading terms will asymptotically describe the Hermitian polynomial 
structure of nodal sets as t ^ 1~. For PDEs with local nonlinearities, this is done in 
a standard manner as in [671 § 4]; in the nonlocal case, this seems can cause technical 
difficulties. However, the DS (13.21) looks (but illusionary) as being obtained from a problem 
with local nonlinearities. In other words, the nonlocal nature of the NSEs is hidden in 
(13. 2p in the structure of the quadratic sum coefficients {(ia,7,/3}, and this do not affect the 
nodal set behaviour for some class of multiple zeros. We will check this as follows: 

We consider a "resonance class" of multiple zeros. Namely, let us assume there exist a 
multiindex subset B and a function h(r) such that 

Cfl(T) ~ h(r) as r ^ +00 for any 3 E B, 

(3.10) I^K J K J 

|c/3(r)| ^ |h(r)| as r ^ +00 for any (3 ^ B. 

In other words, only the coefficients {cpi^r), /? G i3} are assumed to define the nodal set 
via (13. ip . and other terms are negligible as r — +00. Under the natural assumption of 
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a strong enough convergence of the quadratic sums in fl3.2p (this can be expected not to 
be the case for singular blow-up points only), taking the ODEs from fl3.2l) for each (3 E B 
yields, for r ^ 1, 

(3.11) = (A/3 - |)c^ + o(c^), where C/3(r) ~ h(r). 

Hence, the asymptotic balancing of these equations must assume that as r ^ +oo 

(3.12) h~(A/3-i)h =^ c/3(r) ~ h(r) ~ e(-|-5)- and \P\ = k, 

where we may omit lower-order multipliers. Thus, there exists a k > 1 such that \f3\ = k 
for any (3 E B. One can see that for such "resonance" multiple zeros, the nonlocal 
quadratic term in (13. 2p is not important. Thus, in the resonance zero class prescribed by 
(13.101) . as r ^ +oo, on compact subsets in similar to Stokes' problem. 



(3.13) the nodal set of u(y, r) is governed by some solenoidal Hermite polynomials. 



Note that the conclusion that, locally, for any zero of finite order at (0, 1), 

(3.14) nodal sets of u(a;,t) are governed by finite-degree polynomials 

is trivially true for any sufficiently smooth solution. Indeed, this follows from the Taylor 
expansion of such solutions 

(3.15) u(x, t) = ^ (0, 1) (1 - tr + R(x, t), 

where R stands for a higher-order remainder. Translating (13.151) via (I2.18P into the 
expansion for u{y, r) yields some polynomial structure, so (13.140 is obviously true. Thus, 
the principal feature of ( 13.130 is that the Hermite polynomials count only therein. 

Obviously, for the nonlocal problem (I2.2ip . there exist other non-resonance zeros. In- 
deed, let (0, 1) be a zero of u(x, t) of a finite order m > 1, i.e., as x — 0, 

(3.16) u(x, 1) ~ x"^, with \a\ = m. 
We now use the following expansion: 

(3.17) u(x, t) = u(x, 1) - uj(x, 1)(1 - t) + |r u«(x, 1)(1 - tf + ... , 
where, by (12.210 . all the time-derivatives Dfu(a;, 0) can be calculated: 

(3.18) Ui(x, 1) = Au(a;, 1) + (P(u ■ V)u)(a;, 1) ~ s'"-^ + (P(u ■ V)u)(a;, 1), 

with a natural meaning of Ax" ~ x'^^^. If the nonlocal term is negligible here and for 



other time-derivatives, i.e.. 



Ui(x, 1) ~ x" ^ Uit(x, I) x" ^ 

then according to (13.171) this leads to a Hermitian structure of nodal sets. In fact, this cor- 
responds to the pioneering zero-set calculus performed by Sturm in 1836; see his original 
computations in [66^, p. 3]. 
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In general, the nonlocal term in (13.181) is not specified by a local structure of the zero 
under consideration, so, obviously, can essentially affect the zero evolution. For instance, 
as a hint, we can have the following zero: 

(3.19) ut(0, 1) = C ^ =^ u{x, t) ~ X- + (1 - t) ~ e--(l +z''), z = 

so this nodal set is governed by the rescaled variable which us different from the 
standard similarity one y in (I2.18p . Of course, due to the nonlocality of the equation, 
many other types of zeros can be described. Actually, such non-resonance zeros can be 
governed by sufficiently arbitrary polynomials as the general expansion (I3.15P suggests. 

Finally, the proof that zeros of infinite order are not possible (and, as usual in such 
Carleman and Agmon-type uniqueness results, this occurs for u = only) is a difficult 
technical problem; see an example in [671 § 6.2]. For analytic in y solutions of the NSEs 
(see references and results in ^9} I224j ). this problem is nonexistent, and then in (13.131) 
the degree of the solenoidal vector Hermite polynomials is always finite, though can be 
arbitrarily large. 

Note another straightforward consequence of this analysis that this gives the following 
conventional unique continuations result: let (13. 4p hold, (0, 1) he a resonance ^erto, and 
at least one component of the nodal set ofu.{y,T) does not obey (13.131) . Then 

(3.20) u = everywhere. 

Of course, this is not that surprising since the result is just included in the existing and 
properly converging eigenfunction expansion (13. ip under the assumption (I3.10p . 

For elliptic equations P{x, D)u = 0, this has the natural counterpart on strong unique 
continuation property saying that nontrivial solutions cannot have zeros of infinite order; 
a result first proved by Carleman in 1939 for P = —A + V, V E L'^^, in [32]; see 
\5U\ 1205] for further references and modern extensions. 

Thus, this is the first application of solenoidal Hermitian polynomial vector fields for 
regular solutions of the NSEs. We expect that, due to the DS (13.20 . some "traces" of 
such an analysis and Hermite polynomials should be seen in the fully nonlinear study of 
u{y,T) at the singular blow-up point (0,1), where, instead of (13.40 . we have to assume 
that, in the sense of limsup^^^, 

(3.21) |u(0,l)| = +CX). 

3.2. Burnett equations. For (11.60 . the blow-up scaling (12.180 is replaced by 

(3.22) u(x,t) = (l-t)-tu(y,r), y = ^^^^ 
so that the rescaled system (12.210 takes a similar form 

(3.23) u^ = H(u) = (B*-|/)u-P(u- V)u in x M+. 



Indeed, this is hard to check; for PDEs with local nonlinearities, this assumption is not needed, so 
that such a unique continuation theorem makes full sense, [67] . 
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The spectral theory of the given here adjoint operator 

(3.24) B* = -A2-iy.V, where a(B*) = {A^ = -f , |/?| = 0, 1, 2, ...} 

with eigenf unctions being generahzed Hermite polynomials is available in [51]; a solenoidal 
extension in the same lines is needed. Therefore, under the same assumptions, the poly- 
nomial structure of nodal sets is guaranteed for the corresponding Stokes-like and Burnett 
equations (and for an arbitrary 2mth-order of the viscosity — (— A)™u therein). 

4. Second application of blow-up scaling: on convergence to the EEs 

We now present the results of a general application of another related blow-up scaling 
for establishing a connection between blow-up in the NSEs and singularities in the EEs. 
We use the blow-up scaling in the form of [80, § 2], which received other applications 
and extensions in [681 ESl [H] , etc. As usual, such a rescaling near blow-up time leads to 
ancient solutions in Hamilton's notation [91] , which has been a typical technique of R-D 
theory; see various form of its application in [1871 [79] and others. In [71], this scaling 
technique applies to the NSEs and Burnett equations in to present a simple treatment 
of the corresponding Leray-Prodi-Serrin-Ladyzhenskaya regularity L^''' criteria and other 
estimates. Ancient solutions of the 3D NSEs allowed recently to get new non-blow-up 
results for axi-symmetric flows; see |118j and also Section I5.12[ 

We begin with a definition, which settles an "evolution" concept of entropy solutions for 
the EEs. In blow-up R-D theory, such concepts of "entropy- viscosity" come from extended 
semigroup theory, where proper (blow-up or singular) solutions are only those, which 
can be obtained by regular approximations. For problems with the MP, such extended 
semigroup theory leads to the unique continuation (partially or completely unbounded) 
of any blow-up solutions beyond blow-up time t > T, [66l § 6.2]. 

The questions of the vanishing viscosity limits in the NSEs to get the EEs are classical in 
fluid mechanics and lead in general to a number of fundamental open problems; see a clear 
statement of such questions in [SI p. 422], where necessary topologies of convergence are 
prescribed. For general (not necessarily bounded) solutions, the estimates (12. 2p suggest 
the weak-* topology of L°°(M_, L^(M^)), while for bounded solutions this can be improved 
(the reason for the uniform boundedness is associated with the blow-up scaling to be 
applied). However, we are not obliged to deal with a specific convergence, especially since 
the necessary minimal topology for in what follows is still unknown. 

Definition 4.1. A function ("distribution") u(-,t) is said to be a bounded NS-entropy 
(i.e., Navier-Stokes- entropy) solution of the EEs i \1.31\i if it is obtained as the limit 

(4.1) Ufc ^ u as k ^ oo weak-* in L°°(M^ x M_) 
of a sequence of bounded classical solutions {uk{x,t)}, 

(4.2) |ufc(2;,t)|<l and ||ufc(-, t) Ha < C for k > 0, 
of the NSEs with vanishing viscosity coefficients, 

(4.3) Ufc : Ut + P(u ■ V)u = 5fcAu, where 6k —* 0^ as k ^ oo. 
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The given in (14. ip weak-* topology of convergence follows from the fact that, by fl4.2p . 
the sequence {u^} is bounded in L°°(M^ x ]R_). Note that we do not specify in which 
sense then u satisfies the EEs, 

(4.4) ut + P(u ■ V)u = 0, 

and we even do not demand u to be any kind of weak solution of fl4.4p . Note that this is 
not always the case even for nonlinear dispersion equations (NDEs) such as [69l [75] 

(4.5) Ut = {uUx)xx, Ut = {uUxx)xi Ut = uUxxx, ^ud Ut = —{uUx)xxxx, etc. 

On the other hand, the regularity concepts are not always an option in singularity theory: 
it is known that sometimes even analytic extensions of solutions are not a proper one, 
i.e., proper (extremal) solutions have finite regularity; see an example in [66| P- 140]. 

Thus, according to (14. ip . the concept of the NS-entropy for the EEs includes the (nat- 
ural, indeed) way of regular approximations of its solutions. In this sense, this is quite 
similar to conservation laws theory, e.g., for the ID Euler equation 

(4.6) ut + uux = in M x M+, 
with L^-data uq. The parabolic approximation is as follows: 

(4.7) u = lime^o+ Ue, where Us : Ut + uUx = euxx, u{x, 0) = uo{x), 

with the convergence in L^. Of course, classic entropy concepts apply directly to (14. 6p (say, 
in the sense of Oleinik and Kruzhkov), but the entropy regularization description (14.70 is 
completely self-consistent; see |198j for details. Note that conservation laws such as (14.60 
are natural zero-level models for the 3D EEs, where both features of the divergence and 
the L^-control are available. Then Tartar-Murat's compensation compactness approach 
could be also natural, but, unfortunately, seem not applicable to singularly perturbed 
EEs dO]) . 

Further "distributional" properties of NS-entropy solutions of the EEs are unknown, 
and these compose the core of the problem; see below. 

4.1. The NSEs. Thus, we assume that there exist sequences {tk} —>■ < oo, {xk} C 
M.^ , and {Ck} — > +oo such that the solution u{x,t) of (II. ip becomes unbounded: 

(4.8) sup |u(a;, t)| = |u(a;fc, tfc)| = Cfc — +00 as k ^ oo. 

R^x[0,ifc] 

As in [801 § 2] , we then perform the change 

(4.9) Uk{x,t) = u{xk + x,tk + t) = CkWk{y,s), where x = akV, t = als, 
where the sequence {ak} is such that the L^-norm is preserved after rescaling, i.e., 

(4.10) \\Mt)\\2 = ||wfc(s)||2 =^ ak = C'J' ^ 0. 

Taking the NSEs in the nonlocal from as in (I2.2ip . we then obtain the following rescaled 
equations for w = Wfc(|/, s)\ 

(4.11) w, + 5k P(w ■ V)w = Aw, where 4 = c|. 
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Next, after time shifting, s t— > s — Sq, with a fixed arbitrarily large Sq > 0, the solutions 
and data satisfy the uniform bounds: for all A; ^ 1 

(4.12) |wfe(s)|<l and ||wfc(s)||2 < C for all sG[-So,0]. 

The principle (and obvious, otherwise global existence would trivially follow) fact is that 

(4.13) 6k = +00 as k oc (actually meaning that 7^ L^), 

so that (14. lip is a singularly perturbed problem to be analyzed as follows: 

We divide the equation (14.111) by Sk and introduce the new time s = 6kS to get for 
Wfc = Wfc(s) the equation 

(4.14) W5 + P(w • V)w = ^ Aw. 

In view of (14.121) . the sequence is converging weak-* in L°°(]R'^ x R_): 

(4.15) Wj!c(s) ^ w(s) as k 00 

(the convergence also takes place in better topologies). We then need to pass to the limit 
in the rescaled NSEs (14.141) . where the right-hand side vanishes in the weak sense in view 
of the sufficient regularity (14.121) . Concerning the quadratic convection term, in view of 
its divergence (12.41) . one needs extra assumptions to get in the limit a weak formulation 
of the resulting EEs (locally, for sufficiently smooth solutions, this is well-known due to 
pioneering results by Kato (1983), Temam and Wang, etc., see survey 0, § 4], [217J, and 
references therein). Using Definition 14. 1^ as a standard conclusion from (I4.14p . we obtain: 

Proposition 4.1. Under the above hypotheses, the following holds: 

(i) w in l \4-i5^ is a bounded ancient NS- entropy solution of the EEs 

(4.16) w, + P(w-V)w = for s < 0, sup^e^a |w(i/, 0)| = 1; 

(ii) w is not a classic regular solution of the EEs in x M_ . 

Concerning (ii), let us mention that the smooth EEs (I4.16P are a gradient (and sym- 
plectic) system with the positive definite Lyapunov function 

(4.17) iilM' = (<0), 

which is monotone on evolution orbits. Therefore, the omega-limit set in the topology of 
Cioc(IR^) of smooth orbits consists of equilibria, on which V{w) = J |wp vanishes locally 
on compact subsets, so w = 0. In other words, we use the fact that smooth L^-solutions 
of EEs decays to zero in time uniformly on compact subsets. Therefore, such solutions 
cannot satisfy the last normalization condition in (I4.16p . Actually, we do not need a 
deeper discussion here about these rather obscure aspects of EEs theory, since, moreover, 
writing (I4.14p without an extra time-rescaling, i.e., keeping the variable s, 

(4.18) iw, + P(w- V)w = iAw. 

yields as — 00 that w is a weak stationary bounded NS-entropy solution of the EEs. 
Finally, one can replace (11.20 by the following formal: 
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Corollary 4.1. // the only bounded stationary NS-entropy solution of the EEs with the 
additional to ( 14. ip uniform convergence on compact subsets is trivial, then ( \1.2\\ holds. 

The convergence fl4.15p embraces larger compact subsets in x than the self-similar one 
according to i \2.7^ : 

(4.19) \x-Xk\ = 0{C,^^), while \x - x^l^ = 0{C,'). 

Admissible solutions of the problem (14.161) are generally unknown. Thus, (I4.15P shows 
a general relation between blow-up in the NSEs and NS-entropy solutions of the EEs: 
special bounded singular solutions of the EEs must create and support blow-up for the 
NSEs, which eventually must evolve as t ^ T~ on compact subsets that are smaller than 
those in ( 14.19^ for the EEs. It general, the "almost trivial" case (cf. Corollary 14. ip 

(4.20) w(?/) =0 in \ {0} 

cannot be excluded. Obviously, for (I4.20p to be treated seriously, one needs an extra 
stronger 'W-topology" of convergence in (14. ip . which should be adequate to the both 
NSEs and EEs and is unknown. On the other hand, in |7^, we present almost "explicit" 
blow-up infinite energy patterns, where nontrivial limits like (I4.15P take place. Some 
principles of formation of Type II blow-up patterns are discussed in Section 17. 6[ 

4.2. Burnett equations. This is similar, consult [71] for final estimates. 

5. Construction of blow-up twistors on a 2D linear subspace in 

In this section, we demonstrate the fact coming from reaction-diffusion theory that the 
Navier-Stokes equations can admit blow-up behaviour generating a blow-up swirl (ll.24p 
at a stagnation point with accelerating and eventually infinite angular speed as t — ^ T~ . 
Such a blow-up behaviour is associated with the so-called logarithmic travelling waves in 
the tangential angular direction, which are group-invariant solutions that occur for some 
nonlinear diffusion-combustion PDEs. We explain such a behaviour in the cylindrical 
geometry, where the axis of rotation is fixed to be extended to more realistic spherical 
geometry in Section [HI To this end, we will perform a partially invariant construction 
(understood not in a standard way of invariance under a group of transformations). 

Even on the 2D subspace W2 given in (ll.23p . the analysis of existence of such blow- 
up twistors is very difficult and leads to involved nonlinear systems that we are not 
able to study rigorousljjf^. In view of this, we will show another hypothetical way of 
constructing other related non-similarity patterns of unknown rotational-like singularity 
nature in conjunction with the multiplicity claim (ll.29p . 

We begin by noting that the Navier-Stokes equations (ll.ip inheriting for given data 
a strong spatial symmetry (for instance, cylindrically axisymmetric irrotational flows), 
actually reduce their "effective" dimension by one (a very formal and rough issue) and 
then get into the global existence and uniqueness case of classical solutions for N = 2. 



'This blow-up swirl behaviour is naturally attached to the Poincare-Bendixson-type conclusion (|2.46p . 
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This idea goes back to Ukhovskii-Yudovich |211j and Ladyzhenskaya |130j concerning 
axisymmetric geometry; see [193^ I222j for proper detailed definitions of such global smooth 
flows and recent results and [147j for helical symmetries. Formally speaking, the present 
patterns with some features of axisymmetry could be not complicated enough to generate 
singularities with finite kinetic energy, so we will need further mental efforts to remove 
such a restriction to blow-up; see Section [HI On the other hand, global regularity of 
general axisymmetric flows is also not knowio, so that blow-up patterns may be possible 
even in this simplified geometry of the linear dependence on 2;, but, probably, after extra 
precession-type or other necessary manipulations with this axis of rotation; see below. 

5.1. Navier— Stokes equations in cylindrical coordinates. We introduce the stan- 
dard cylindrical coordinates (r, (^9, z) in M^, where 

(5.1) r"^ = x'^ + y'^ and e^ = (sin</), cosy^, 0), = {— cos ip, simp, 0), e2 = (0,0,1). 
Denote the corresponding velocity field and the pressure as follows: 

(5.2) u = {ur, v^, Wz) = (f/, V,W) = Uer + Ve^ + We^ and P, 

where V stands for the swirl component of velocity. Then equations (11. ip take the form 



(5.3) 



( iU-lV' = -Pr + AU-^V^-^U, 
iiV+'-UV = -lP^ + AV + ^U^-^V, 
i-^W = -P, + AW, 



^Ur + lU + }:V^ + W, = 0. 
Here, the full time-derivative ^ and the Laplacian A are given by 

^ = Dt + UDr + i VD^ + WD,, 

/r /|'\ At ' r ^ ^' 

being partial derivatives, and A2 the Laplacian in the polar variables (r, yj). 

5.2. On axi-symmetric flows: towards global existence. A flow (or a vector field) 
is called axi-symmetric if it is invariant under rotations in ip. The principal advantage of 
such flows is that the \^-equation in (15.31) takes the form of a standard parabolic one, 

(5.5) Vt + UVr + WV, = A,y - (1 f/ + ^)V, 

so that some of the MP, Harnack's inequalities, Nash-Moser iterations, De Giorgi-Nash 
estimates, and other tools of classic parabolic theory can be expected to take a part. For 
instance, the standard change for the 3D Laplacian gives the divergent equation: 

(5.6) rV = V =^ Vt + UVr + WV, + f K = A,1/, 
and hence by comparison, via the control on the parabolic boundary, 

(5.7) \V{r,z,t)\<C =^ \V{r,z,t)\<^. 



'Though, it is indeed an advantage to know that, by ^55^, singularities may occur at the z-axis only. 
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A more involved use of the MP yields some important results; e.g., to prohibit Type I 
blow-up mSlEB]: 

|u(x, t)| < J^_^ =^ no blow-up at t = T~. 

Note again that axi-symmetric flows with no swirl, i.e., with V = 0, are regular; see 
Ukhovskii-Yudovich |211] . Ladyzhenskaya [130] . and more recent new results in |193[ 
\222\ \118\ [38] (the regularity remains under the presence of helical symmetries |147 ] ) . 

Thus, for axi-symmetric settings, the V^-equation is simple, ( 15. Sp . while the third In- 
equation in (15. 3p gives the pressure, 

(5.8) P = J-{iW-AW) (W^L_=0). 
Finally, we calculate W from the last div-equation, 

(5.9) = r ^ irU)r. 
Plugging all this into the [/-equation yields 

Ut + UUr + WU, = AU+lV'-^U 

+77 m i^U)r)^ - (i (rf/,).)^ - [U{1 (r[/).) - [i {rU)4 (rf/)..] J. 

2 2 

As we have seen in (12.211) . the first two integral linear terms on the right-hand side in 
(15.101) can be incorporated into the main derivatives AU and Ut, since these have the 
necessary good signs, so that the equation reads 

(J + A)f/t + uUr + wu, = (A + C2)u + }y^-^u 

^'•''^ -77{[U{HrU)r)X+[HrU)rHrUUl], 

2 2 

where £i > and £2 < (in the metric of L^) are pseudo-differential operators with 
easily computed symbols (these are directly related to the projector P) 

(5.12) A = -Dr{\Dr{rI)) = -Ar + ^ /, A = ^r(Ai D,(r/)). 

Let us present some typical RD-like speculations concerning blow-up. Assume that, in 
the {U, l^)-system (15. lip . (15. 5p . a blow-up occurs at the origin r = 0,2 = 0ast^ T^. 
Then, in view of the L^-boundedness, this must be single-point blow-up. Then, the 
linear terms in the system are not essential as t — >■ T~ in comparison with a number of 
quadratic ones. This also related to the linear nonlocal terms, where integration over the 
necessarily rescaled spatial variables diminishes their rates of divergence. Therefore, the 
quadratic pointwise and integral terms only in the above (f/, y)-system can be responsible 
for blow-up. We then need to consider two cases: 

Case I: nonlocal quadratic terms are negligible in (15. lip . Then, in the limit t T~, the 
^/-equation becomes asymptotically pointwise in main terms, and the system reads: 

'Ut + UUr + WU, = AU + l V^, 
Vt + UVr + WV, = AV~\ UV. 
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(5.13) 



It then follows by a standard comparison (barrier) arguments that the growth of all the 
localized solutions are controlled by the ODE system for sup ersolut ions, 



(5.14) 



* r •' ' dU _ _V 



which thus assumes no blow-up at all. 

Case II: V is not essential for blow-up. Then, performing a preliminary passage to the 
limit as r = Tfc + s — s> oo (r = — ln(r — t)), we arrive at the system with V = 0, which is 
known to admit no blow-up. 

Case III: quadratic nonlocal terms reinforce ^ V"^. Thus, this is the only possible case. 

However, a careful analysis of those two integral terms in (15. lip shows that both have 
the wrong sign to do that. Without pretending to any rigorous conclusions, we speculate 
about the sign of the integrals in (15.151) using standard clues, which look naive, but very 
often give correct answers in many RD-type problems. 

Indeed, integrating over a small neighbourhood of the maximum point in U (say, x = 0), 
where single-point blow-up occurs, these terms have the signs of the following integrals 
(others have a similar nature and can be checked out analogously): 

(5.15) ~ - t/,A,f/ and - AM^U.. 

Thus, we assume that the main sign-dominant part in these integrals is delivered by 
integration over sufficiently small neighbourhood of the maximum point. Then, since 
there Ur < (we think that the internal area, where Ur can be positive, is not dominated 
in the integral), ArU < 0, and Uz < 0, these integrals cannot be essentially positive on 
shrinking compact subsets to create an additional new type of "nonlocal" blow-up. In 
other words, the quadratic integral operators have the tendencies to be non-positive on 
such typical axi-symmetric flows and actually assure their extra stabilization. 

A difficult scrutinized analysis is necessary to fix such a non-blow-up conclusion. How- 
ever, this is a good sign for us, since we see that, for blow-up, a special (^-acceleration as 
t —>■ T~ is crucially needed to create extra positive quadratic integral terms to generate 
blow-up, where the TW dependence (11.251) to be introduced is the simplest opportunity. 

5.3. Positive integral quadratic (/^-dependent operators: a rout to blow-up. We 

briefly review such an opportunity. Thus, (15. 9p reads (we box Z^^-dependent operators) 



(5.16) W = m{rU)r + 

This adds into the right-hand side of the [/-equation (15. lip the following extra operators: 

Ut + uUr + lvu^ + wUz = Au + lv^-...- {\u(Jv^, 



(5.17) 



+ 



The l^-equation, due to the pressure term —^P^, also gets a new complicated operator, 

i VV^ + WVz = AV-}. UV 
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(5.1^ 



Vt + UVr 



^.n^w-Aw),. 



Therefore, plugging (15.161) yields several yj-dependent quadratic operators listed below: 

Vt + UVr + \ VV^ + WV, = AV-lUV- ... 



(5.19) 



Thus, the control of those new swirl-type operators in (I5.17P and (15.191) becomes much 
more difficult and this settles the core difficulty of the blow- up/non-blow-up problem 
for the NSEs. Now the (/^-derivatives are essentially involved into evolution, including 
even the third-order one in the integral term containing ^ VV^^^p. In particular, the 
simplest way to introduce the yj- dependence is the proposed blow-up swirl mechanism in 
(11.251) . which actually has a self-similar form and corresponds to periodic omega-limits. 
Thus, it seems plausible that some of the integrals in (I5.19P can be divergent enough 
to generate blow-up, unless all of them have the opposite negative sign, or their mutual 
interaction prevents singularities. Recall that each integration in z, according to the 
similarity variable y in (I2.18P yields the multiplier a/T — t, which is not powerful enough 
to compensate other singular scalings in (f/, V, r) there. Note again that, since U +oo, 
the regularizing anti-blow-up term ~ UV in the \^-equation must be defeated by some of 
the divergent nonlocal terms with a (y9-swirl. As in R-D blow-up theory, all such similarity 
of approximate self-similar balances eventually lead to complicated open problems for 
nonlocal systems, which later on we will attack and explain by using other, more local, 
blow-up approaches. 

On the other hand, for purposes of proving global solvability, the above {U, y)-system 
(I5.17p . (I5.19P (it is a simplified one) looks most reasonable. Then a careful step by step 
analysis of nonlinear mutual interaction of all the operators in (I5.19P is unavoidable. 

5.4. Generalized von Karman solutions on a partially invariant linear subspace: 
first term of general expansion. We look for velocity field and pressure of the form 

(5.20) U = U{r,ip,t), V = V{r,ip,t), W = zW{r,ip,t), and P = P{r,ip,t). 

This dependence models a class of dynamical stretched 3D vortex fiows including Burgers ' 
vortices (1948) [24j. The structure of such solutions also corresponds to the earlier (1921) 
classic von Kdrmdn swirling flow solutions of (II. ip exhibiting typical linear dependencies 
on two independent spatial variables x and y |110[ 1112] , 



(5.21) 



u = f'{z)x - g{z)y, 
V = f'{z)x + g{z)y, 
w = -2f{z), 
p=-2[f'{z) + p{z)], 



where functions / and g satisfy a system of two nonlinear ODEs, 



(5.22) 



g" + 2fg' - 2f'g = 0. 
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These solutions have been apphed to various problems of fluid dynamics; see Berker |13j . 
Note that, unlike von Karman similarity solutions fl5.21l) that induce ODEs fl5.22p . the 
reduction f l5.20p leads to a more complicated system of nonlinear PDEs, which anyway is 
simpler than the original model. Stationary solutions of the Euler equations (11.311) with 
a linear ^-dependence of w as in (I5.20p were also used by Oseen (1927) [TTlJ (cf. the 
structure (12.150 for Oseen's vortex); see more details in [87] . 

Actually, (I5.20p means that the vector field u and p belong to a 2D linear subspace, 

(5.23) u, peW2 = Span {1,2}. 

Then some of the equations (15.250 below (projections of (15. 3p onto W2) express the partial 
invariance of W2 with respect to the nonlinear operators in (15. 3p : see [76!, Ch. 7] for further 
details and examples. Here the invariance of an /-dimensional linear subspace Wi under 
a given nonlinear operator A is understood in the usual mapping sense, 

A{Wi) C Wi. 

Partial invariance means that only a subset of Wi satisfies this inclusion. Similar 2D 
restrictions on W2 = Span {1, y} exist for the Navier-Stokes equations in dimension N = 
2; see [761 P- 34] and references therein. The corresponding blow-up solutions are described 
in [791 Ch. 8]. 

It is key to note that the third velocity component W in (I5.20p becomes the first term 
in the asymptotic expansion of general solutions, for which 

(5.24) W = zW + E(,>2) ^'W,. 

Then, as usual in asymptotic expansion theory [107j . the first term is governed by the 
most nonlinear system (to be studied), while the rest of the coefficients solve "linearized" 
systems (that, in view of equations for other components, can be also rather difficult, at 
least first ones). In general, with a great luck, analytic expansions such as (15.240 might 
lead to solutions that are spatially localized in z; see further discussions below. 

Thus, substituting into yields 

(Ut + uUr + }yu^-lv^ = -Pr + A2U-^v^-j,u, 

Vt + UVr + ^VV^ + ^UV = -lP^ + A2V+^U^-^V, 
Wt + UWr + I VW^ + W^ = A2W, 
Ur + hU + hV^ + W = 0. 



(5.25) 



The last equation in (I5.25p easily determines U in terms of W and V (cf. (15.160 ). 

(5.26) U = -^ J zW{z, ip,t)dz-} J V^{z, if, t) dz. 
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As customary, the first equation in fl5.25p is a pressure one that a posteriori is going to 
define the pressure. Excluding P from first two equations yields the following PDE: 

[Ut + UUr + lVU^-',V'-A,U + ^V, + j,U]^ 
= [r {Vt + UVr + l VV^ + IUV-A2V-^U^ + ^ V)] ^. 

It is not difficult to see that this awkward equation is just a pseudo-parabolic PDE that 
also can be reduced to a more standard semilinear nonlocal form. Therefore, (15.271) em- 
phasizes the fact that the NSEs (11. ip admit a unique local semigroup of smooth solutions. 
Thus, taking the third equation, 

(5.28) Wt + UWr + i VW^ + W^ = A2W (with the constraint flOG]) ). 

we arrive at the system (15.271) . (I5.28P for two unknowns W and V. For convenience, we 
now return to the original system (I5.25P for performing further blow-up scaling. 

5.5. Blow-up twistor variables and rescaled equations. We introduce next the fol- 
lowing blow-up rescaled independent variables, where, for convenience, the blow-up time 
is reduced to T = 0: 

(5.29) U=-^u, V = -^^v, W = j^,w, and P = j^,p. 



-["'1 " — ^' " — {-t) "^^^ ~ i-t) 

The rescaled dependent variables are given by: 

(5.30) y = (f = fi — a \n{—t) , and r = — ln(— t), 

where a 7^ is an unknown parameter, which cannot be scaled out and plays a role of a 
nonlinear eigenvalue for the future "stationary" problem. As usual in blow-up problems, 
the new time variable gets infinite, 

(5.31) T = - In(-t) +00 as t ^ 0". 

It is crucial that, according to (I5.30p . on compact intervals in the rescaled angle n, 

(5.32) {p = fi — a \ia{—t) = fi + a r 00 as r — > -|-oo (i.e., t — 0~). 

Therefore, in the rescaled angle variable, the dependence ip = fi + a t refiects a travelling 
wave angular behaviour, where a is the wave angular speed. The angular dependence in 
(15.320 corresponds to infinite acceleration of rotation of all the velocity components about 
the 2-axis, and therefore we call special pattern solutions (I5.29p . associated with scaling 
(I5.30p . blow-up twistors. 

Remark 1: on non blow-up flows with spiral symmetry. Such solutions of the 
NSEs (ll.ip in the cylindrical coordinates with the usual TW-type angular dependence 

(5.33) (p = fi + at 

are well known and were studied by Bytev [25J in 1972; see also the results of the group 
classification of the NSEs in [2J. Unlike (15.290 . the standard invariant of the group of 
translations associated with (I5.33P does not allow blow-up. Indeed, logTWs and simply 
TWs are different group-invariant solutions, and the former ones assume extra invariance 
relative to a group of scalings. 
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Remark 2: on tornado-type blow-up in complex NSEs. Tornado-type blow-up 
behaviour of the complex version of the NSEs, for which blow-up was established in Li- 
Sinai [138] ■ was confirmed in |5j by numerical methods (though the achieved numerical 
evidence concerning the structure of such singularities still looks rather unsufficient). 

In view of (11 ■4p . let us now check the nature of the scalings in (I5.29p . (15.301) . According 
to full similarity rescaling in (12. 7p . the variable is 

(5.34) C = ^^ =^ ^W = ^^Cw, 

and hence all these scaling factors are self-similar. However, (II. 4p is not a warning for 
us, since logarithmic angular TWs do not belong to the framework of standard similar- 
ity patterns (in fact, these are related to periodic orbits). In addition, we are going to 
use nonstationary evolution governed by the PDEs (15.350 including patterns without a 
similarity stabilization. On the other hand, for axisymmetric flows, it is known that the 
blow-up must be of Type II, i.e., with the blow-up rate faster than self-similar [llSj [38]. 

Substituting (I5.29p . (15.301) into the general system (I5.25P yields the following nonsta- 
tionary rescaled equations: 



(5.35) 



,,2 Un 2 V, 



' Ur + \ yuy + \ u + au^ + uuy + ^ vu^ - ^v"^ = -py + A2U -^v^- 

Vr + \yVy + \v + aVf, + UVy + ^VV^ + ^UV = -'^P^ + A2V + 

Wr + I ywy + W + aw^ + uwy + ^ vw^ + w'^ = A2W, 
Uy + - u + - Vi, + w = 0. 



The pressure and the constraint (15.261) interpretation remain the same as for (15.251) . so 
(I5.35P is a nonstationary nonlocal system of two equations for unknowns w and v. Recall 
that the 2-Laplacian A2 now takes the full form 

(5.36) A2 = D^^ + ];Dy + ^Dl. 

5.6. Remark: log-TWs in reaction-diffusion systems. Using blow-up logarithmic 
travelling waves comes as a fruitful idea from reaction-diffusion theory; see examples in 
|187l pp. 105, 308, 411]. For instance, the classic quadratic porous medium equations with 
source and convection in ID (a canonical combustion problem with regional blow-up, |187[ 
Ch. 4]) 

admits the following blow-up logarithmic travelling waves: 

uix,t) = j^f{y), y = x-a\n{-t) =^ f + a f = {f f')' + f f + p. 

The scaling group-invariant nature of such logTWs seems was first obtained by Ovsian- 
nikov in 1959 [172j , who performed a full group classification of the nonlinear heat equation 



Ut = {k{u)ux)x, 
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for arbitrary functions k{u). In particular, such invariant solutions appear for the porous 
medium and fast diffusion equations for k{u) = m", n ^ 0: 

nt = KnJ. =^ 3 uix,t)=t-v.f{x + a\nt), where - i / + af = (/"/')'• 

Blow-up angular dependence such as in (15.321) was studied later on in [[7j, where the 
corresponding similarity solutions for the react ion- diffusion equation with source 

(5.37) Ut = V ■ (u'^Vu) + u'^ in M=^x(0,T) (/5 > 1, a > 0), 

were indicated by reducing the PDE to a quasilinear elliptic problem (it seems, there is no 
still a rigorous proof of existence of such patterns). For parabolic models such as fl5.37p . 
that are order-preserving via the MP and do not have a natural vorticity mechanism, 
such "spiral waves" as t — T~ must be generated by large enough initial data specially 
"rotationally" distributed in M^. For the fluid model (II. ip with typical vorticity features, 
this can be different. 



5.7. On some auxiliary properties of twistor structures. As usual, these are sta- 
tionary solutions of (15.351) . which are independent of the time- variable r, i.e., in the 
original variables. 



(5.3^ 



U{r,(p,t) = -^u(?/,^), V{r,ip,t) = {)(?/, /i), 
W{r,(p,t) = ^w{y,fj,), P{r,ip,t) = jz^p{y,fi). 
The rescaled profiles u, v, w, and p solve the corresponding stationary system, 

^yVy + lv + av^ + uVy + ^vv^ + ^uv = -^Pf, + A2V + ^ - ^ i), 
i yWy + W + aWfj. + uwy + ^ vw^j, + = A2t&, 

Uv + 7.U+7.Vu + W = Q. 



(5.39) 



y ■ y ^ 

We now present the full system for w, v. First, we perform the reflection 
(5.40) w t— > —w. 

Second, we have from the last equation in (15.391) 



(5.41) 

Then the i&-equation reads 



u 



- J zw dz — - f dz. 



A2W — I yWy — w + w'^ 



(5.42) 



j ZWdz)Wy+\{ j V^dz)Wy-Hw^ = {). 
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Thirdly, the (m, ■0)-equation (with fl5.4ip ) takes the form 
(5.43) 

= [y{iy^y + l^ + ov^ ^uvy^\ vv^ + ^uv - A2V - ^u^ + ^v)]y. 

Thus, we arrive at the system of two PDEs fl5.42p . fl5.43p . with the nonlocal constraint 
(15.411) . which, as we have seen, actually means the presence of an extra first-order PDE. 

This system with the real parameter a 7^ is indeed complicated. We note that the 
Laplacian A2 in (I5.36P includes second-order derivatives in both radial y and the angular 
yU variables. The latter dependence is crucial for such similarity twistors. We do not 
intent and do not plan to study this system somehow rigorously, though some results, in 
particular, about asymptotic distributions a.s y ^ +00 can be obtained and will be of 
importance in what follows. 

Obviously, in view of (I5.20p . these similarity twistors as special kind of blow-up vortices 
are not spatially localized in the 2;-direction, so they have infinite energy always, and, for 
their existence, a whole M^-space should be taken into account (we skip at this moment 
their expansion meaning (I5.24p ). But as typical in reaction-diffusion equations such as 
fll.36p . (I1.54p . fll.56p . and many others, such similarity or approximate similarity blow-up 
structures can appear from local finite and well-spatially-localized data. Here there appear 
hard questions of their evolution stability (in which intermediate sense?) or behaviour 
close to centre and/or stable manifolds associated with their partial similarity space-time 
geometry. These questions are addressed to the full non-stationary system (I5.35p . 

As we have mentioned, in the case of nonexistence of such "stationary" swirl patterns, 
some time-dependent perturbations of these structures can lead to various non self-similar 
blow-up twistors that evolve close to invariant manifolds associated with the scaling (I5.38p . 
Actually, precisely this always happens to the semilinear scalar curvature equation fll.60p 
for < 6 and for the Frank-Kamenetskii equation (I1.36P for A^ = 1 or 2. Unfortunately, 
this will lead to much more difficult mathematics than for a simpler purely stationary 
similarity one (that does not look easy at all, of course). 

5.8. On extension beyond blow-up for t > T. This is connected with Leray's blow- 
up scenario (II. 5p . We do not discuss such difficult issues here in any detail and just 
mention that, as typical, the extension after blow-up is a secondary question relative to 
the blow-up ones on the behaviour as t ^ T~ , which still remains mysterious. 

Let us also note that, concerning blow-up twistors with the angular variable (15.320 . a 
natural way of extensions of possible blow-up solutions is to use the forward variable 

(5.44) (f = fi + a\nt for t>0 

and — t I— s> t in (I5.29p . where a ^ a, in general (but signer = signer to keep the direction of 
rotation), is an extra parameter of matching at t = T = 0. Principles of proper matching 
of a blow-up flow for t < with a more regular one for t > for (ll.ip are not clear. For 
simpler nonlinear PDEs such similarity matching is possible; see [TOj for (ll.57p . 
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5.9. Discussion: on auxiliary properties of blow-up twistors. Let us look for some 
auxiliary properties of the "stationary" profiles u, v, w as solutions of fl5.4ip - fl5.43l) . To 
reveal a first key feature, assume for a moment that these profiles are of changing sign, 
so that rescaling f l5.29p would mean non-uniform blow-up, where components diverge as 
t — ^ 0~ to ±cx) on some subsets. 

We next need to point out another curious possible feature of such solutions that is of 
importance for general understanding of the nature of such singularities. Let us demand 
that the similarity profiles have angular periodic behaviour with zero mean as ?/ — >■ oo. 
For instance, assume for the component w{y,fi) that 

(5.45) t&(?/, /i) = + ... as y ^ +00, where jC(/i)d/i = 0. 

^ 

Then, in view of the behaviour of the rescaled variables in fl5.30p near blow-up time, under 
some extra hypotheses, this can imply, by the blowing up rotational behaviour, that 

(5.46) W{r,ifi,t)^0 as t ^ Q- , 

in the weak sense on the corresponding subsets in the rescaled variables. Indeed, in view of 
the zero mean in fl5.45p . the equality (15.460 is then a manifestation of Riemann's Lemma 
from Fourier transform theory. As more usual and typical, under weaker assumptions, 
such solutions locally represent an oscillating "building block" (see [T7]) for x ~ 0, t — »• 0~ 
in terms of the corresponding vorticity = V x u, where u is then defined by the Biot- 
Savart law, 

(5.47) u(a;, t) = J j§E^ x uj{y, t) dy in ^ 

The property fl5.46p (and also relative to the vorticity, which is easier) then represents 
special type of singularities that are almost "invisible" (in fact, efficiently "nonexistent") 
close to the blow-up time in the natural integral (weak) sense. Recall that this could be a 
key feature since all the differential equations of fluid mechanics are derived from kinetic 
equations with integral operators in collision-like terms by approximating typical integral 
kernels involved by kernels with pointwise supports (Grad's method in Chapman-Enskog 
expansions). Therefore, the integral, average meaning of coherent macro-structures and 
micro-singularities are of practical importance only. 

Another principal question is as follows: can, in the non-stationary setting, the "zero 
mean oscillation property" such as fl5.46p in a neighbourhood of the origin affect and 
diminish the asymptotic behaviour of the corresponding solutions as y ^ 00 to get them 
into finite energy class? It seems that this is not that essential for the present cylindrical 
geometry, but can be key for the spherical one with the unique single point at the origin; 
see Section |Hi 

5.10. Discussion: on asymptotic stability of twistors. We begin with the following 
first observation concerning evolutionary significance of similarity scaling: 

Assume, not taken into account possible energy characteristics, that a suitable non- 
trivial solution of the stationary system fl5.42p . fl5.43p exists. Then the next and a more 
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difficult step is to study its asymptotic stability in the framework of the full non-stationary 
system fl5.35p . This stabilization phenomenon then becomes a difficult open problem. 
Observe that there is no chance to prove that f l5.35p is a gradient dynamical system in 
any admissible metric. Recall that the negative result in [103t Th. 1, § 2] establishing the 
self-similar ban is proved under the a priori assumption on convergence to a fixed rescaled 
profile (then necessarily this proffie must be a stationary solution and hence it is zero by 
|162 ]): and finally, the exponential decay via spectral and semigroup characteristics of the 
rescaled infinitesimal generator B* in L^, (R^) implies that the solution u(x, t) is bounded 
at t = 0". 

In this connection, it is worth mentioning that, as known from reaction- diffusion the- 
ory, regardless bad energy and other global properties of twistor profiles, good rescaled 
solutions can stabilize to them in a local topology, e.g., in view of interior regularity, 

(5.48) uniformly on compact subsets 

in the rescale variables. Of course, uniform stabilization in M.^ is then impossible in view 
of the energy discrepancy. 

5.11. Discussion: on linearized blow-up patterns about a constant equilibrium, 

(I). Assume next that either (i) any suitable non-trivial steady profiles satisfying fl5.42p . 
fl5.43p are nonexistent, or (ii) there exists a non-trivial solution (both possibilities are 
suitable for us). What kind of other behaviour do we then expect of the nonstationary 
system (15.350 ? In other words, is there any hope to get a kind of entirely non-self-similar 
blow-up twistor behaviour? 

The idea of such "linearized construction" is as follows. We ffist need to fix a family 
of simpler local solutions or "almost" (say, slightly perturbed) solutions of the dynamical 
system under consideration. As the next step, we linearize the flow about this family and 
use the orbits on the corresponding stable or centre manifolds to match them with the 
surrounding orbits of the necessary behaviour and regularity. These matching procedures, 
though being rigorously justifies for some simpler parabolic PDEs such as (I1.36P or (ll.60p 
(we have presented the references), remain open for many other quasilinear and higher- 
order parabolic reaction-diffusion models, and, surely, this will be the case for (II. ip . 

We show how then a countable family of the so-called linearized blow-up patterns (un- 
like the above similarity patterns that are nonlinear eigenfunctions) can be constructed. 
Consider the nonstationary t/;-equation from (15.350 bearing in mind the reflection (l5.4Up . 

(5.49) Wr = A2W — 1 ywy — w + — uwy — aw^ — ^ vw^. 
We see that (I5.49p admits the constant equilibrium 

(5.50) = 1, i.e., iy(r,v3,t) = p^, 

which corresponds to uniform global blow-up as t — > 0^ in the whole space. This is 
the easiest exact solution, which we are going to linearize the flow about. It is known 
from reaction-diffusion theory (the proof is obvious for systems with the MP) that in 
order to blow-up at the fixed time t = 0, the rescaled solution w{y, r) should always 
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be "sufficiently" close (in a certain metric) to the constant blow-up profile (15.501) . since 
otherwise, essentially smaller solutions must exhibit larger blow-up times, so these become 
exponentially small as r — > -\-oo after rescaling (15.291) : see a spectral justification later on. 
Thus, we perform the standard linearization about (15.501) by setting 

(5.51) w = l + Y 

to get the "linearized" equations (we have taken into account the ^-representation (I5.4ip ): 

Yr = ^2Y -yYy + Y- oY, -\vY,^\ {j,'v,)Yy + B{Y), 
^^•^^^ where B{Y) = Y^ - I [J^' zY)Yy, 

is a quadratic perturbation SisY—>-0. Consider the radial part of the linear operatoi0 in 
(15.521) excluding at this moment its tangential angular first- and second-order operators 
including that in A2, 

(5.53) B* = A2 - yDy + I. 
Changing the radial variable yields 

(5.54) y^f^ =^ B* = 2(A2-iyZ^, + i/), 
where we observe the adjoint Hermite operator (cf. (I1.70p ) 

(5.55) B* = A2-^yDy, with B = A2 + lyDy + I in L'^{R^). 

On the other hand, B is self-adjoint in the weighted space L^(R^), with the weight 

(5.56) p{y) = e 4 , 

so that this falls into the scope of classic theory of linear self-adjoint operators, [T5] . 
In particular, the adjoint operator B* in the adjoint space L^*(M^), where 

(5.57) p*(^y) = ^ = e-4, 

has the discrete spectrum (here we take into account the restricted radial symmetry) 

(5.58) a{B*) = {\2k = -k, A; = 0, 1, 2, ...}, 

and the eigenfunctions are normalized Hermite polynomials [T5| p. 48] 

(5.59) i)lj^{y) = C2kH2k{y), = fc = 0,1,2,.... 

As was already pointed out, in the general non-radial setting in M^, all these polynomials 
are obtained by differentiating the rescaled Gaussian F[y) of the fundamental solution of 
the heat operator Df — A^, 

(5.60) Fiy) = j^e-^ (iV = 2), 
so that the following generating formula holds: 

(5.61) ^*(y) = ^cpD^F{y)^cpHp{y), where = ^ 

■^^The non-radial eigenvalue problem with a ^ is studied in Section [6l 
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and P = (/3i, ...,Pn), with = Pi + ... + Pn, is a multiindex. This set of polynomials 
include all the angular-dependent eigenfunctions, so these are complete and closed in the 
whole weighted space L^*(]R^). Note that, unlike Section [221 "we do not need solenoidal 
test for eigenfunctions involved. 
It then follows from (15.541) that 

(5.62) a(B*) = {X2k = l-2k, k = 0,1,2, ...}. 

The first unstable mode with k = and Aq = 1 corresponds to the unavoidable unstability 
of the scaling (I5.38P with respect to small perturbation of the blow-up time, when we 
replace ^— T ^ 0. On the manifold of solutions with the same blow-up time, this 
unstable mode plays no role and is excluded. 

Thus, the first actual mode takes place for k = 1, with 

(5.63) A2 = -l and V^^(i/) = C2 (l - i y^). 

This corresponds to patterns with the following behaviour for large r ^ 1 in the inner 
region that is characterized by arbitrarily large compact subsets in y: 

(5.64) W2iT) ^l + e~^C^;{y) + ... {C > 0), 

where the tangential operator in (15.521) should be now also taken into account. This 
makes the spectral theory more involved (other equations are also included) and can lead 
to complicated computations. For the purely radial case, analogous computations will be 
performed in the next section. 

Meantime, we deal with the conventional expansion (15.641) . which should be matched 
with the outer behaviour for y ^ 1. The corresponding rescaled variable in this outer 
region is also seen from (15.641) by using the quadratic behaviour of H2{y) as y — > oo. 
Indeed, this gives the first term of the stationary expansion in the new rescaled variable: 

(5.65) W2{t) ~ 1 — C2C e~'' y'^ + ... ~ 1 — C2CC^ + ... for y > 1, where C, = e^^^y. 

In a similar manner, we find other patterns with the inner behaviour governed by other 
stable ID eigenspaces of B*, where 

(5.66) W2kir)^l + e^^''^Cr2k{y) + -, where A; = 1, 2, 3, ... , 

and ip2kiy) other higher-degree non-monotone Hermite polynomials. In general, this 
can lead to a countable set of different blow-up patterns; see details in [70] . 

For the non-radial case, where all the tangential operators should be taken into account, 
the matching procedure with the outer region can be very complicated and cannot be jus- 
tified rigorously for the full model. This is done in |7D] for simpler blow-up axisymmetric 
jets that are admitted by the Navier-Stokes equations in M^. 

In general, on matching, we arrive at a special countable family of blow-up patterns 
that are characterized by a stronger and wider propagation than the (existing or not) 
self-similarity scaling (I5.29P suggests. For instance, for the scalar curvature equation 
(11.601) . the first generic blow-up pattern (ll.62p . which is characterized by a slow drift on 
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a centre manifold corresponding to the similarity rescaled variables (11.641) . shows a faster 
propagation in the x-direction with the extra logarithmic factor (here T = 0) 



~ = \/| ln(-t)| +00 as t ^0 . 

Such an expanding blow-up twistor can be not that adequate to the nature of single-point 
blow-up for the Navier-Stokes equations (at least, what we could expect). 
Therefore, we begin discussing the second, probably, more reahstic situation. 

5.12. Discussion: on linearized patterns about singular equilibria, (II). We now 
will use other equilibria of the equations (11. ip . which are not constant and hence more 
and better localized about the stagnation point x = 0. 

Singular radial equilibria. The first candidate is singular stationary solutions that exhibit 
strong singularities at x = 0, and are possibly not any weak, very weak, or mild solutions 
at all. This is not that important for us, since we are going to use them just for linearizing 
and next remove the singularity by matching with a regular bundle at the singularity point. 

As an easy illustration of the type of calculus to be performed later on, let us look for 
simple homogenuity —1 separable stationary solutions of the original system (I5.25p . 

(5.67) f/(r,y,) = ^, V{r,^) = ^, w{r,^) = ^, P(r,^) ^ ^ 
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These separable structures are invariant under the scaling group in (I5.29p . so that the same 
singular patterns in the variables (15.300 can be used for the study of twistor behaviour in 
the rescaled system ( 15.35p . 

Substituting (15.671) into (I5.25P yields the following system on the unit circle in M?: 



(5.68) 



A' - BA - 2B' + ^2 + 52 + 2D = 0, 
B" - BB' + 2A' - D' = 0, 
C" - BC + 4^ + 2AC - = 0, 
5' + C = 0. 



The second equation is integrated once that, on using the last one, gives D as a quadratic 
function of other unknowns, 

(5.69) D = 2A-\B'^ - C + Mq (Mq = const.) 
On substitution into (I5.68p . we obtain an easier system 

( A" - BA' + A^ + AA + 2Mo = 0, 

(5.70) \ C" - BC + 4C + 2AC - = 0, 

[ 5' + C = 0. 

Substituting C = —B' from the last equation yields two equations for functions A and B, 

(A"-BA' + A'' + iA + 2Mo = 0, 
\b"' - BB" + {B'f + 2{A + 2)B' = 0, 
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which still remain a difficult fifth-order ODE on with an arbitrary parameter Mq G M, 
which can play a role of a nonlinear eigenvalue. In particular, it is easy to check that 
there are no explicit trigonometric solutions, where 

A,B eW3 = Span {1, cos /v?, cos/y?} [l e M). 
Finally, we note that there exists simpler "irrotational" (no swirl) solutions (15.671) : 

(5.72) C = and + + 2D = 0, 

about which first linearized analysis can be naturally began with. 

It is worth mentioning here that Sverak [204] proved that all singular stationary ho- 
mogenuity —1 equilibria of the NSEs in \ {0} are equivalent, up to orthogonal trans- 
formations, to Slezkin-Landau's solutions; a stronger result was obtained in [155] . see 
Section [71 Since singularities of (I5.67P are concentrated at the z-axis {r = 0}, we cannot 
use the result of [204j . though it does give a hope to avoid to scrutinize the general system 
(I5.70|) . which looks not that easy at all. In addition, there is a strong nonexistence result 
in [118 , § 5] for axi-symmetric ancient solution^ of the NSEs: 

(5.73) \u(x,t)\ < , 'r , in M3x(-oo,0) =^ u = 0. 

It seems this does not directly apply to equilibria f l5.67p and similar others, since these 
are singular at r = (i.e., are not bounded ancient solutions). Anyway, it seems that 
some singular equilibria can be eventually ruled out but not all of them. 

Thus, by U we denote a certain singular equilibrium of the NSEs, not necessarily given 
by the homogeneous formulae fl5.67|) . In general, the description of all the possible singular 
stationary orbits at x = (or at r = 0) leads to very difficult study of ill-posed elliptic 
evolution equations as explained in Section 18. 1[ 

Construction of blow-up patterns: in need of non self-adjoint theory. Thus, we are not 
going to study the system (I5.7U|) in detail here, and will concentrate on the principles 
of construction of blow-up twistors using linearization about the manifold corresponding 
to such singularities. Denoting by u the 3- vector of the rescaled variables u = {u,v,w), 
we write ( I5.35p as a dynamical system for the rescaled variable (cf. ( 12.211) ) 

(5.74) u' = H(u) for r > 0, 

where, as in fl2.21|) or fl2.23|) . the pressure variable is excluded by projecting the solution 
space onto the kernel of the gradient operator. Thus, we have fixed U(|/,/x) ^ as a 
solenoidal stationary singular solutions of the rescaled equation fl5.74p . 
We next perform the linearization by setting 

(5.75) u(r)=U + Y(r) 

"'^I.e., defined for all t < 0, [Ml. By scaling (|2.6p . any L°° blow-up solution u{x,t) generates a non- 
trivial ancient uniformly bounded one v{x,t) = ^ u(xfc + -§^,tk + -^), where Ck = sup^ |u(a;,tfe)| = 
|u(xfc,tfc)| +00 (on the applications, see [118j ). This reflects typical scaling tools of R~D theory [79| : 
cf. Section [H and applications to solvability and bounds for nonlinear parabolic equations in [80 , § 2] . 
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to get the linear equation for Y, 

(5.76) Y = H'(U)Y + D(Y) for r > 0, 

where D denotes a nonhnear perturbation, which is quadratic as Y ^ 0. 

It follows from the general steady structure such as fl5.67p that the first step of such 
a construction is to find the point spectrum of the linear integral (pseudo-differential) 
matrix operator B* = H'(U) for vector- valued functions with periodic in /i (and possibly 
singular at the origin) coefficients. Thus, we define the Inner Region / as a family of 
special space-time subsets, where the linearized equation holds asymptotically: 

(5.77) Inner Region I: Y = H'(U)Y + ... for r > 1. 

Recall that, by construction, H'(U) is assumed to act in a solenoidal vector field. 
A number of accompanying questions arise. E.g., one can pose the following: 

Question (i).- B* = H'(U) is not self-adjoint in no weighted L^-spaces; 

Question (ii).- what is the domain of B* (hopefully a kind of Sobolev space if^, (M^))? 

Question (in): what is the space and domain of the adjoint operator B? 

Question (iv).- it is not clear why B* and B could have enough or at least some real 
(or complex) eigenvalues, compact resolvent (by compact embedding H^* C -^^p.? which 
conditions at the singular origin?), and bi-orthogonality property of bases (if any); 

Question (v).- why the eigenf unctions of B* can be at least approximately associated with 
some structures of finite polynomials (possibly separable in the angular direction)? 

Question (vi) on solenoidal eigenf unctions and classes. Though this is a natural part of 
the eigenvalue problem for the pair {B*, B}, let us recall the key: assuming the — 1 
homogenuity (and if not?) of the equilibrium such as fl5.67p employed, we are obhged to 
perform a "solinoidalization" to get a sufficient amount of (possibly even some polynomial 
where again the angle separation is assumed?) eigenfunctions ip'^ , which will generated 
eigenspaces of solenoidal fields as in f l2.3ip : ... , etc. 

For the singular S-L solutions (17. ip . this will be continued in Section [71 where we also 
discuss in greater detail the questions of matching of various asymptotic inner and outer 
regions. 

Nevertheless, the parabolic experience of doing linearized theory for 2mth-order PDEs 
such as f ll.56p shows that the corresponding linear operators B and B* are also non self- 
adjoint (the only self-adjoint case is for m = 1; cf. fl2.24p ). but admit real point spectrum 
only, [51j. However, the matrix spectrum problem for the pseudo-differential operator B* 
in (15.761) and for the corresponding formally (in the topology of L^) adjoint operator B 
are much more difficult and remain open; see Section [7] for further details and related 
important comments in [204] . 

Thus, as for an illustration, assume that B* has a real point spectrum 

(5.78) a{B*) = {\k, k=\p\ = 0,1,2,...} (/3 is a multiindex), 
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where the first positive eigenvalue, 

(5.79) Ao > ( = 1, it seems), 

reflects the natural unstability with respect to the change of the blow-up time. We also 
assume that, as a standard fact from linear operator theory (that also needs a difficult 
approving), the eigenfunctions {v^} are assumed to be bi-orthonormal to the correspond- 
ing adjoint basis {v^}, and both bases are complete, closed, and form Riesz-type bases in 
the corresponding solenoidal spaces. 

We are interested in eigenvalues A^ < 0. The most interesting case occurs when A^ = 0, 
since it corresponds to a behaviour close to the centre subspace of B*. Note that proving 
existence of a centre or stable invariant manifolds for the full nonlinear problem (15.741) is a 
very difficult problem. As a formal hint, assuming that the kernel of B* is one-dimensional 
spanned by an eigenfunctions (for any dimension, the analysis is similar with a system 
to occur), we are looking for solutions in the form0 

(5.80) Yo(r) ~ C/3(r)v; + w(r) for r > 1 (v^ G cv* = {civl, Cat;^, CsvD), 

where the remainder w(r) is supposed to be orthogonal to kerB* in the dual metric 
between the spaces for B and B* denoted by (■, ■). Substituting into (15.761) and multiplying 
by yields the following asymptotic ODE system: 

(5.81) Cfs = (D(c^v;), v^) + ... = Q^(c^) + ... , 

where Q/3(c^) is a vector quadratic form in M.^. Looking for a standard solution with a 
power decay yields a quadratic algebraic system on the coefficients, 

(5.82) c^(r) = ^ + ... =^ Q^(A^) = -A^. 

Assuming that the matrix equation in (I5.82p has a solution A^ ^ 0, we obtain the 
following asymptotic pattern: 

(5.83) Yo(y,r)~^v;(y) + ... as r ^ +oo (A^ = 0). 

In this way, we observe a slow logarithmic drift, with r = — ln(— t), relative to the centre 
subspace of H'(U). 

For Afc < 0, with k = we arrive at more standard stable subspace patterns 

(5.84) Yfc(r) ~e^'=-C^v; + ... (A^ < 0). 

Both (15.801) and (I5.84p describe for r 1 asymptotically small perturbations in (15.751) of 
a steady singular solenoidal field U, so we have to choose solenoidal eigenfunctions. 

On outer matching: Outer Region. The main idea of matching of these patterns with the 
outer region, which is most remote is explained in [TOj, where we are assuming that v^(?/) 
have a polynomial decay, for some coefficients 6k > 0, k = \(3\, 

(5.85) \^r*Jy)\^bky'^'' + ... as y +oo. 



'Recall (|2.3ip for such eigcnfunction expansions. 
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Here we arrive at a simpler first-order matrix operator, since in the new rescaled variables 
of the outer region, all the Laplacians form asymptotically small singular perturbations 
simply meaning that Euler's equations are dominant here. Of course, the resulting singu- 
larly perturbed dynamical system is difficult to tackle to pass to the limit r — >• +00. 

Note the following crucial property of such an extension into the next outer region 
II, where solution is much less: we are not now obliged to keep there the self-similar 
power-like decay at infinity. Moreover, we now can match this blow-up behaviour with 
the exponentially decay associated with the heat semigroup e^^"^ for the equations such as 
(15.491) and others. This matching is not a key difficulty of the blow-up twistor construction 
and deals with not that singular orbits. Again, we take into account patterns (I2.3ip . 

On principles of singular inner matching: Inner Region II. The last crucial step is to 
match the centre (I5.83P and stable (I5.84p subspace patterns with the regular and bounded 
solenoidal flow in the singular inner region close to the origin y = 0. This is the most 
involved matter, where the divergence-free asymptotics satisfying something like (12.311) 
become key. The singular inner matching in a countable number of various cases is re- 
sponsible for existence or nonexistence of solenoidal blow-up patterns and is the main open 
problem. 

Assuming for deflniteness that such matching exists, we predict some other evolution 
features of the resulting blow-up patterns. Here we formally suppose that the presented 
asymptotic expansions go along the solenoidal classes like in (I2.3ip . so that such spatial- 
temporal structures can be used for some preliminary estimates of blow-up evolution. 

Assume flrst that in the stable subspace representation (I5.84p . the eigenfunction is 
singular at ?/ = with the following behaviour (7^ > 1): 

(5.86) |v;(y)| ~4r^^ + ... as y^O. 

Therefore, the pattern (I5.84p has the asymptotic behaviour in the singular inner region 
of the type (we omit here non-essential constants and neglect other minor multipliers) 

(5.87) Ufc(y, r) ~ f - e^'^-Cpy-'^ + ... . 

Calculating as in (I1.46P (see [47j) the maximum point in y of the function in fl5.87p yields 

(5.88) supy>o |ufc(2/,T)| ~ e ^fc-^ ^+00 as r +00 (A^ < 0). 

Under the above assumptions, this is about a right bound on regular hypothetical con- 
tinuation of the pattern from an exponentially small neighbourhood of the singular 
state U to smooth bounded solutions ioi y ^ Note that the behaviour of the pattern 
as 1/ — i> can be very complicated including blow-up swirl-like features. This gives a 
preliminary estimate of the rate of this non self-similar blow-up in the u-variable 

(5.89) ||ufc(-,t)|U ~ (-t)~^^^ as t^O- (A; = |/?| > 1). 

Similarly, for the centre subspace pattern fl5.83p (if any), this very rough estimate imply 
the following perturbation of the self-similar rate: 

(5.90) ||uo(-,t)|U ~ ^|ln(-t)|^ as t ^ 0". 
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See Section 17.61 for further necessary properties of such a matching. 

Remark: back to reaction-diffusion. For single parabohc equations such as (]1.36p . 
fll.601) (for > 6), or a more general combustion model 

(5.91) ut = Au + \u\P'^u in x M+ (p > 1), 

such examples of creating non-similarity asymptotics by matched asymptotic expansion 
techniques are known and have been rigorously justified; see [HI EH [721 [771 EH], and 
references therein. Even for these reasonably simple parabolic PDEs, in the supercrit- 
ical Sobolev range as in (11.341) . the set of singular stationary solutions can be rather 
comphcated that affect the evolution properties of blow-up and other asymptotics; see 
|146l 11561 1157] and Section 17.101 For solutions of (15.911) of changing sign (that are always 
taking place in our case), the blow-up patterns for p = ps become more involved and their 
actual existence is still not fully justified rigorously, [59] . 

The reaction-diffusion equation (15.911) does not assume any conservation of energy, 
but, however, there are examples of the so called incomplete blow-up, when the solution 
u{x,t) becomes unbounded in L°°(]R^) at a single moment t = T only; see [77] and 
earlier references therein concerning equation (I1.36P , as well as |157l 1158] for more recent 
extensions and achievements. For the model (II. ip . blow-up moments are incomplete, in 
view of the a priori bound (12. 2p . 

It seems this quite a bounded mathematical experience cannot be directly translated 
to the Navier-Stokes rescaled equations (I5.74p . 

We have discussed the second principle to generate another countable set of blow-up 
twistors. No doubt that there are other ways to get various countable families of such 
patterns, which we are still not aware of and cannot even imagine. 

5.13. Discussion: on a very formal way to create a vertex and branching to lo- 
calized smooth blow-up twistors; a twistor ring. Here we discuss another important 
issue concerning the twistors on the 2D subspace Span {1, z}. 

(i) A vertex is by truncation of the VT-component in (15.201) via the positive part, 

(5.92) W = {z)+W, 

so that = for ^ < 0. This leads to a weak solution of (II. ip . so a certain effort is 
necessary to check whether this gives a Leray-Hopf solution in the sense of the inequality 
in (12. 2p for all t < 0, i.e., before the blow-up occurs. On the other hand, in view of the ex- 
istence of a local semigroup of smooth bounded solutions, existence of such solutions with 
a vertex assumes performing a certain smoothing at the stagnation point at a; = 0. This 
cannot be done in the above self-similar manner, so the proof of such a "smooth branch- 
ing" (if any) from the given non-smooth blow-up twistor represents a typical difficult open 
problem of modern theory of nonlinear PDEs. 

(ii) Similarly, there occurs another, not simpler, open problem of "smooth branching" 
from this unbounded twistor of other, partially and asymptotically self-similar blow-up 
twistors that might be sufficiently spatially localized as t — 0~. As we have mentioned, 
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a possible way to check a possibility of such branching consists of using asymptotic ex- 
pansion theory to create solutions via series like fl5.24p . Obviously, this leads to difficult 
nonlinear (first ones) and further linearized systems for the expansion coefficients and 
represents another face of this open branching problem. 

(iii) There is another formal geometric way to create blow-up patterns by deforming the 
axis of their symmetry. For instance, let us assume, in view of our pretty local analysis 
of such blow-up behaviour for small 2; ~ 0, that the 2;-axis can be bend into a finite ring 
such that at some points similar blow-up structures oriented in the tangential 2;-direction 
occur with a certain periodic behaviour along the ring generatrixlf^. This leads to "twistor 
rings" with a very difficult open accompanying mathematics. 

5.14. On twistors in Euler equations. Similar blow-up twistors, with analogous con- 
cepts of their swirl-like and other extensions, can be formally detected for the Euler 
equations fll.3ip . where the calculations are supposed to be simpler but actually they are 
not. Without the "curse of local smoothness" for the Navier-Stokes equations (11. ip . for 
the Euler ones (ll.3ip containing first-order differential operators only without a typical 
interior regularity, the necessary branching phenomena are easier to justify, so there is a 
hope of doing these more rigorously. However, one can see that, neglecting Laplacians A2 
in the systems (I5.25P and (15.390 . indeed simplifies the analysis but the remaining PDEs 
are still very difficult to understand rigorously that demands new concepts of solutions 
and entropy regularity. The stationary system (I5.4ip - (l5.43p with no Laplacian operators 
A2's takes the form 

' — I yWy — w — aWfj_ + w"^ 

-\ ( /cf dz) + ^ ( /o" ^^)'^y - ^ = 0- 

with the same nonlocal constraint (15.410 . It is still a difficult system of a first and a 
second-order nonlocal evolution PDEs, 

i.e., actually, is a system of higher-order equations. 

As we know, regardless its solvability, the linearization of this system about constant, 
singular, or other equilibria defines the linear operators that may generate other countable 
families of non-self-similar blow-up patterns with the evolution on centre and/or stable 
eigenspaces, as in Sections 15.111 and I5.12[ 

It is important to mention that, following the lines of this construction, we necessarily 
arrive at the non-stationary system corresponding to (I5.93p . Since odd-order operators are 
dominated here, a proper "entropy" setting for solutions with weak and strong discontinu- 
ities will be necessary to check evolution consistency of the blow-up patterns constructed^. 



For the Euler equations (ll.3ip . see [180j . 

'Recall that there is still no a successful notion of weak solutions of Euler equations in 3D. 
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This is expected to be a hard problem that, with a clear inevitability, accompanies this 
pointwise blow-up analysis of the Euler equations. 

Finally, we believe that general concepts of swirling rotations, axis precessions, and 
vertex motion of periodic or quasi-periodic nature developed in Section [H] can be apphed 
to the Euler equations fll.3ip . but we will not develop these here and return to the Navier- 
Stokes ones (II. ip : see |180j for further refreshing and rather exotic ideas. Recall that 
typical "roUing-up mechanisms" for appeared bubble caps and other swirling features of 
formation of blow-up singularities were observed numerically even in the axisymmetric 
setting; see [9ll 1181] . etc., though these scenarios of blow-up seem remain still under 
scrutiny, while some of them have been ruled out in other works. 

6. On non- radial blow-up patterns on W2: eigenfunctions with swirl 

In this section, we show how to extend some ideas coming from the simpler model 
proposed in Ohkitani [168] with a countable set of blow-up solutions constructed in [70] . 
We return to the NSEs restricted to the subspace W2 = Span{l,2} in Section EH 

6.1. The system on the subspace W2 with swirl and angular dependence. For 

convenience, we replace 

(6.1) W^W, 

so that the system of three equations from Section 15.41 takes the form 

' Wt = A2W - UWr - i VW^ + W\ 

[/ = i /; zW{z, ^, t) dz - i /; V^{z, ^, t) dz, 
[Ut + UUr+'- VU^ -Iv^-^^u + ^V^ + ^U]^ 

. = [r{Vt + UVr + lVV^ + \UV-A2V-^U^ + ^V)\^. 

Indeed, this is a very difficult system. As we have mentioned, its axi-symmetric version 
admits further, more rigorous study, |168l[70] . Here we demonstrate other aspects of such 
solutions. 

Thus, we apply the same scaling as in (15.291) . (I5.30p to the general system (16. 2p to get 



(6.2) 



(6.3) 



Wr = A2W — I yWy — w — aw^ — uWy — ^ vw^ + w 



2 



U 

y 



= [y{vr + (Jv^ + uVy + \^ vVf, + ^uv - A2V + \yvy + \v - ^u^, + 



y ■ y - ■ -z-^ y ■ -Z y- ■ y^ 'Ay 

The linearization w = 1 -|- F as in ( 15.5ip yields the linearized equation (cf. ( ]5.52p ) 

(6.4) n = B*F + F - aY^ -^^y-'y + Dl^'), 

where D is as in (15.521) and B* is the adjoint Hermite operator (I5.55p . We next use the 
second equation for u in (16. 3p to get 

(6.5) u = l^\P^zYdz-\ilv,dz. 
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Assume now that 

(6.6) Y{y,fi,T) and v{y,^,T) are exponentially small for r ^ 1. 

Then by (16.51) m = | + ... up to exponentially small perturbations, so that small solutions 
of (16. 4p are still governed by the linear operator (15.531) 

(6.7) B*Y + Y -aY^ = B*Y -aYf, + ... as t ^ +oo. 

Consider first the full operator with the rotational part with the adjoint one in L^(M^) 

(6.8) Ll = B*-aD^ and U = B + aD^ {B = A2 + yDy + 3l) . 

Similar to B* (similar properties are available for 2mth-order counterparts [51]), it can 
be shown that L* is a bounded operator from iJ^, (M^) to L^*(]R^, has compact resolvent, 
and a point spectrum only. The only, but a key difference, is that, it seems, for cr 7^ 0, 
(16. 8p does not admit a proper self-adjoint setting, so that the whole point spectrum is not 
necessarily real. Writing down in detail the spectral problem in L^*(]R^), 

(6.9) Lir ^ryy + lry + -\yry + r- ^% = 

we see that the essential use of the angular operators makes it rather difficult. Obviously, 
any radial Hermite polynomial remains an eigenfunction for any a G M. Finding a suitable 
point spectrum and non-radial eigenfunctions with essential //-dependence for (16. 9p can be 
a difficult problem. One can see that a standard separation of variables is non-applicable 
for M . 

Thus, a principal question is whether the operator (16. Sp has enough real eigenvalues for 
construction of linearized blow-up patterns. Let us show that in this sense (16. 8p is quite 
suitable. As an illustration of a local approach to (16.91) . we apply the classic techniques 
|212j to trace out branching of eigenfunctions from of B* at cr = 0. Let the kernel of 
B* — \pl has the dimension M > 2. Then looking for the expansion 

(6.10) = Efli Ck{(T)tplk + + - and A(a) = + as + ... , 
we obtain the following problem (A^ are as given in (16.131) ): 

(6.11) (B* - Xpi)^ = E(fc) cfc(o) [srp,k + irpM ■ 

This yields M orthogonality conditions of solvability 

(6.12) <E(fc)C.(0)[sV';,, + (^;,,)'J,V^ft,>=0 for j = l,2, ...,M, 

that define the eigenfunctions of B*, from which branching is available. Then (16.111) 
yields the unique solution (p that shows the branching evolution in (I6.10p (as usual, the 
dimension of the kernel then can play a part). In other words, (16.81) in the inverse integral 
form can be treated as a compact perturbation of a self-adjoint operator, and classic 
perturbation theory applies }114 ] to calculate the deformation of the real point spectrum. 

This shows that there exists branching for small angular speeds \a\ of eigenfunctions 
with swirl from standard Hermite polynomials for a = 0. For compact inverse integral 
operators involved, these are continuous a-curves are indefinitely extensible and can end 
up at other bifurcation point only or can be unbounded; see [23]. In other words, the 
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swirl-dependent operators (16.111) have enough real eigenf unctions for using in necessary 
construction of linearized blow-up patterns. Some of the related questions of linear oper- 
ator theory for (16. lip remain obscure; e.g., eigenfunction closure for such non-symmetric 
cases (but indeed this looks doable though can be rather technical). 

It is curious that even in the case a = 0, we can find essentially non-radial patterns, 
which can give an insight into a swirl structure of blow-up patterns that have a clear 
(/^-dependence. Namely, we now take into account all the non-radial Hermite polynomial 
eigenfunctions (I5.6ip of B* with the spectrum 

(6.13) a(B*) = {A;3 = 1-1/31, |/?| =0,1,2,...}. 

Therefore, there exist solutions of (16.40 with the asymptotic behaviour as r ^ +oo, 

(6.14) w = l- e^P^t{jp{y, fi) + ... for any \l3\ > 1, 

provided that (16.60 holds on compact subsets in y. Obviously, the third equation in (16.31) 
admits exponentially small solutions v for u being an exponentially small perturbation of 
I according to (16. 5p . Then patterns (I6.14p in the inner region make sense. 
The extensions of such patterns into the outer region is similar. We have : 

(6.15) ^/3(l/,/u) ~ / = |2/P'/j(/i) as y^+oo, 

where the first entry y means (1/1,1/2)^ and fj are homogeneous harmonic polynomials 
being the eigenfunctions of the Laplace-Beltrami operator = D"^ on the circle C M^, 

(6.16) AJ, = -ff, on §1. 
Therefore, (I6.14p yields the outer variable 

(6.17) «;~l-|C|l^l/,(/i), where C = yeP^% = ^ < 0. 

Finally, rewriting the first equation in (16. 3p in terms of the new variable C > by 
the change corresponding to the region governed by the new spatial rescaled variable 
C, we recast the original equation by setting (this change is explained in detail in [70] ) 
^(CjT) = w{y,T). This yields the following perturbed Hamilton- Jacobi (Euler) equation: 

= Af,{d) + e-'P^-A^e + ePP^ [\ ( V,) - vw^] , 

(6.18) 

where A^(^) = -{pp + |)C^^c -\{^, " ^ + 

Thus, eventually, on passage to the limit r ^ cxd, we obtain the steady problem 

(6.19) -[pp + lYh^-l^[j^^zh)h^-h + h^ = Q, MC)~l-C^/,(/i), 

For each fixed /i G [0,27r), this is an ODE, which was the object of a detailed study 
in [To]. Therefore, for any /j(/i) > 0, such a profile always exists and is compactly 
supported with non-radial support. Unfortunately, for fj^jj) < 0, the ODE (I6.19P gives 
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a monotone increasing solution which moreover blow-up in finite (. Therefore, such non- 
radial patterns can be constructed in any connected sector 5*^ = {fi : fj{fi) > 0} with zero 
Dirichlet conditions at the boundary rays. There holds: 

1/31-2 

(6.20) ( = x{—t) 2|/3| ^> separable standing-wave blow-up for \(3\ = 2. 

This case includes the radial pattern already studied and others, which are no symmetric. 
For any \(3\ > 2, the patterns exhibit global blow-up as usual in the fixed angular sectors. 

We hope that using non-trivial rotation a ^ will supply us, via the eigenvalue problem 
(16. 9p . some extra new eigenfunctions that allow us to get a non-radial blow-up pattern 
with swirl in the whole space. 

7. On blow-up patterns concentrated about Slezkin-Landau singular 
solutions of a "submerged jet": first example of linearization 

We now discuss some ways to construct necessary blow-up patterns to be applied in 
greater detail in Section [8] by using the spherical geometry. Namely, we introduce impli- 
cations of using the following classic singular solutions of the NSEs: 

7.1. Singular homogeneous equilibria and their properties. In 1934, Slezkin [196] 
(see also comments in [197]) showed that (11.11) admit special stationary solutions with 
the singularity and spatial decay of the velocity field ~ ^ by reducing the problem to 
a linear hypergeometric-type ODE; see more details in Appendix C. In 19440, Landau 
|132j found a family of explicit solutions of that type describing steady fiows induced 
by a point source, which leads to the setting of a submerge jet that is oriented along 
the positive part Oz of the 2;-axis. This gives the following one-parameter family of the 
explicit Slezkin-Landau singular stationary solutions {usl(x),psl(x)}, of (11.11) : 



(7.1) 



2(cz—r)x 
{cr—z)'^r ' 
2{cz—r)y 
{cr—z)^r ' 
2(cr^-2zr+cz^) 

{cr—z)^r ' 
A{cz—r) 
{cr—z)^r ' 



MIX) = 7^ , 



w(x) 
p(x) 



where r"^ = x"^ + y'^ + z"^ (= |xp), and c G M is a constant such that |c| > 1. In view of a 
strong singularity at a; = 0, the existence of such a solution demands an extra force at the 
origin, so actually (17. ip is a fundamental solution of the stationary operators of the NSEs 
satisfying, in the sense of distributions (this computation was already done in Landau's 
original work [Ml P- 300]; see also P pp. 2-9], [nSl P- 182], or [2HI p. 250]), 

(u- V)u-h Vp- Am = 6(c)(5(x)j, V ■ u = 0, where j = (0, 0, 1)^ and 

(7.2) 

Kc) = 3(|^[2 + 6c-3c(c2-l)ln(^)] =i^-|f + ... as c^+oo. 



'^'^In Landau-Lifshitz's monograph [1331 P- 81], this solution dated 1943, with no reference shown. 
Possibly, this reflects the year, when the solution was actually constructed by Landau, but not published. 
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Thus, physically speaking, this means that 

(7.3) the steady S-L solution demands a permanent fluid injection at the origin x = 0. 
Thus, usL satisfies (this will be used in the rescaled blow-up variables y): 

usl(x) ~ ^ as X — s> 0, and, more precisely, 

(7.4) usL(x) = iuo + 0(^) as c +oo, where uo=(^,^,f)^, 

|usl(x)| ^ oo as c ^ 1+ on the semiaxis Oz^ = {x = y = 0, z > 0} . 

To illustrate the last property of unbounded steady profiles, consider the first component 
in (17. ip for a fixed 2; > (for z = 0, the estimate is similar with 4 2 at the end): 

(7.5) m(x)| = , = - '"^y,^^^ -.00 as x^ + y^^O. 

This is rather impressive (and promising for blow-up around): it turns out that arbitrarily 
large vector fields (of special structure) can be locked in a singular steady pattern. 

In the spherical coordinates, the flow (17.11) . being symmetric about the polar (Oz-) 
axis, reads (this is the actual Slezkin [196] and Landau form [132j : see also |133[ p. 82]) 

a\ — l+cos2 e-2ccos6» 



V'"^^ r(c-cos6»)2 ' 

(7.6) <!^(^^) = K^' ^ = 0' 

/ n\ _ 4(ccos6>-l) 
PV^^J — r2(c-cos6»)2- 

Then, in accordance with the divergence (17. 5p . we have 

(7.7) v(r,e)\ = 2 00 as ^^0. 

^ ^ ^ ' ^\c=l r sm(|) 

The singular S-L solutions are locally, but not globally: 

(7.8) USL e Ll^iM.^), but USL ^ L\R^) for all |c| > 1. 

It is also worth noting for further blow-up use of (17. 6p that the total mass flux through any 
closed surface around the origin is equal to zero jl33i p. 83]; cf. the "vanishing oscillatory 
property" (15.460 . It turns out that (17. ip are the only possible stationary homogeneous of 
degree —1, regular except the origin (0, 0, 0) solutions of ( 11. ip |204j : see also earher result 



proved under the z-axis symmetry assumption. One can see that (17.10 are invariant 
under the scaling group in fl2.18|) . Note that the whole set of possible steady singularities 
of the NSEs is still not fully known; even in the class 

(7.9) |u(x)| < where > is a constant. 

However, an essential first step in this direction was done in Miura-Tsai [155j . who proved: 
if u is any very weak steady solution of the NSE in M.^ \ {0}, then 

(7.10) u(x) satisfies (17. 9p for some small C* > =^ u = usl for some |c| > 1, 

where, by (17.50 . (17. 7p . c is assumed to be large enough. In any case, the S-L solutions 
(17. ip should play a crucial role, since these are expected to be isolated from other singular 
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ones (if any); see further details in |204l 1155"] . In addition, there are no smaller singular 
equilibria (see the result and earlier references in [117] ): if u, p G C°°{Br \ {0}), then 

(7.11) |u(x)| = o(|^) as X => is removable and u G there. 

A similar removable singularity theorem for the nonstationary NSEs [125J, where (17.111) 
is assumed to be valid for a weak solution u(x, t) in some set Br x (0,to); includes a 
smallness condition on u in L°°(0,to; L^)- 

Note that in addition to the symmetry of these solutions about the z-axis, fl7.ll) also 
implies that such a flow does not exhibit any swirl, i.e., the angular component w = 
= in fl7.6p . This, in view of global existence results for cylindrically axisymmetric 
(irrotational) flows [211^ I13U] (see |222j for details), indeed, reduces the chances to get a 
reasonable blow-up pattern moving along the corresponding quasi-stationary manifolds (a 
finite-dimensional centre, if any, or a stable one of infinite dimension). Nevertheless, the 
situation with the steady manifold induced by (17. ip is not that hopeless. For instance, 
instead of the "standard" similarity swirl given by the log-law (15.320 . one can use a slower 
rotational mechanism by setting 

(7.12) ^ = + (^f^it), where (-t)/t'(t) ^0 as t ^ 0~ . 

Then, the stationary term cru^ in (15.350 is replaced by the asymptotically vanishing one, 

(7.13) au^^a{-t)K'{t)u^^O (e.g., /^(t) = for t ^ 0", with 5 > O), 

on bounded smooth orbits. Therefore, for a suitable class of solutions, passing to the 
limit along a sequence {rk} +oo by setting t \—>- Tk + t will fix us the previous limit 
irrotational equations admitting the "stationary" singularity (17.10 . Of course, this leads to 
a delicate matching procedure of connecting the slow whirling flow characterized by (17.120 
with Slezkin-Landau's "quasi-steady" solutions (17. ip . which actually should determine the 
function k{t). We return to other extensions and use of (17. ip in Section [HI 

7.2. On some hypothetical extensions. On the other hand, it is principal for blow-up 
patterns to know whether the rescaled NSEs admit other solutions with swirl. In general, 
classification of singular states for the NSEs is a difficult open problem (see [204j ). though 
it is a necessary step for understanding of existence/nonexistence of finite energy blow-up 
patterns for such dynamical systems. 

For instance, it is crucial to check existence of singular solutions with the TW depen- 
dence in the angular direction, with 

(7.14) iph^ip + at. 

This gives the nonlinear eigenvalue problem (we assume that the singularity has not been 
essentially changed, so we may keep the same its ^-interpretation, for simplicity) 

(7.15) au^ + (u- V)u + Vp- Au = 6(a)5(x)j, V ■ u = 0. 

For (Jo = and fe(O) = 6(c), this gives the S-L solutions satisfying (17.20 . Do other 
eigenvalues {ak G M} exist? A partial negative answer is available (see |155] as a guide), 
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but it does not cover the whole range. The corresponding singular states {ugL} then can 
be used for constructing various blow-up patterns. 

In the rescaled NSEs (12.211) . f l2.23p . one needs to know existence of regular states with 
the blow-up angular dependence f l5.30p (for cr = 0, no solutions [162] ) 



Here u{y) may be assumed to be bounded at y = 0, so that it is natural to suppose that, 
along some subsequences or continuous branches (see Section EH for Type II using), 



uniformly on compact subsets in \ {0}. In other words, the singular S-L solutions 
could serve as "envelopes" of the regular ones (this is obscure and questionable). Then, it 
seems, at least a two-parameter family of such regular profiles might be expected. In any 
case, the extra nonlinear eigenvalues {ck} (or other perturbation-like mechanisms) can 
essentially affect the structure of those hypothetical regular swirling states {u^} solving 
(17.161) . We again state that these speculations are made under the clear absence of any 
clue on existence of those singular quasi-steady patterns with various blow-up swirls. Both 
clear mathematically existence or nonexistence conclusions are desperately needed here. 

7.3. On possible blow-up patterns: towards "swirling tornado" about S— L sin- 
gular equilibria. Inner Region I. We now perform first steps of the blow-up strategy 
according to the blow-up scenario developed in Section [L6] now applied to the S-L singular 
stationary solution. So we are going to check whether a "blow-up swirling tornado" (not 
a twistor, since we do not always apply directly the logTW mechanism, though it is not 
excluded) can appear by approaching the L°°-singularity at x = by drifting as t — T~ 
by "screwing in" about this singular equilibrium structure. There is a clear suspicion that 
this can hardly happen without, according to (17. 3p . a permanent fluid injection, which 
is not available for bounded L^-solutions. Note that the hypothetical blow-up tornado is 
going to occur very fast, during a miserable time scale, when the total injection is neg- 
ligible. Anyway, this physical negative motivation should find some clear mathematical 
issues in support or not. 

Recall that the behaviour in Inner Region I is assumed to be characterized by the 
linearized problem (I5.77p . Thus, following the same standard principles of matching, 
linearization of the general equation (I2.2ip about usl via (I5.75P yields, in Inner Region 
I, the equation (15.760 with quite a tricky linear pseudo-differential (integro-differential) 



operator obtained from (12.230 . which in C(j"(M^ \ {0}) reads 

H'(usl)Y = (B* - 1 /) Y - PCY, where CY = div (usl ® Y + Y ® usl) 
(7.18) = (^SL • V)Y + (Y ■ V)usL, so H'(usl)Y = (B* - i /) Y - CY 



(7.16) 



au^ + H(u) =0 in M^ u G ivf„^(M^) (u ^ L'^{R^) in general.) 



(7.17) 



u(0) = C and u usl(c) as |C| +oo (|c| > 1) 
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Here, Ji includes all the local differential terms, while J 2 the integral one. Note that, in 
general, according to fl7.16p . the above linearized operator must include the operator 



(7.19) ... — crD^Y + ... , so fl7.18p contains two free parameters: c and a, 



with a nonlinear eigenvalue a coming from the stationary problem already containing the 
parameter c > 1. In a most general setting, the linearization is performed relative to 
any unknown S-L-type singular steady states delivered by fl7.16p . For simphcity, we will 
continue to deal with the standard S-L solutions, actually meaning that, in addition to 
other branching ideas, one needs also (cf. Section [6TT1) : 



(7.20) to develop branching of eigenfunctions from the logTW speed a = 0, 



though possibly this is not an issue since demanding too much from the operator (17.180 
(A = must be in the spectrum). On the other hand, even if (I7.20p makes sense, it 
is then also plausible that such a branching is available (or not) for |cr| small only, and 
proper eigenfunctions may occur via "saddle-node" bifurcations for large \a\, which are 
very difficult to detect. In any case, further detailing of these arguments seems excessive, 
since, even without (17.200 and/or others, the study gets very complicated with already 
too many open problems to appear. 

In all the cases, we observe in (I7.18P singular unbounded coefficients as ?/ ^ 00 as well 
as y ^ 0. A proper treatment of ?/ = 00 is settled via the space L^, (R^), as usual. The 
singular point y = is much worse. 

Local operator: discrete spectrum by improved Hardy-Leray inequality . 
Namely, using (17.40 yields that the singularity of the potentials in the local terms, 

(7.21) CY ^ (usL ■ V)Y + (Y ■ V)usl ~ (^ ■ V)Y + (Y ■ V)^, 

is well covered by Hardy's classic inequality at least for all large c > 0, since both terms 
"act" like the inverse square potential. Thus, according to (ll.43p . we need 

(7.22) J,Y^AY+^ as ,-.0 ^ - ^ < ^U^3 = i 

Note that, for axi-symmetric solenoidal flows, the Hardy optimal constant ch = -^^ — 
becomes better, and, as shown in Costin-Maz'ya [44], can be replaced by 



I'-^JJ f-H.axi-sol. — 4 Ar2+2Ar-4hv=:^ ~ 68 ~ U.OUm... ^ 



4 Af2+2Ar-4 I Ar=3 

In other words, the operator AY + in the axi-symmetric solenoidal field admits a 

proper setting in for x ~ 0. Note that the inequality 

(7.24) i /m3 ^ < /m3 IVYP (AT = 3) 

already appeared in the same Leray's pioneering paper in 1934 [137] . Thus, ^ in (I7.24p can 
be replaced by || that improves the admissible range of c's according to (17.220 . Hence, 
we obtain the first important and even principal conclusion: 

(7.25) for any c ^ 1, the local linearized operator Ji is well-posed. 
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Therefore, this local part of the linearized operator has a compact resolvent, a discrete 
spectrum in the same weighted L^-space, etc., though several questions including com- 
pleteness and closure of the eigenfunction set and others remain open and difficult in this 
non self-adjoint case. Moreover, it follows from f l7.18p and f l7.2ip that for c ^ 1, with a 
proper functional setting, Ji has eigenvalues that are close to those in (12.251) for B* — | /. 
In addition, by classic perturbation theory (see Kato |114j ). it can then be shown that 

eigenfunctions \'v%} of Ji can be obtained from generalized 

(7.26) I fiJ ^ to 

Hermite polynomials {v|,} via branching at c = +oo. 

This justifies existence of infinitely many real eigenvalues of Ji, at least, for c ^ 1 (global 
extension of bifurcation c-branches to finite values of c > 1 is a difficult open problem). 

Since in (I2.2ip P is a projector in {L'^Y, ||P|| = 1, it seems reasonable to expect that 
some part of discrete spectrum may be governed by the local differential operator Ji. 
Consider the full linearized equation 

(7.27) = (B* - i /) Y - PCY + D(Y) = H'(usl)Y + D(Y), 

where D is a quadratic perturbation. Assuming that there exists a proper eigenvalue 
A/3 e crp(H'(usL)) 7^ 0, with ReA/3 < 0, we expect the behaviour Y(r) ~ ce^^'^v^ for 
r ^ 1, if Re A/3 < 0, or (as in Section [5.121 Y(r) ~ C/3(r)v|, close to a centre subspace, if 
Re A/3 = 0. [Both to be matched with more regular flows around.] Substituting this into 
(I7.27P yields the eigenvalue problem, where the nonlocal term implies a tricky setting, 

(7.28) (B* - I /)v; - PCv; = A/3v; (v^ e H^p^M.^), at least.) 

Again, it should be noted that the angular operator (17.190 must enter this eigenvalue 
problem. We claim that this extra technicality can be tackled by separation in the angular 
variable, though this leads to more difficult calculus and a proper adaptation of further 
conclusions are not easy but doable. However, this shows a natural opportunity to be 
retained within local spectral theory for Ji. Since, as we have seen, PC — as c ^ +oo, 
this demands solving the following problem for generalized Hermite polynomials {v^ G 
<|.*3, A; = 1^1} as in Section ED 

(7.29) 



c ^ 1: for which v^, the linearized term Cv^ is or "almost" solenoidal? 



Obviously, (17.280 shows that if Cv^ is fully divergence-free, then P = J on it, so the local 
operator Ji takes the full power, and this settles the spectral problem for c ^ 1. As we 
expect, this will allow us to start a series-like expansion in terms of the small parameter 
^ of the necessary div-free eigenfunction of the full problem (I7.28p . 

We did not check if such eigenfunctions are actually existent, but, in a whole, the 
problem (I7.29P looks doable. Indeed, we do not expect to find such polynomials 
satisfying (17.290 for sufficiently small |/3| = 0, 1 or 2, etc. The first such / = for which 
(17.290 holds with necessary accuracy, would show a complicated geometry (with remnants 
of swirl, axis precession, and so on), which would eventually affect the global structure 
of a possible blow-up pattern. Once we have found such solenoidal generalized Hermite 
polynomials satisfying (17.290 . this starts the not less simpler procedure of constructing 
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those "twisters" by matching with more regular flows close to the singular point y = and 
outside the pattern blow-up region (the latter is supposed to be simpler). Of course, it is 
too early discussing any seriously this nonlinear matching matter (but we will, noting also 
that we have already mentioned some formal principles of such procedures that are well 
established in blow-up R-D theory). But indeed the linearized spectral analysis about 
the rescaled S-L solution is assumed to show some basics of formation of such mysterious 
blow-up patters including their strong angular motion (the eigenfunctions are supposed to 
be essentially non-radial), "precession" features of their axis of rotation (to be discussed 
in greater detail), and even their possible merging into a closed path of several individual 
twistors (hopefully, not infinite). 

Nonlocal operator: an obscure functional setting . Anyway, despite specu- 
lations concerning the more or less standard local operator Ji (with basics related to 
Hardy-Leray inequalities), one cannot avoid analyzing the nonlocal integral operator J2. 
The spectral questions for the nonlocal operator in (17.181) become more involved. Consider 
J2 in C(f (M'^ \ {0}), where, on integration by parts, 

J2Y ^ C3 / ^ E {^.Yi + uI^^y:^ 

(7 30) ^'''^ 

We obtain singular integral operators with the kernels, formally exhibiting the behaviour 

(7.31) Kiy,z) = K^{y,z) + K,{y,z), where lir^l ^ 1 ^ ^, jir^l ~ ^ ^ ^. 

First view of Ki : too singular, a condition at needed. The kernel Ki is standard polar 
relative to the first multiplier having ~ j^rip) with fi = 2 < N = 3. However, the second 
multiplier, at first sight, exhibits a non-integrable singularity j^p- at 2; = in demands 

special conditions at the origin e.g., of the type (otherwise, the integrals in f l7.30p ought 
to be understood in a generalized v. p. sense, see below) 

(7.32) \Y{z)\ = 0{p{\z\)) ^0 as z 0, where /g^dr<oo. 

On the one hand, it seems clear that any such condition (17.321) on the behaviour at the 
origin can reduce a hope to find a proper eigenfunction of fl7.18p . On the other hand, it is 
well known that eigenvalues of pseudo-differential operators can be of infinite multiplicity 
(e.g., for those with constant symbol), i.e., the admitted behaviour at singularity points 
can be extremely various (unlike the differential operators, for which such a behaviour 
is restricted by their orders and asymptotics of the coefficients). However, even for the 
local differential operators, imposing such conditions is not that hopeless, as the following 
example shows: 

Example: differential operator with a locally non-integrable potential. Consider 
for simplicity a scalar differential operator associated with Ji in (17.181) . with an extra 
inverse cubic potential (^ LfQ^(R'^)) as in (17.311) : 

(7.33) E = B*-i/-C-^/ (A^O). 
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Obviously, looking for possible eigenf unctions, the behaviour of solutions close to the 
singular origin is defined by the principal part, so that, as z — > 0, 

(7.34) Eip* = Xi)* =^ Ai;* - j^ip* + ... = =^ ip*{z) uj{-^) \y\-^e^\'\'' , 

where u on the unit sphere §^ in M'^ makes the angular separation. Substituting the last 
WKBJ-type expansion term into the operator (ip*)" + ]f| (V^*)' ~ j^s'^* — ^ yields 

(7.35) a = -h S=z, and = AA =^ a± = ±2y/A for A>0 

{A > corresponds to "stable" potential, i.e., monotone principal operator). Thus, close 
to 0, there exist two asymptotic bundles of solutions with essentially different behaviours: 

(7.36) (2:) ~ |z|"t e~^^/^ and ip^z) ^ \z\~^ e^"^^ ^ 00. 
Therefore, posing the condition (cf. fl7.32p ) 

(7.37) Y{0) = 

does not spoil at all the eigenvalue problem for E and just eliminates those singular be- 
haviours at (half of all of them) making the spectrum discrete. In linear operator theory, 
the same is usually and naturally done by assuming Y G L'^{Bi), which nevertheless has 
come true after certain asymptotic analysis presented above. For the principal radial 
ordinary differential operator, this means that the singular point |z| = is in the limit 
point case (an index characterization is also available for the elliptic setting). 

On the contrary (and this case seems more correctly describes the nature of such an 
"unstable" potential in (17.311) . though not that obviously), if A < 0, then (17.351) yields 
a± = ±2i a/|v4|, so that both bundles are singular and oscillatory: as z — * 0, 

(7.38) C(^) ~ cos (2^^), ^;(;,) ^ |^|-! sin (2^^), r± e Ll^; 

the limit circle case for the operator A — I. Then the discrete spectrum of E is also 
achieved by posing special conditions at the singularity, |160] . and then any eigenfunctions 
exhibit the singular behaviour (I7.38p . It is another hard problem to check whether (17.381) 
allows to match such an oscillatory bundle with smoother L°°-flow near the origin z = 
for r ^ 1; see below. Therefore, in the case of total oscillatory bundle (I7.38p . the 
condition such as (17.371) makes no sense and, evidently, destroys the eigenvalue problem 
giving cr(E) = 0. 

Note finally, that the oscillatory- type behaviour similar to (17.381) . in general, is not an 
absolute obstacle for getting from such "eigenfunctions" a proper blow-up pattern. We 
discuss this for some R-D equations in Section [?.10[ Then, for the behaviour like (I7.38p . 
the integral operator as in (17.301) should be understood in a (canonical) regularized sense, 
see below. 

Singular kernel K2 is of Calderon-Zygmund type. The integral operator with the kernel 
K2 in (I7.3ip is singular. However, in view of the divergence form of the operator C 
in (I7.18p . there is a hope (not yet fully justified) that this kernel falls into the scope 
of Calderon-Zygmund's classic result (1952) [28j| saying that such an operator can be 
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bounded in Lp(R'^) (in the local sense, meaning that, as usual, we cut-off the infinity 
by the appropriate weight p) for some p > 1; see |2U6l Ch. 5, § 5] for modern overview 
and references. Most of the text-books on pseudo-differential operator theory quote such 
fundamental results as being its origin; see |151l p. 278]. Note that the conditions of the 
fundamental Calderon-Zygmund result in the standard form |206l p. 16] do not directly 
cover the singularities of K2, so that an extra hard work is essential. In addition, one 
surely needs extensions of such boundedness results to L^-spaces of vector- valued functions 
Y (such a study was already initiated by Calderon himself in 1962); these questions being 
well understood; see references and a survey on further operator- valued issues in |105] . 

Back to Ki: the operator at in a c.r. sense. Indeed, such a possibility to "overrun" (if 
possible) the very restrictive condition (17.321) is to take into account the specific divergence 
part of the jjp-terms in (I7.30p . To explain this, let us fix i = j = 1 in the first term, 

where for z ^ and y ^ 0, for G C^, we have that the singular part in the regularized 
value sense vanishes: 

(7-39) ~C73^rH0)4KL),^,^d. = 0, 

since the second-order zi-derivative, according to (17.11) . is an odd function in zi. Namely, 
the corresponding part of the integral operator in (I7.30p acts like the standard distribution 
with the regularization 

^1 

(7.40) (p^,y,) = J-^^(£H^(^^WMds for ipeC^{R), 

which is canonical, c.r. (i.e., this regularization keeps the linear properties of the function- 
al, as well as the differentiation). A similar c.r.-property is observed for the whole Y^- 
term appeared in (17.301) that consists of three members: '^■[{uli^zizi + iusL)ziz2 + ("^31)2123] • 
Note that all of them are odd relative Zi, so the indefinite integral is even. It seems that 
direct L^'-theory is not applicable to such c.r.-integral operators that thus deserve further 
study. Elsewhere, in other, less singular terms, the standard v. p. -sense c.r. 

^ zi' ' I J — 00s JO s 

in the Cauchy sense can occurs, leading to more standard integral operators with typical 
singularities of the Hilbert transform in zi (a Calderon-Zygmund operator in M), etc. 

Continuing the fruitful idea of eigenfunction branching at c = +00 as in the local case 
above, we obtain an open problem for the whole pseudo-differential operator (I7.18P : 



(7.41) In which operator topology, does H'(usl) ^ B* — ^ / as c ^ +cxd? 



Or the limit always contains a certain "singular" part? The answer is principal. 

Note that any extra condition such as (I7.32p can be inconsistent with admitted "weakly 
singular" behaviour corresponding to the local operator Ji composed from B* and C 
in (I7.2ip provided that the latter ones are "dominant". Then this case becomes rather 
standard and, as in blow-up R-D theory (see Section [L6|) . is governed by improved Hardy's 
inequalities discussed above. 
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However, if the nonlocal operator J2 is "dominant" close to the origiro, which seems 
unavoidable and represents the main difficulty of blow-up in the NSEs (cf. the uncertain 
polynomial micro-structure of multiple zeros detected in Section [3] by the same reason), 
then the conditions such as fl7.32p can be natural, giving rise to the necessity of differ- 
ent spectral theory, which is entirely unknown and represents an open problem. We do 
not know whether the linearized pseudo-differential operator (17.181) with such a strongly 
singular kernel (and hence with a nonregular full symbol; see an overview in |206l p. 45]) 
has a nontrivial discrete spectrum in a proper functional setting involving conditions at 
the singularity and Lorentz-Marcinkiewicz or Zygmund-type spaces (see |206[ p. 37] with 
applications to nonlinear PDEs in this Taylor's volume). 

Note that, plausibly, by some mysterious reason, such a pseudo-differential operator 
may admit just a single proper eigenfunction (for further use in blow-up matching, with a 
positive or negative conclusion, it does not matter), which cannot be detected in principle 
by any general advanced spectral theory. This is a difficult problem, where numerics for 
such involved nonlocal operators, can be key, though a definite justified answer: "yes" 
or "no" to existence of proper eigenfunctions and hence centre-stable eigenspaces, can be 
very questionable. 

Overall, spectral properties of the pseudo-differential operator fl7.30p are of great de- 
mand (its symbol is not still well understood) and can be key for existence/nonexistence 
of proper blow-up patterns on centre-stable manifolds created by the S-L exact steady 
singular solutions. Actually, even the negative result on nonexistence would play a role 
for concentrating on other more involved scenarios of blow-up to be focused on. Of course, 
as usual, since the S-L solutions (17. 6p are axi-symmetric with no swirl, our analysis is 
assumed to include blow-up swirling mechanism as in (I5.32p . so that the resulting op- 
erators contain typical angular terms like —cfD^ as in (I7.19p . with, possibly, sufficiently 
large angular speeds a (this assumes complicated separation angular techniques). Thus, 
Questions (i)-(vi) from Section 15.121 do deserve further study in this case, and, with a 
certain luck, will give a ffist insight into blow-up singularities for the NSEs. 

Let us note that, after matching of the linearized patterns on the stable (centre) man- 
ifold for Ji 2 with more regular flow close to the origin (this is even more difficult open 
problem, which makes no sense if spectral theory is still unavailable; Sections 11.61 and 
future rr6] may be consulted for an idea), a typical shape of the resulting patterns (ac- 
cording to fl7.32p or other restrictions) will have a form of swirling "tornado" about a S-L 
singular steady solution, which is self-focused onto the origin x = as t ^ T~ = 0~. 

Final conclusion: spectral results are expected by branching at c = +00. Despite various 
difficulties and open problems already detected, our final conclusion is not fully negative. 
Though the extra condition such as (I7.32p looks rather restrictive and even frightening, 
the branching idea at c = +00 correlated with this well, since: 



(7.42) there exist infinitely many Hermite polynomials satisfying f l7.32p . 



^OActually, as seen from P in (|2.2ip , both local and nonlocal operator parts are of a "similar power" . 
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For instance, fl2.32p implies that those are v^^^, v^g; ^^3, V24, V25, V2g, V27, V2g, etc. 
However, obviously conditions such as (17.321) can mean that u.{y, r) is always singular at 
y = 0, so f l7.29p is indeed more preferable. 

Thus, within the positive reflection of the hard demand f l7.4ip and in view of the desired 
feature (17.420 . we expect that 

(7.43) for 1, H'(ugL) may have several real eigenvalues close to \p = — 

and the total number of those gets infinite as c — +00. Indeed, in this limit, we observe 
a "convergence" to the Hermite operator B* — | /. In other words, the above analysis 
makes it possible to start a real procedure of checking whether at least one from this 
countable set of linearized structures admits a proper matching to get a finite energy 
blow-up pattern for the NSEs (or all of them are hopeless, which nevertheless would not 
prove nonexistence of blow-up since there are other blow-up scenarios). In other words, 

(7.44) S-L solutions with c 1 are not still forbidden for blow-up evolution around, 

while extensions of the branches to finite c > 1 is even more promising, but indeed 
extremely difficult. Then the still mysterious matching procedure gets principal, which 
itself can cross out all the previous "spectral" speculations and illusive achievements. 

7.4. First application of spectral theory: towards periodic blow-up patterns. 

As a first elementary application of the above spectral discussion, returning to sim- 
pler periodic blow-up orbits inducing (I2.46p . we pose a straightforward question on the 
Andronov-Hopf (A-H) classic scenario of bifurcation of periodic orbits from the singular 
equilibrium ugL- Namely, we state: 

to check if an A-H bifurcation can occur at some cr = (Tah; c = Cah 
^ ■ ^ for the operator fl7A8|) . dZHD, i.e., 3 iuj e (T(H'(ugL) + (TarD^), ^ 0. 

The asymptotic behaviour of the corresponding eigenfunction v^(?/) as ?/ — is also of 
crucial importance to get a periodic pattern by matching with the regular bounded fiow 
for y 0. In case of both positive answers, this would lead to a periodic twistor blow-up 
pattern for a ^ a ah that gives rise to the cu-limit set (12.461) . Respectively, existence of 
multiple eigenvalues iui,..., iun may lead to quasi-periodic orbits obeying (I2.47p . These 
are hard open hypothetical questions. 

Let us finish this brief discussion with a negative result concerning bifurcation from 0: 



Proposition 7.1. For the rescaled equation {\2.21\) with the swirl operator ( 7.19), an A-H 
bifurcation from is impossible. 

Proof. Consider the eigenvalue problem for the linearized operator about 0: 

(7.46) (B* - I /)v* - aD^w* = Av* in L^.. 
Using the symmetry (I2.24p and divergence of on = (0, 27r), we get 

(7.47) ^||v*||r2 =-||Vv*||r2 - ^ \W* \\ r2 ^ Re A < for auy (7 G M. □ 
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7.5. Inner Region II: matching with smoother solenoidal flow near the origin. 

Assume that the previously posed spectral problem has been solved successfully, so we 
have found the actual Inner Region I, where f l5.77p holds. For simplicity, we assume a 
stable subspace behaviour: 

(7.48) Inner Region I : u{y, r) ~ Usl(z/) + e^^^v* (y) for r > 1, 

on a certain "maximal" set y G T(r), where such an expansion is applicable. Roughly 
speaking, it is given by 

(7.49) T(r) ~ {r(r) < \y\ < R{t)}, 

where R{t) characterizes the outer matching with smooth and almost regular flow, while 
much smaller r(r) — is the crucial sphere, on which the inner matching with the flow 
bounded at the origin is supposed to occur. Thus, in particular, 

(7.50) set T(r) essentially depends on the eigenfunction v*(?/) behaviour as ?/ — » 0. 

Indeed, as a proper illustration, T(r) is different for the first monotone patterns as in 
(17.361) and for the singular oscillatory ones in (17.381) . For other types of possible eigen- 
functions, the actual matching as in (I5.86p - (l5.88p defines the inner boundary of T(r). 

Concerning the outer boundary of the radius R{t), which is not that essential, it is 
estimated as follows: since v*(|/) is assumed to be close to a Hermite polynomial for 
y ^ I oi order k large and hence A^(c) = ^j^, with |7| = k, this radios is characterized 
by matching with zero: for r ^ 1, 

(7.51) usL(y)+e^--v;(y)~0 =^ ^-e-^-|y|'=~0 =^ i?(r)~e5, 

as should be via the basic kinetic energy estimate (12.21) in the rescaled variables (I2.18p . 
However, as we will see, this outer region cannot be important for the inner matching. 

In Inner Region II, we have to return to the full original rescaled equation (I2.2ip . where 
for convenience we set P = / — Q, 

(7.52) = (B* - I /)u - (u • V)u + Q(u ■ V)u, 

where the convection term (u • V)u is negligible in comparison with the linear B*-term: 
this is the actual definition of Region II. The non-local term keeps be influential and even 
leading therein (we have seen a similar nonlocal phenomenon in the study of polynomial 
micro-structure of multiple zeros in Section [HI where the integral term was shown to 
essentially deform the Hermitian solenoidal fields). 

Thus, in Region II, the following asymptotic equation occurs: 

(7.53) = (B* - |/)u + f*(u) + ..., where f*(u) = Q(u ■ V)u as r ^ +oo. 

It follows (or this can be a key assumption to be checked) that, in view of the expansion 
fl7.48p . the forcing term in (I7.53P can be estimated as follows: 

(7.54) f*(u) ~ f*(usL) ~ Qtm(usl ■ V)usl, 

where in the last term we assume integration in fl2.23p over T(r) only, which represents 
the leading-order expansion term of the nonlocal term. If the latter is not true, one needs 
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to include also the nonlocal portion that depends on the still unknown solution expansion 
in Region I, that makes the analysis more complicated (but formally doable). The Outer 
Region, where the solution is much smaller, is assumed to produce no essential influence 
on the integral. 

Thus, the possibility of matching of Regions I and II depends on existence of a smooth 
bounded solution of the following limit problem: 

(7.55) = (B*-l/)u + f*(^)(r), where f*(^)(r) = Qxm(usl) for r > 1. 
It is convenient to use the generalized Hermite polynomials to describe the solution: 

(7.56) u(r) = Ec/3V^> where = (A/j - |)c/3 + (f*(^)(r), v^). 

The first condition (by no means, a necessary and/or sufficient)) of proper matching 
reads then as follows: 



(7.57) the auxiliary problem f l7.56p has an L°°-solution for all r ^ 1 



(otherwise, uq ^ L°°). As a second one, for a rough checking of matching, consider the 
equation for the leading Fourier coefficient for /? = 0, with Aq = 0: 

(7.58) Co = -|co+(f^(,)(r),vo), 

where vq is a constant vector (a polynomial of degree zero). Integrating this and assuming 
a slow growth divergence, one can expect a "quasi-stationary" behaviour given by fl7.58p : 

(7.59) Co~2(f;(,)(r),vo) for r > 1. 

Assuming that r(r) (otherwise, no blow-up as r ^ +oo), matching the stable 
manifold behaviour about the S-L profile msl assumes, at least, that, for r ^ 1, 

(7-60) M(0,r) ~Co(r) ~msl(x)||^I^^(^) =^ |co(r)| ~ 

Involving other expansion coefficients {c^(r)} will lead to a similar but more complicated 
relation with r(r). Altogether, (17.591) and (17.601) define a complicated nonlinear integral 
equation for the expansion coefficient Cq{t) and the matching radius r(r). On integration 
in the non-local term given via fl2.23p . i.e., 

(7.61) fr^l^) ~ Ir(r) Mk,.,^^,.,' where ~ ^, 

it follows that the RHS in (17.591) behaves as ~ ^ as r — 0, which is satisfactory for 
a possible matching purpose (at least, an obvious contradiction for matching is not an 
immediate option). Surely, these are just rough matching estimates, and further more 
difficult matching study along the lines of (I7.56p . (17.590 is necessary. 

Thus, due to (17.501) . the "quality" of the eigenfunction v*(?/) as ?/ — > is essentially 
involved into (17.571) via the set of integration T(r) for r ^ 1. This somehow involves 
verifying a certain "integrability" (sufficient integral convergence) over T(r), where the 
asymptotics of the eigenfunction v* (y) as ?/ — is key. In other words, this eigenfunction 
must belong to a "proper functional class" to make the above computations meaningful. In 
general, similar to the construction in Section [L6| the actual blow-up structure in Region 
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II will indeed depend on the behaviour of v*(?/) as ?/ ^ (unknown), and hopefully will 
look like a smoother "quasi-stationary" (driven by A only) evolution. On the other hand, 
even for very bad, singular and/or highly oscillatory eigenfunctions v*, there is still some 
plausible hope that the integrals in (17.561) still properly converge in such a manner that 
(17.571) holds (this can be checked), and the matching with a bounded flow for y ^ can 
be purely "accidental" . This means that (I7.57P can be valid for some special and possibly 
very singular eigenfunctions v*, so that the functional setting for the eigenvalue problem 
for the occurred pseudo-differential operator can be treated in a wide sense. 

Further matching conditions, involving also the already obtained Region I expansions 
(and the Outer Region if necessary), may be important, which eventually give the matched 
solutions in terms of a converging bounded functional series. Finally, in view of the 
precaution (17.31) . another fruitful idea is to diminish the fluid injection by assuming that 

(7.62) c = c(r) +00 as r — +oo, 
which is an individual subject of Section [7771 

7.6. A scaling view to formation of Type II blow-up solutions: heteroclinic 
orbits are necessary. We then need to consider the full rescaled equation taking into 
account the angular operator (17.191) : 

(7.63) u^ + au^ + P(u ■ V)u = Au - i y • Vu - ^ u. 

Assuming that the blow-up is faster than the self-similar one, i.e., is of Type II and 

(7.64) supy \u{y,Tk)\ = Ck oo as {rk} oo, 
we apply the Cfc-scaling technique as in Section H] by setting 

(7.65) u = CfcW, y = yk + a-kZ, t = + a|s, /i al^, where = ^ — > 0. 
The resulting equation for w^ takes the following perturbed form: 

(7.66) w, + (Tw^ + P(w ■ V)w = Aw - ^ ■ Vw + iw). 

k 

In view of uniform boundedness and further regularity of the sequence {w^.} on compact 
subsets in x M, we can pass to the limit A; — > oo in the weak sense in (I7.66P to conclude 
that Wfc(s) must approach a regular solution W satisfying the NSEs 

(7.67) W, + (xW^ + P(W ■ V) W = AW in x M, ||W(0)||oo = 1- 

Note that W(2;, s) ^ is an ancient solution, which is defined for all s < 0. At the same 
time, by construction, it is also a future solution, which must defined for all s > 0. Indeed, 
one can see that if W(s) blows up at some finite s = S~ > 0, this would contradict the 
Type II solution w(?/,r) is globally defined for all r > 0. Thus, by scahng of the Type II 
blow-up orbit (17.651) . we arrive at the problem (I7.67p . which defines: 

(7.68) { heteroclinic solution W(2;, s) ^ 0} = {ancient for s < 0} U {future for s > 0}. 

Note that constant stationary solutions Wq of (I7.67P are the simplest possibilities, 
which being rather suspicious (too elementary) have not still been ruled out. We do 
not know whether or not non-constant or non-steady global bounded regular solutions of 
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(17.671) exist for some cr 7^ 0. Note that in the stabihzation problem to a constant solution 
Wo, which appear by linearization: 



(7.69) W(s)=Wo + v(s) =^ V, + av^ + P(Wo ■ V)v = Av, 

the linear operator does not have a clear discrete spectrum, unlike the one appeared in 
Region I. Therefore, the actual exponential stabilization rate in Region II is inherited from 
the spectral problem in Region I, assumed to be properly solved beforehand. In any case, 
this analysis shows a quasi-stationary nature of formation of Type II blow-up patterns on 
shrinking spatial subsets around blow-up points. Of course, the actual details of such a 
behaviour can be much more complicated involving various extra singularity mechanisms. 

Regardless a clear lack of rigorous arguments here on existence of Type II blow-up 
patters, a negative conclusion occurs (surely, not surprisingly): 



(7.70) 



if the hypothetical heteroclinic patterns W{z, s) from (I7.67P do not match 
for 1^1 ^ 1 typical spatial structures in (17.11) or those in (17.161) . 
tornado-type blow-up around the S-L singularities is not possible. 



7.7. Blow-up drift on the c-manifold of S— L solutions. As we know from R D 
theory, the above blow-up scenario does not exhaust all the types of possible singular 
patterns. As a new evolution possibility, we again consider the S-L solutions, but currently 
we assume a centre-subspace-like evolution on their whole manifold. 

We begin with the first linearized procedure to calculate a first approximation of such a 
behaviour. Namely, following [72], we assume that slow blow-up evolution occurs "along" 
an unknown functional dependence c = c(r) in (I7.18p . i.e., via (17.621) . 

(7.71) u(r) = UsL(c(r)) + Y(r) G Lj,(M^), where c(r) ^ +00 as r ^ +00. 

Then, instead of (I5.76p . using the second asymptotics in (I7.4p . we substitute (17.710 into 
the NSEs with the corresponding extra force on the right-hand side as in (17. 2p of the 
form = -^S(y)j + ... . Omitting higher-order terms, we obtain the following perturbed 
inhomogeneous PDE: 

^2) - ^ uo = H'(usL(c(r))Y + D(Y) + ... 

= (B* - i /) Y - i [(uo • V)Y + (Y ■ V)uo] - i C3 /^a (Ol^o, Y) + ... . 

In the last term, we mean the nonlocal operator in (I7.18P defines (still formally) at the 
constant vector field Uq, with a special truncation at the origin. 

We next looking for a solution governed by the leading perturbation in (17.720 . 

(7.73) Y(r) = if{T)^*{y) + w(r), where (/?(r) = ^0 as r ^ +00, 

and we assume that w(r)±\I^* in the dual metric, i.e., with a proper definition of the 
adjoint element (a linear functional), 

(7.74) (w(r),vl>)^0 (vf G ^^(R^)). 
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Substituting (17731) into (17721) yields 

(7.75) 0** + ^uo = <^(B* - i/)vl/* - f [.](uo, VI/*) + ... , 

where [■](uo,^*) denotes the hnear operators in (I7.72p at uq (with possibly a necessary 
truncation at 0). 

Countable set of exponential patterns. Formally, it then follows that (I7.75P admits a 
countable set of linear solenoidal blow-up patterns (here = — |) 

(7.76) Yk{y, r) = ak{T)\l{y) + ... , where - (uq, v^) = (Afc - \)ak. 
In particular, we have the following special family of functions {cA;(r)}: 

(7.77) Cfc(r) = e(5-^fe)- =^ ^^(r) = (i - A^) r e^^""^)" for r > 1. 
Power decay pattern. Assuming now that in (17.731) . 

(7.78) |0(r)| < |(^(r)| for r > 1, 

we mean that c(r) is a slow growing function not of an exponential form. Then the only 
possible way to balance the terms in (17.751) is the vector ^* to satisfy: 

(7.79) (B* - i/)^* = Uq. 

Since the operator is invertible, there exists a unique solution \E'* constructed by an 
eigenfunction expansion via solenoidal Hermite polynomials, so is also solenoidal. It 
follows from (17.41) that this linear procedure gives a similar singularity at the origin (cf. 
below with a nonlinear one), 

(7.80) l^*(l/)l~^ as y^O. 

The further balance in (17.751) is then obtained, as usual, by multiplication by the adjoint 
orthonormal element ^ that yields an asymptotic ODE for ^{r) for r ^ 1: 

(7.81) = 7o| + ..., where 70 = (^(uo, ^*), ^). 

We must admit that actually, 70 is not supposed to be a constant, since in this dual 
product, a further matching of the pattern with a bounded flow for y k, should be 
assumed. This matching (a cut-off for \y\ <^ 1) procedure makes the integral in (17.81!) 
finite for all r ^ 1, but then we conclude that 

(7.82) 7o = 7o(r)^oo as t ^ +00 (e.g., 70 ~ In r), 

though the divergence turns out to be slower than that for v5(t), so this does not essentially 
affect the leading term in (17.811) . Thus, we arrive at the equation for (^{t): 

(7.83) (^)' = ^o^ + ... =^ c(r)~-^r + ... as r +cx) 

(recall that we do not exclude the case, e.g., c(r) ~ rlnr for r ^ 1). 

Hence, up to lower-order multipliers, we get the following expansion of such solenoidal 
patterns, which is convenient to write in terms of a series of the form (again, lower-order, 
logarithmic-like factors are not taken into account): 

(7.84) u(y,r) ~ iuo + :;3_u^ + J_U2 + ..., 
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where for convenience we put ui ~ To get the coefficients {ufc} of this nonhnear 
expansion, one needs to use the full rescaled equation (I2.2ip with the operator (12.231) . It 
then follows that the equation for the third coefficient U2 takes the form 

(7.85) -ui = (B* - i J)u2 - (uo, V)ui - (ui, V)uo - J^,{-){uo,ni), 

where in the last term we again assume a certain evolution cut-off procedure near the 
origin. Then we obtain the singularity for U2(y): 

(7.86) _(uo,V)ui~^ =^ My)-1^ as y^O, 

i.e, we have found that the third term is more singular at the origin than two previous 
ones. Therefore, similar to the procedure described in our ffist Type II blow-up structure 
in (I1.45P for the Frank-Kamenetskii equation, we observe that on small compact subsets 
in y, a model one-sided (say from above) behaviour of the rescaled vector field can be 
roughly estimated as 

(7.87) |u(y, r) I ~ U{y, r)~i^ + J,^ + ^MM + .... 

Taking into account the ffist and the third terms, the absolute positive maximum of the 
scalar function U {y, r) is attained at 

2 t2 

(7.88) \y\ ~ e"'^ — > 0, and hence snpyU{y,T) ~ ^ for t 1. 

Once we have known the radius (17.881) of truncation of the singularity structure of U at the 
origin, all computations can be redone to see its actual influence on the final asymptotics. 
In view of scaling (I2.18p . overall, this yields a formal expansion (I7.84p of the solenoidal 
vector field for a quadratic NSEs system, which is singular at the origin, but can admit 
some kind of critical surface inflection-like behaviour at the level set (T = 1) 

(7.89) H^^t)\-70^\ 

In other words, values (17.891) for t ^ T~, the singular vector field u{x,t) is expected 
to exhibit certain special transitional behaviour, which possibly can be used for further 
necessary matching and branching to create more reasonable blow-up patterns. 

The pattern with special level sets as in (I7.88P is obtained in the simplest situation, 
where no "resonance" between different terms of the expansion (I7.87p . which itself can 
generate logarithmic factors, is assumed. In general, for quadratic dynamical systems as 
(I2.2ip . there can be several types of asymptotic centre- manifold-type expansions, which 
are very sensitive and depend on resonance conditions of various terms involved; see invari- 
ant manifold theory in |145j . The validity of such conditions are very difficult to check, 
especially in the presence of extra outer matchings involved. However, in our opinion, 
(17.880 correctly characterizes some features of the blow-up behaviour of such patterns 
up to some extra exponents and/or slower growing factors. A justification of such a 
behaviour requires a quite involved matching analysis of fully using eigenfunction expan- 
sions on the vector solenoidal Hermitian polynomials and related delicate and technical 
matching procedures involved. 
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7.8. On a possibility of blow-up on regular equilibria. We now follow the discus- 
sion in Section II. 9[ These structures are more exotic but also reasonably well-known in 
reaction-diffusion theory; cf. [59] and [72]. For simplicity, we take a = in (17. 161) (or 
(17.151) . then a new perturbation occurs) and study a possibility of blow-up behaviour 

(7.90) u{y,T) = (y9(r)u((/?(r)y) + Y, where ipi^r) ^ oo as r — +oo. 
To this end, we first introduce the blow-up variables 

(7.91) u(?/,t) = V9(r)v(z, s), z = (p{T)y, ip"^ (t) dr = ds 
where v solves the following rescaled non-autonomous PDE: 

(7.92) v, = H(v)-p(s)[(^-V)v + v], where p(,) = + 1 J,) (r). 

If (T 7^ 0, we introduce the angular TW dependence p = /i + s, so that the extra term 
— crv^ should be put into the right-hand side (cf. (I1.32p ). which does not change the 
concepts of matching. It then follows that we may estimate v(s) as follows: 

(7.93) v(s) = u + p{s)W, provided that ^ ^ 0, s ^ +oo 

(i.e., p(s) is not exponentially decaying and is algebraic), where W is a proper solutions 
of the linearized inhomogeneous equation 

(7.94) Il'{u)W -[{z-V)u + u] = 0. 

In order to get possible acceptable families of functions {fkiT)}, we return to the 
original rescaled variable y and use (I7.90p in (I2.2ip . (12.230 to get for Y the equation 

(7.95) Y. = (B* - I /) Y + Ui^iT)ui^iT)y)) + ... , 

where we omit higher-order terms. By the assumption (I7.17p . we have from (17.21) that the 
main inhomogeneous term in (I7.95P can be estimates as follows: 

(7.96) H(^(r)ii(^(r)y)) = 5(y) + ... =^ = (B* - i /) Y - S{y) + ... . 
We then balance the terms similar to (I7.77P : 

(7.97) Yfc(r) =e(^'=-5)^v* + ... and Cfc(r) = e^^-^fe)^ for r > 1. 

Comparing (I7.90p . (I7.97P and (I7.9ip . (I7.93P yields the matching condition, which we 
take in the simplest form to catch the exponential factors only: 

(7.98) p(^) = ^ + i^~^e(^-i)^ =^ ^,(r)~(i-A,)e(^^'=)^-oo. 
Let us first check the consistency of the expansion (I7.97P for (17.910 : 

(7.99) s ~ e(*^+i)^ =^ pfc(s) ~ i for s > 1, 

so that each Pk{s) has the desired algebraic (and non-integrable) decay at infinity. 

Recall that these are blow-up patterns constructed via a hypothetical evolution on the 
manifolds of "steady" (or "TW-swirl" for a ^ 0) solutions, which includes both families 
of the singular S-L profiles and the regular ones as in (17.160 . Including the angular 
eigenvalue a ^ will make the formal analysis more complicated to say nothing about a 
rigorous justification (a finite energy interpretation of such blow-up patterns is also hard). 
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7.9. Towards Burnett equations. Finally, we note that, for the Burnett equations 
(11. 6p . a similar formal analysis can be performed relative to the singular solutions with a 
different behaviour near the origin, 

(7.100) uss(x)~^ as x ^ 0, 

though proving existence of such steady structures is a difficult problem, as well as a 
rigorous mathematical justification of the expansion and matching procedures for both 
families of singular and bounded states. Concerning the linearization approach as in Sec- 
tion 17.81 it can be performed in similar lines by using spectral theory and the generalized 
Hermite polynomials for the operator (13.241) . [51j. The possible ways of "interaction" with 
the spectral characteristics of the nonlocal operator J2 are out of question here. 

On the other hand, it is not excluded that (11.61) can admit a purely self-similar blow-up 
with various multiple patterns (a finite or countable set? - seems should be finite as 
higher-order parabolic flows suggest [22]; however, the solenoidal restriction can indeed 
easily spoil Leray's-type "similarity party" even here). Of course, the nonexistence proof 
based on the MP ideas from [162] is not applicable here. Thus, (II. 6p can admit much 
more complicated families of blow-up patterns (it seems that numerics can help here 
being however very difficult), though the classification problem (ll.20p remains a principle 
difficult issue, which, as expected, will never be completely solved (too difficult, and not 
that essential?). 

7.10. Complicated oscillatory singular equilibria for an R— D equation. Here, 
as a key illustration to some of our speculations, we briefly consider (ll.34p . In radial 
geometry, such Emden (1907)-Fowler (1914) equations have been most carefully studied 
in ODE theory. First detailed classiflcation of solutions were obtained in Gel'fand [SI] 
(this ODE section is known be written by Barenblatt) and Joseph-Lundgren |108] : see 
further references and applications for blow-up in [66, § 6.5]. We briefly comment on 
still difficult and seems not completely well-understood oscillatory properties of singular 
solutions. 

First of all, (ll.34p admits the homogeneous singular stationary solutions (SSS): 

(7.101) f/(r) = ±ar-^ where /i = ^, = [fi{N - 2 - fi)]^^ {p > j^) . 
To describe others of changing sign, we introduce the oscillatory component ip by 

(7.102) n(r) = r-'^(^(-lnr) : + (2/i + 2 - iV)(^' + ^(/i + 2 - A^)(^ + V = 0. 

Figure [T] shows the non-oscillatory character of singular solutions in the subcritical 
Sobolev range p < p^hy shooting from s = (r = 1). Both Figures [T]( a) and (b) explain 
that as r — »• 0, i.e., s +cxo, the oscillatory component stabilizes to constants ±C^ as in 
(I7.10ip . excluding countable set of regular patterns {uk, k > 0}, which are bounded at 
r = and have a fixed number of sign changes. Here (a) shows shooting the first profile 
Mo(r) > with no zeros, while (b) corresponds to shooting of ui{r) with a single zero 
on r G (0,1) (s G (0, -|-oo)). Surrounding those regular Mfe(r) are continuous families of 
singular solutions exhibiting the behaviour (17.1010 as r ^ 0. 
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Shooting first reqular solution: m=1 , p=4<Pg=5, N=3 Stiooting second reguiar soiution: m=1 , p=4<Pg=5, N=3 




10 20 30 40 50 60 70 80 20 40 60 80 100 

s s 

(a) shooting of uo{r) (b) shooting of ui{r) 

Figure 1 . Shooting the first (a) and the second (b) regular equiUbria via ODE 
in (m02|) for p = 4 < PS = 5, TV = 3. 

Next Figure [2] shows the crucial appearance of infinitely many oscillatory periodic orbits 
(p{s) at the critical value p = Ps, where fl7.102p takes the autonomous variational form 

(7.103) ip" - ^-^^ ^ + \^\P-^^ = 0. 

Finally, in Figure [3], we show how the oscillatory character of singular solutions dramat- 
ically changes for p > ps, where nonlinear spiral out behaviour replaces the bounded 
periodic one. In (a), we take p = 6 > ps = 5 for = 3. In (b), the highly oscillatory 
behaviour is shown for p = 2 and = 17, i.e., above the uniqueness [77] critical exponent 

(7.104) p* = i + __|-^ for iV>ll {p* = l<2 for iV=17), 

which plays a role in a number of blow-up problems for this R-D equation; see [77] and 
|146j for further details. Thus, the spiral-type oscillations for large p > ps get arbitrarily 
large, as r — > 0, so singular equilibria stay arbitrarily far from the homogeneous one 
fl7.10ip (is this possible for the steady NSEs (11.331) relative to the homogeneous S-L 
solutions fl7A|) ?l 

Thus, the singular equilibrium manifold for (ll.34p is rather complicated even in the 
radial case. Of course, in the radial geometry, the blow-up analysis of the corresponding 
parabolic equation (11.630 is essentially simplified by using the Sturmian argument of 
intersection comparison (the number of intersections of different solutions cannot increase 
with time; C. Sturm, 1836), and this can prohibit some scenarios of blow-up; see various 
applications in [BSl [77[ 1146] . In the non-radial geometry, even for (ll.63p . this advantage 
is almost nonexistent, the singular equilibria can get much more complicated and are 
unknown. On the other hand, the MP remains then in place, but its application to 
control, simphfy, or prohibit possible blow-up structures is unclear. 

Thus, the non-radial geometry in such parabolic problems can result in existence of 
new blow-up patterns with evolution "close" to such singular manifolds, where specific 
"swirl-torsion" blow-up phenomena may occur. Of course, the matching conditions such 
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Oscillatory component: m=1 , p=p =5, N=3 




Figure 2. Periodic oscillations in the ODE in (17.102P for p = 5 = ps, N = 3. 



Non -periodic oscillatory components: m=1, p=6>p =5, N=3 




J, ig''Non -periodic oscillatory components: m=1, p=2>p =9/5, N=17 




(a) PS <p<p* 



(b) p>p* 



Figure 3. Oscillatory singularities of the ODE in (|7.102p for p 
N = 3 (a) and for p = 2 > = |, iV = 17 (b). 



6 > PS = 5, 



as (17.571) for such singular equilibria may prohibit, via obvious divergence, almost all types 
of such blow-up patterns (the sum in (I7.56P is not from L°°), but we cannot always rule 
out a possibility of "accidental convergence" of the corresponding integral over T(r). 

We suggest that this R-D singular equilibria experience can be taken into account 
for the NSEs, and then the open problem (11.331) turns out to be key (does it admit 
local oscillatory singular solutions that have nothing to do with the homogeneous S-L 
one (I7.1l) ?-a principal question). Recall that regardless the presence of the standard 
Laplacian A in both (11.331) and (II. ip . the MP arguments cannot be somehow essential for 
such nonlocal problems (cf. (I2.2ip ). [All speculations are modulo (I7.57p .] 
Remark: singular equilibria in the bi-harmonic case and open problems. The 
analogous elliptic problem for the Burnett equations (11.61) is the bi-harmonic one. 



(7.105) 



u 



0, 
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which admits a similar to fl7.102p substitution with fi = and the ODE: 

_y,(4) _ j^^'" _ Bip" - Cif' -Dip+ |y?|P~V = 0, where 
A = 2(2^ + 4-N), 5 = ^ 18/i + 11 + (iV - l)(iV - 9 - 6/i), 
^' ^ C = 2[2^3^9/i2 + ll/i + 3 + (A^-l)[(iV-3)(/i+l)-3/i2-6/i- 
D = ij{fi + 2)[(/i + l)(/i + 3) + (AT - 1){N - 5 - 2/i)]. 

This is a harder equation than (I7.102p . and its complexity shows how difficult proofs of 
global or blow-up bounds on solutions of the corresponding parabolic PDE (11.561) . m = 2, 
can be for p > ps = with > 4. The homogeneous SSS of (I7.105P is 



^7 ^Q^^ Uir) = ±a r-^ where /i = ^, = > 0, 

existing for p > N > 4, or p < > 2. 



Global radial solutions of (I7.105P for p > ps are obtained in [83] , where further references 
can be found. 

Again, blow-up patterns can be created by the manifold of singular equilibria with 
completely unknown non-radial structure (with no traces of the MP); the analogy of the 
matching (I7.57P is supposed to be taken into account; see [Tl] for a discussion. Any 
general estimates of Type I or II blow-up for (11.560 for p > ps are absent. Moreover, 
we claim that a full and complete description of all the (non-radial) blow-up patterns for 
(I1.56P will never be achieved. Similar difficulties concerning existence and nonexistence 
of various blow-up patterns occur for supercritical nonlinear Schodinger equations such 
as (II. 9p (see key references on this subject in Merle-Raphael |15U] and Visan [215] ) and 
for many other important higher-order PDEs (see the list around (11.90 ) and systems of 
the twenty-first century PDE mathematics/application. 

8. On complicated blow-up patterns with swirl and precessions 

Using the previous ideas of possible concepts of blow-up in the equations (II. ip . we 
continue the construction of more refined structures of a full complexity. We now use 
spherical coordinates that allow us to fix more complicated singular stationary solenoidal 
fields, in a neighbourhood of which some non-steady blow-up phenomena may occur. 

8.1. Basic solutions with swirl in spherical coordinates. Consider the equations 
(II. ip in the spherical polar coordinate system, with u = {ur,uo,u^) = {u,v,w): 

' Ut + UUr + - VUq H wu^ — - — - nP' 

II / I J, 17 I rs,\nti V r r 

- r, J- A ^, 2 „, 2 „, 2cot6> „, 2 

- —Pr + — Z2U — -2V0 I 



sin 

^.1) { = -^pe + A3V - + ^ue- ^^w^, 

I • J. r-p I o r'^ sin y ° r^smS "i"' 
Wt + UWr + - + - VWe + VW H ^ WU 

II I ' y. Ij. 171^ Ij. gjjj ^ 

2 „, I 2 cot ( 



: P^ + A3W 2 -la W + 2 a Uu, + 2 a V, 



Ur + -u + -ve + —v + = 0. 
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Here A3 denotes the spherical M^-Laplacian: 



(8.21 



(sm9he)e + 



h 



Homogeneous singular stationary states. Similar to fl5.67p . we begin with simpler singular 
stationary solutions of (18.1 p of the homogenuity —1, 

(8.3) {u,v,w) = l{u,v,w){(p,e), p=^p{ip,e). 

This solution representation means that the operators of the stationary equations (18. ip . 
in view of their homogenuity, perform the following mappings of linear subspaces: 

(8.4) Wi = Span{i} ^W2 = Span{i, i}. 
Substituting yields the following PDE system for the four unknowns (u, v,w,p): 



1.5) 



■0^ - w'^ 



+ U00 + cot 6ug + 



1 



sin2 "-W 



vve + 



wv^ — cot 9w^ 



-Pe + vee + cot 9vg + 



1 



2u — 2vo — 2 cot 6v - 
- ^v + 2ue-^ 



sine ^V' 



VW0 + cot 6vw + ^ ww^ 



-^eP^ + ^oe + cot ewe + -^,w^^--^ 



■ w 



sine ^"P + sine 



u + V0 + cot 6v + 



Sine 



0. 



Indeed, the system looks rather frightening for studying in general. As we have mentioned, 
fortunately, it was proved in [204j that, up to an isometry, the only non-trivial (— 1)- 
homogenuity stationary solutions in \ {0} are the Slezkin-Landau ones (17.10 (it was 
also conjectured there that branching from usl is impossible). Further extensions, showing 
the exceptional role of the S-L solutions, are obtained in ^15'5j . In any case, (18.50 shows 



the range of typical difficulties concerning systems that inevitably occur while searching 
for other types of stationary or non-stationary singularity manifolds for linearization. 

Nevertheless, according to our "linearized strategy" of further construction of blow-up 
patterns (and for answering questions (II. 2p or (II. 3p ). 



(8.6) 



all properties of all solutions of (18. 5p should be known in detail. 



Indeed, singular stationary solutions (18. Sp . (18.50 can be key for a possible successful match- 
ing to create a blow-up pattern. The second step is then to answer: 



8.7) Questions (i)-(vi), § 15.121 for operators hnearized about all solutions of (18. 5p . 



To underline the complexity of the problems (18. 6p and (18.70 . we will stress below the 
attention to some known examples of particular solutions of that type. 

On general singular stationary states: elliptic evolution. Of course, the system (18.50 does 
not contain all the necessary singular stationary states for (18. ip . which might be important 
for blow-up constructions. The general representation of such solutions takes the form 

P = 72Piv,d,s), 
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where s = — Inr — > +00 as r ^ 0+ plays the role of the time- variable. Then 

(8.9) Ur = --^{us + u), Urr = :^{uss + ^Us + 2u), etc, 

so that the stationary system ( 18. ip . according to the variable s, takes the form of "elliptic 
evolution" equations for s ^ 1, 

(8.10) {Uss = ---, Vss = ---, Wss = ---, Us = .... 

In general regularity linear PDEs theory, such "blow-up" scalings lead to complicated 
spectral theory of pencils of linear operators, whose spectrum and root functions classify 
all possible types of singularities occurred. We refer to seminal Kondrat'ev's papers in 
the 1960s [1211 1122] and monographs by Maz'ya with collaborators [1231 1124j (further 
extensions via MathSciNet). 

As is well-known (since Hadamard's classic example), such an elliptic-like evolution is 
ill-posed and almost all of the orbits are destroyed before reaching the singularity point 
s = +00 (r = 0), but anyway the rest of the orbits that are defined for all s ^ 1 
describe all possible singular steady states for (18. ip . Such an analysis has been effectively 
implemented for a number of single semilinear elliptic problems, even in the case of non- 
Lipschitz nonlinearities; see e.g., a machinery and a full list of references in |T3]. Indeed, 
for the stationary system of four PDEs as in (18. ip . the problem of identifying all possible 
global evolution trajectories reaching s = +00 is extremely difficult. Recall again that 
each such a non-trivial singular stationary state with possible inclinations and precessions 
of the swirl axis (see further developments below) can be key for construction of a blow-up 
pattern by a some kind of a linearization and matching procedures. Of course, eventually 
we cannot escape the problem (18. 7p for each of those singular stationary solutions (18. 8p . 

In other words, the singular stationary problem (I1.33P plays a first key role for under- 
standing of the blow-up mechanism of the NSEs, to say nothing about a similar problem 
with rotations to be represented later on. 

8.2. Scaling and introducing blowing up angular singularity with precessions. 

We return to the general system (18.10 in the standard rescaled blow-up variables 

(8.11) u = ^u, p=jl-^p, y=^, r = -ln(-t), 

where, as one of the possibilities, we first assume the similarity swirling mechanism of 
accelerating rotation in the angle ip, 

(8.12) (p = fi-aln{-t)=fi + aT (a ^ 0). 

As we have seen, using the variable (18.120 introduces into the system an extra parameter 
cr e M being a nonlinear eigenvalue that increases the probability of successful matching 
of fiows in various regions. 
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Substituting flHUD . (KT2^ into yields the following system: 



' Ur + \ yuy + I M + aUf^ + uuy + ^ -O-ue + 

2 cote 



(8.13) 



2 



A,n- 



u 



Vg - 



V — 



y sm 
y'^ sin 9 



1^)2 



1 



sin 9 



y 

2 cote 



y sm e A* 



y 



cote ^2 



le "^M' 



t&T- + I + I w + crw)^ + uwy + ^ -uw) + ^ -fi + vw + ^ 



ysmf) A* 



1/ sm I 



■ + A3W 



j/2 sin-^ I 



y' sm( 



f, I 2 cot 



i/^sinS 



W,. 



0. 



V y y y y y sm f 

The rescaled Laplacian is now 

(8.14) A,h = ' {y'hy)y + 1 (sin 9 he 



1 

»2 sin^ e '"/^M- 



Recall that, in view of the ban (11.41) . for a = non-trivial self-similar, i.e., stationary, 
solutions of fl8.13p are in fact non-existent. More generally, we are supposed to perform 
a matching asymptotic expansion construction of blow-up patterns for (18.131) using some 
already known quasi- stationary manifolds. It should be mentioned the possibility of taking 
0" = for the (^-independence solutions 



(8.15) 



u = u{r,6,T), p = p{r,9,T). 



The long history of difficulties in proving global existence of solutions (18.151) [193] suggests 
that a blow-up pattern can be reveled even in this restricted geometry. 

However, it seems that a most rehable approach to blow-up patterns should comprise: 

(i) either the blow-up swirl mechanism (18.121) with a 7^ 0, or 

(ii) the asymptotically slowing down mechanism such as in (17.121) . (17.131) . formally 
corresponding to the case 



(8.16) 



a = a(r) ^0 as r ^ +00 in (Km . 



Both cases seem can suit the linearized construction about the Slezkin-Landau solu- 
tions (17. 6p or others more regular, which assumes a difficult spectral analysis of the linear 
operator as in (I7.18p . where the linearization is now performed relative to the nonlinear 
operators in (18.130 . Note that for the asymptotically cr = in any similarity annulus 
{ ^ < C\f^ }, with arbitrary C > 1, the construction remains the same and in- 

cludes general eigenfunctions of (17. 18p having a (^-dependence. Again, Questions (i)-(vi) 
from Section [5. 121 appear. Assume that we can answer these questions for some particular 
setting of singular quasi- stationary manifold. For instance, then as a by-product, this 
would actually mean that we may look for a point spectrum of B* with the eigenfunc- 
tion behaviour of the third w-component of the linearization (15.750 about the singular 
equilibrium, say, U = usl, 



(8.17) 



,Ai.r 



V^^(r, 6^, (y9) + ... as r — > +00, 



where G crp(B*), k = |/3|, is such that Re < (or simply A^ < for a real eigenvalue 
that also can be expected). In other words, the swirling blow-up behaviour then will 
occur on smaller shrinking subsets {r = ©(v^^)} as t — >• 0~, where one can expect both 
scenarios associated with f l8.12p or f l7.12p . f l7.13p . or with a constant, independent of ip 
rotation as for solutions (18.150 . 

Let us mention the following aspect from Section 15.91 can the possible zero mean 
oscillation property similar to (15.461) (in the original {x, t)-variables) affect the asymptotic 
behaviour of solutions as x — oo in such a way that the patterns will attain a finite 
energy? Actually, this is related to the hard problem of the asymptotic behaviour for the 
NSEs in a compliment of a bounded smooth domain, on the boundary of which blow-up 
swirl rotations with the zero mean as in (I5.46P is prescribed. Possibly, this could affect the 
asymptotics of the solutions which become better localized in the L^-sense. The difficulty 
of this question dramatically increases if a possible axis precession (see below) is taken 
into account. 

Thus, it seems that, according to the above scenarios, there is no chance to study 
possible admissible solutions of the resulting systems such as (18.130 in a reasonable and 
reliable generality and mathematical strictness. So we will continue to develop some 
necessary formal arguments in an attempt to give further hints for understanding such 
potential complicated blow-up patterns. As in Sections |5. 121 and [71 we will assume that in 
some intermediate region the blow-up behaviour goes along the quasi-stationary manifolds 
of singular steady solutions of (I8.13p . which, hopefully, have asymptotically the form (18. 3p . 
with r \—>- y, i.e., do not exhibit somehow essential dependence on the yj-torsion (which may 
get crucial for smaller y). Then, as we have seen, a complicated linearized operator such 
as in (15.761) and (17.181) occurs. The main difficulty is not studying its spectral properties in 
suitable weighted topology, but the matching procedures with bounded orbits for y ^ 0. 
For large y 3> 1, we always assume that there exists a possibility of matching the inner 
region with a properly deformed rescaled kernel of the fundamental solution of the operator 
Dr — A3, since all the evolution equations in (18.130 contain the necessary counterparts. 
Observe that the resulting patterns are not supposed to be of a simple self-similar form, 
so they do not exhibit any uniform homogenuity as in (18. 3p for steady profiles and/or 
L°^-boundedness as non-stationary similarity solutions after scaling (18. lip . 
An axis precession mechanism. In order to further increase the probability of such a 
matching of various manifolds for systems like (18.130 . in order to suit (I2.47p . it is neces- 
sary and natural for such swirling/vortex flows to introduce an additional axis precession 
mechanism for these blow-up twistors. This cannot be described explicitly or by a system 
on a lower-dimensional subspace. However, there is a standard asymptotic approach. 

Thus, we first need to plug into system another parameter in order to change the ax- 
isymmetric (with swirl) geometry of the solutions involved. As it has been mentioned, any 
strong symmetry (or quasi-symmetry) constraint in the Navier-Stokes equations reduces 
the dimension of the space, so approach them closer to the regular M^-case. Therefore, 
thinking about creating even more complicated geometry and not hesitating to perform 
rather weird transformations with such unknown solutions, we will show a way how to 
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perform a necessary perturbation of those solution structures. As a first simple formal 
illustration, let us take into account an extra slow motion in the ^-angle by setting 

(8.18) e = p + £T, where je] < 1. 

This will mean introducing into the system fl8.13p extra 0(£:)-operators, i.e., 

(8.19) ... + e-Up + ... , ...+eVp + ..., ...+ewp + ... into the LHSs of equations. 
Using in all the equations the corresponding e-expansions such as 

(8.20) sin6' = sin(p + £:r) = sinp + £:rcosp + ..., = st^ + ..., 

^ ' vr ' / r ' r ' 5 siny smp smp ' 

etc., gives extra entries of e into the system. According to asymptotic expansion theory 
(see Il'in |107j for typical difficult methods and further references) and not taking into 
account at this moment singularities introduced by f l8.20p into the PDE system, this 
makes it possible to look for solutions in the form of formal expansion^ 

(8.21) {u,v,w,p) = {uo,vo,wo,po) + e{ui,vi,wi,pi) + ... . 

By a standard procedure, being fully developed, the series f l8.2ip gives a unique formal 
representation of a certain formal solution. The main step is the first one, for e = 
that leads to the above nonlinear systems (which were assumed can be solved). The next 
terms {uk,Vk,Wk,Pk) are obtained by iterating some linear systems, that is easier, and 
existence of a solution can be checked by standard (but not always simple) arguments. 
The convergence of such series, as usual, is not required, and is often extremely difficult 
to guarantee even for lower-order PDEs, where the rate of convergence or asymptotics 
are also hardly understandable. Typically, such series serve as asymptotic ones, i.e., each 
next term correctly describes the behaviour of the actual solution as e ^ 0. 

Thus, fl8.2ip is now assumed to describe branching of suitable solutions from the un- 
perturbed one at e = with no fast swirl. Here, we refer to classic bifurcation-branching 
theory for equations with compact nonlinear integral operators, [T2[ HG], 11261 1212j . etc. 
Proving the actual branching is a deadly difficult problem, especially, since we are in- 
terested in quite special solutions only. Therefore, we do not concentrate on branching 
phenomena and continue to reveal possible features of such blow-up twistors. 

Note that fl8.18p indicates a tendency of the desired precession with potentially un- 
bounded deviation of the swirl axis. To avoid such an unpleasant (and seems non- 
realistic?) pattern, we assume that the actual precession is governed by (cf. fl7.12p ) 

(8.22) 9 = p + eK{T,e), where k(t, e) is bounded, 

so that fl8.18p contains the first term of the r-expansion, and extra difficult asymptotic 
matching theory occurs. 

Periodic and quasi-periodic axis precession. However, even the axis precession according 
to the simplified dependence f l8.18p can give insight into the actual behaviour of such 

'^"'^We discuss this simplified version of an e-expansion for convenience. In general, depending on the 
kernel of the linearized operator and other factors, the expansion can depend on the small parameter 
e^/', where integer Z > 1 is defined from the solvability of the corresponding nonlinear algebraic systems 
on the coefficients; see general bifurcation-branching theory, e.g., [212j . 
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matching blow-up patterns. In particular, as a formal illustration, ignoring at the moment 
a periodic-like behaviour in r, consider both angular dependencies of the swirl axis on the 
unit sphere §^ C ™^ 



(8.23) 



9 = eT. 



Since, as in the classic representation of periodic and quasi-periodic motion on a torus in 
M^, the angle behaviour is understood modulo 27r, so we can have both scenarios on the 
sphere Namely, assume for a moment that a matching procedure has turned out to 
be successful for a given pair of eigenvalues (cr, e). Then f l8.23p indicates that the periodic 
scenario of the axis evolution associated with fl2.46p (but not entirely; see below) takes 
place provided that 

(8.24) f G Q is rational, 
and, vice versa, we have a quasi-periodic precession if 

(8.25) I is an irrational number. 

The periodic scenario (I8.24p looks being the simplest one to create a blow-up singularity 
in the dynamical system (12.211) provided stabilization to a point is forbidden by (11.41) 
(if to follow the principle of Occam's Rasoi@)- But of course, this does not rule out 
other patterns, which actually, do not look more complicated and admit entirely similar 
characterization fl8.25p (though we do not exclude the non-precision case e = 0). For 
(18.221) . we can observe a more realistic scenario of very small precession exhibiting periodic 
or quasi-periodic motion of the axis. 

Briefly on precession of the vertex. It is quite natural that including both swirl and axis 
precession mechanisms will also require an extra "precession of the vertex" (as an ac- 
ceptable analogy with a rotating top on a sufficiently smooth surface suggests); otherwise 
such a complicated non-symmetric vortex evolution with swirl and precession having a 
fixed stagnation point would not be possible. This assumes the slow-variable change of 
the original rescaled coordinate system, 

(8.26) x^x + a(r), 

where hopefully a(r) can be a bounded function with also periodic or quasi-periodic 
behaviour as r — » ocO- This will involve into the systems such as (18. lip extra perturbed 
operators according to the change 

(8.27) ^ + (a' ■ V)u, 



W. Ockham's lex PARSIMONIAE: "entia non sunt multiplicanda praeter necessitatem" . 
'^'^For reaction-diffusion equations such as (|1.36p or (|1.60p . the blow-up set is introduced B[uo] = {xq S 
: 3 {xn} Xq and {t„} — > T~ such as |u(a;„,t„)| oo}. It is then proved that B is closed (e.g., 
is a point; see Friedman-McLeod [62j for a pioneering approach), and next blow-up scaling is performed 
relative to an xq € B[uo] by setting y — ^~ , etc. 
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so that these extra terms can be responsible for additional evolution blow-up phenomena 
(even in the asymptotically vanishing case a'(r) 0, since this also might support a quasi- 
periodic-like motion; the integrable case J |a'(r)|dr < oo makes (18.261) non-essential) 
similar to those induced by swirling (17.121) and precession (I8.22p . All three mechanisms 
taking altogether lead to a possibility to attempt to construct a blow-up twistor pattern 
originated in the outer region by the Slezkin-Landau singular steady solutions (17.61) or 
others, which possibly are still unknown. 

In conclusion, let us mention that we do not stress attention to a possible more physical- 
mechanical interpretation of the blow-up patterns under speculations. Namely, we do 
not know and cannot imagine how many actual twistors should be involved in such a 
"hypothetical fluid configuration" to produce such a pattern under the fixed blow-up frame 
and zoom. We just recall that, as a consistent part of our "swirling-like philosophy", we 
mean that complicated blow-up patterns can be locally trapped as t — T~ by the quasi- 
stationary singularity manifold described by the problems (ll.33p . (I1.32p . and others. Such 
a blow-up drift along those singular manifolds will itself define the type of generalized 
swirling that is necessary to support the evolution. 

Mechanical interpretations of various solutions and patterns have always been very 
effective, difficult, and involved techniques of modern applied and mathematical fluid 
dynamics, which can be also efficient after a blow-up scaling, but possibly could fail in 
view of the fact that the blow-up micro-structure of the NSEs might have nothing to 
do with their classic macro-coherent structures studied during almost two centuries. A 
possibly acceptable example is as follows: the local atomic (micro-scaled) structure of a 
desk, where this paper is about to be finished, has nothing to do with its global (macro) 
properties. This is up to the obvious fact that here the micro- and macro-mechanics are 
different: the quantum and Newton's ones. A similar phenomenon occurs for the NSEs: 
the operators of the original PDEs (ll.ip and of the blow-up rescaled ones (12.210 live in 
completely different spaces, L^(M^) and L^, (M'^). 

Again on L°° -hounded or unbounded rescaled orbits {u(r)}. Finally, we add to this formal 
description of possible blow-up patterns of twistor type the following necessary observation 
based on our previous analysis. Namely, according to matching concepts revealed in 
Sections 15.111 [5rT2| and [71 there exist two cases (sub-scenarios relative to the above): 

(i) Type I: the rescaled orbit {u(r)} is bounded in as r +oo. Then 
this corresponds to scenarios (I2.46P and (I2.47p . as before, and 

(ii) Type II: ||u(r)||oo — > oo as r +oo. This can happen according to the matching 
analysis and can be driven by an exponential (I5.88P or power-like as in (15.900 divergence 
(or others) depending on the manifold that was used for matching purposes. Of course, 
this does not ruled out the quasi-periodic behaviour (15.120 of bounded orbits in the 
necessary new rescaled variables. 

In both cases, the blow-up patterns remain exhibiting similar properties of swirl and 
precession but, possibly, on different spatio-temporal subsets. 
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9. Final remarks 



9.1. Micro-structure of turbulence. As the attentive Reader has noticed, the main 
goal of the present essay is not about achieving or even essentially approaching a definite 
answer to the fundamental open problem: fll.2p or (11.31) for the NSEs f ll.ip . Actually, this 
is more about approaching better understanding the Goal (11.201) . In other words, this is 
about presenting some ideas on a description of 



(9.1) existing micro-scaled fiuid configurations appearing from smooth data. 



We recall that even for bounded smooth solutions, revealing such micro-structure of mul- 
tiple zeros at regular points (Section [3]) led us to some difficult problems, though their 
Hermite polynomial solenoidal structure was partially justified. 

Concerning singular points, in (19. ip . the term "micro- scaled" is key meaning to look for 
fiuid configurations that are seen via microscopic "blow-up rescaling zoom" as in (l2.7p F^. 
This means a description and a classification of the admitted "local turbulent micro- 
structures" of the Navier-Stokes equations (II. ip . The problem of "micro-structure" can 
be posed for any linear or nonlinear evolution PDE of parabolic, hyperbolic, nonlinear 
dispersion, etc., types, and gets very complicated even for simple models (cf. [SH § 9] as 
a short introduction to this involved subject). 

As a by-product, solving (19. ip would also mean the negative answer (II. 30 provided that 
the family of those configurations would include some blow-up patterns. Not pretending 
at all to giving any comprehensive insight into the problem (19.10 . we just have shown 
that including a standard similarity "log-torsion" mechanism produces a number of very 
complicated singular stationary (cf. (18.60 ) or evolutionary dynamical systems as some 
lower-dimensional reduction of the NSEs. This also implies that a detailed study of such 
reduced, but still very complicated, dynamical systems and corresponding equilibria are 
necessary and unavoidable steps that should have been passed before even thinking about 
attacking the Millennium Problem (along the proposed lines). 

9.2. A final pessimistic expectation in general PDE theory. As a consequence of 
all the above blow-up speculations, we first state the simplest claim: for the NSEs, 



(9.2) from the side of (II. 3p . checking all blow-up configurations is impossible. 



More precisely, involving the positive part (II. 2p . it can be emphasized that, for a suffi- 
ciently wide class of complicated dynamical systems "M x {N + 1)" such as (12.210 with 
M = 4, = 3 (the numbers of dependent and independent variables involved), with a 



^■^In general, this reflects a small part of fluid dynamic problems of fundamental importance; e.g., cf. 
A.M. Lyapunov Master's Thesis "On Stability of Spheroidal Equilibrium Forms of a Rotating Fluid" 
(S. -Petersburg University, 27th January, 1885; Supervisor: P.L. Chebyshov), which was a forerunner of 
Lyapunov's stability theory and a starting point for his correspondence since 1885 with H. Poincare, who 
also in 1885 obtained a linearized system of Eulcr equations about a rotating fluid pattern as a rigid body 
around the z-axis; see |179| for an extra account. 
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similar mathematics including a global energy control of solutions (not enough to guar- 
antee L'^-bounds by embedding and/or interpolation), divergence of operators, scaling 
laws, and other related and necessary propertie^^. 



(9.3) a definite answer to claims like (I1.3P and (ll.2p . in general, is impossible. 



As we have seen, for the NSEs (11. ip . some ideas of blow-up focusing can be understood 
via standard asymptotic language, though a full justification could also take years or do 
not admit such at all. For the fourth-order bi-harmonic operator as in (11.61) . being also 
parabolic but with no order-preserving and nonlocal properties, a similar proof often can 
be characterized as being completely illusive. Moreover, here a self-similar blow-up is then 
rather plausible as for (ll.56p . m = 2, [22], though will be extremely difficult to prove. In 
other words, there is a huge probability that 



(9.4) Problem: (II. 2p or (11.30 . can be analytically non-solvable. 



Rephrasing the above, problems of such complexity from the side (11.21) can be "accidentally 



solvable'ij, when special tricks associated with these equations and operators only (not 
robust techniques admitting perturbations, i.e., "structurally stable") can rule out some 
essential part of the core difficulties; all modulo (11.221) . But this "accidental" feature 
is rather unlikely: both (11.21) and (11.30 are too much related to each other, and proofs 
should pass through a lot of similar extremely difficult stages. We should accept that, 
in modern nonlinear higher-order PDE theory, many standard results, which had been 
perfectly solved for lower-order counterparts by inventing great mathematical methods by 
great mathematicians, do not admit rigorous setting in principle. In this case, looking and 
searching for results admitting rigorous proofs would be a wrong idea, also contradicting 
Kolmogorov's thoughts; see a footnote in Section 11.41 

To this end, we recall that a definite, positive or negative, answer in particular assumes 
a rigorous checking whether (12.470 is true or not for any n > 2. Let us also remind 
that, unlike (I8.23p . in general, we are talking about a quasi-periodic motion in an infinite- 
dimensional functional space for u(r), so that possible localization of that is a very difficult 
problem. Recall that such quasiperiodic motions can be trapped in a vicinity of singular 
"equilibria" (11.321) or (I1.33P of very complicated unknown structure. In its turn, via a 
matching approach, this means checking whether the corresponding asymptotic bundles do 
or do not overlap for any n = 1, 2, 3, ... . This creates a restriction that cannot be checked 
analytically in general because it depends on unknown and unpredictable conditions of 
matching of various pair of infinite-dimensional local vector bundles. In view of the 
energy control, those matching look like checking if, for a suitable 2D restriction, there 



Just in case, modulo (|1.22p ; recent almost purely "parabolic" and MP-type results in |118l [55] (I 
am sure that further, even stronger papers in this directions will follow soon) for nonexistence of Type 
I blow-up for axi-symmetric flows inspired some optimism, though, as we have tried to show in Sections 
15.2115.31 and [3 this can be just the beginning of a very long road to success. 

^^Or, at least, it then could be not an exaggeration to comment that, posing this Millennium Problem, 
it would be also useful to specify, which millennium it is supposed to be. 
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exists a heteroclinic path connecting two saddles for a given dynamical system on a 2D 
manifold. For general DSs, existence/nonexistence of such a path cannot be predicted, 
and moreover, DSs with such heteroclinic orbits are known to be structurally unstable 
(Andronov-Pontriagin-Peixoto's theorem, 1937-57). It then seems that such a matching 
is not possible for almost all DSs like that. This is true, but since, according to fl2.47p . 
at least a countable (or more than that) number of such possibilities is supposed to be 
checked, this changes the dimension of the parameter space and, eventually, makes the 
problem to be analytically non-solvable. 

If, in reality, the justification of the negative claim (11.31) would have been proved of being 
of a geometric "configuration" as a matching of the type saddle-saddle on some manifold 
for establishing existence of a blow-up pattern with finite energy, it would clearly suggest 
that any enhancements of functional space and corresponding facilities for proving (11.21) 
by more and more refined interpolation-embedding techniques or similar would be entirely 
hopeless. Unless there is a parameter "gap" prohibiting such matchings uniformly in n > 1 
in (I2.47P and for other types of patterns. Hence, the existence problem takes a principally 
other background and becomes the question of blow-up scaling. 

In other words, then the NSEs regularity problem should be classified as being of "point- 
wise sense", i.e., its global solvability is not controlled by any possible a priori bounds 
in or related Sobolev spaces, so it is hopeless to try to derive those. Nowadays, more 
and more nonlinear evolution PDEs (for instance, supercritical nonlinear Schrodinger- 
type or sufficiently multi-dimensional Burnett equations) penetrate into this class, which 
demands principally new mathematics of the truly twenty first century. One should be 
ready to recognize that, in this PDE class, most of typical mathematical problems will be 
analytically non-solvable. This by no means diminishes the role of the pure mathematics. 
On the contrary, this implies that the whole mathematical culture will be needed to built 
a necessary well-organized understanding of the problem, under the pressure that no even 
a hope for any definite rigorous answers exists. 

Eventually, the author quite bravely expresses his personal opinion, which in view of 
the possible feature (19.41) might have some sense: 

(9.5) for (II. ip . global existence is more plausible than blow-up. 

Actually, this is a pure probability based upon above long speculations and discussions. 
However, global existence issues, though expected to be somehow related to blowing up 
ones, were not the subject of the paper. 

Acknowledgments. The author thanks C. Budd for a short but stimulating discussion. 
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Appendix A. English Translation of N.A. Slezkin's paper [196] 

ON AN INTEGRABILITY CASE OF FULL DIFFEREN- 
TIAL EQUATIONS OF THE MOTION OF A VISCOUS FLUID 

N. A. Slezkin 
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[Scientific-Research Institute of Mathematics and Mechanics] 

If the motion of viscous fluid is stationary and the axisymmetry of the flow takes place, then, 
as is known, the stream function satisfies the following differential equation in the cylindrical 
coordinates: 

(A.l) -l^dD± l^dD± 2^^jj^^^jjjj^ 

\ J r oz or r or az az ^ t i 

where 

D = -- — + 

dr '-^ r dr dz'^ ' 

or, on introduction of the conical coordinate p and cos 6 = t, the equation (jA.lh takes the form: 
where 

D = 9^ + i-^'' 

Under a certain assumption on the form of the function -0, equation (jA.2p can be reduced to an 
ordinary differential equation of fourth order, and the latter to a Riccati equation. Indeed, set 

(A.3) ^ = pfir). 

Then we have: 

and the equation (|A.2|1 takes the form: 

(A.4) ff" + 3/7" = - T^)f'' - 4Tf"']. 

The left-hand side of this equation can be represented as 

//'" + 3/7" = 

while the right-hand one as 

u[{i-T')f'+2Tfr, 

and then (|A.4p is rewritten as follows: 

(A.5) (i/2)'" = ^[(l_^2)^, + 2r/r. 

Integrating it three times, we obtain: 

(A.6) i f - - T^)f' + 2r/] = Co + Cit + Car^, 

where Co, Ci, C2 are constants of integration. 
Equation ()A.6P can be rewritten as: 

Thus, we have obtained a Riccati differential equation of the form that is not explicitly 
integrated. On substitution 

this reduces to the linear differential equation of the 2nd order: 



(A.8) ^ + Co±c^i^y = o, 
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The solution of the latter one can be studied by using methods of analytic theory of differential 
equations. 



[Uchenije Zapiski MGU, No. II, 1934.] 



UBER EINEN INTEGRIERBAREN FALL DER VOLLSTANDIGEN BEWE- 
GUNGSGLEICHUNGEN EINER ZAHEN FLUSSIGKEIT 

N. Slioskin 

[Institut fiir Mathematik und Mechanik] 

(Zusammenfassung) 

Es wird der Fall untersucht, in dem die Gleichung fiir die Stromfunk tion sich auf eine 
gewohnliche Gleichung, und diese letztere auf eine Riccatische Gleichung zuriickfiihren lasst. 

[Wissenschaftliche Berichte der Moskauer Staatsuniversitat, H. II, 1934.] 



Appendix B. English Translation of N.A. Slezkin's paper [197J 

Remark on the notes of Yu. V. Rumer, "The problem of a submerged jet" 

and of L.G. Loitianskii, "Propagation of a whirling jet into an infinite spaces filled 

with the same fiuid"'^] 

N. A. S 1 e z k i n (Moscow) 

In both notes, it was pointed out that the problem on a laminar submerged jet was first 
considered by L.D. Landau (Mechanics of Continuum Media^'^l) and that the solution of this 
problem represents a new case of explicit integration of equations of motion of a viscous fluid. 
The note of such kind in the second part is not correct. 

In the note by N. A. Slezkin "On an integrability case of full differential equations of the 
motion of a viscous fluid" W published in Uchenije Zapiski MGU, No. II, 1934, it was shown 
that the full differential equation for the stream function of a steady axisymmetric motion in 
the spherical coordinates r and r = cos 9 

D = + 1-^' 
under the assumption ^ = rf{T) reduces to the equation 

fi _ 1 f 2 _ 2r f I Co+CiT+C2T^ 
■I ~ 2i/(l-r2) J J ^ 1-t2 
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was reduced to the linear differential equation 



Equating the constants Co, Ci, C2 to zero, we arrive at the solution 



/ 



1 



A+T 



which was obtained by L. D. Landau within the study of the problem of a submerged jet. This 
case of integration of full differential equations of the motion of a viscous fluid was referred to 
in L. I. Sedov's book^^l (p. 104) and in Rosenblatt's paper^^]. 
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Appendix C. On some simpler exact singular steady solutions: more history 

For further application to blow-up patterns, we will need to revise some stationary solutions 
of (jl.ip that are singular (unbounded) at the origin x = 0. As above, such singular equilibria 
can be used for attempting to construct different blow-up patterns that are not bounded in 
the rescaled similarity variables as in (j2.7p . Constructing various exact or explicit particular 
solutions of the NSEs is a classic, effective, diverse, and, at the same time, difficult way of 
understanding fluid flows. Some basic ideas go back to Oseen and von Karman. There are 
various other approaches to constructing other, exact or explicit, less and non-singular solutions 
of the Navier-Stokes equations; see monographs [133^ I173j and surveys in [86l I167j . with a 
number of references therein. 

Singular solutions (j8.3p of the homogenuity —1 have a long history. Actually, Slezkin-Landau's 
explicit solutions (jT.ip belong to the same class. As we have mentioned, the first ODE reduction 
was due to Slezkin |196j in 1934. For convenience of the Reader, in Appendix A, we present 
the English translation of Slezkin's short note |196] . while Appendix B contains the translation 
of his further note [197] (1954), so both rather convincingly show that Slezkin's ODE analysis 

^'''clearly, a misprint: should be (1 — t^).-VAG. 
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includes the crucial explicit solution. Later, Landau [132] in 1944 derived a particular explicit 
solution using a deep physical motivation and interpretation associated with submerged jet (the 
derivation in Landau-Lifshitz p^33j pp. 81-83] well corresponds to the that in the original paper 
|132j ). Some other aspects of derivation of Slezkin-Landau solutions can be found in Sedov [190^ 
pp. 99-104] and Batchelor [Ql p. 206]. 

It seems that, before Slezkin, the first attempt to construct exact stationary solutions was 
performed by Strakchowitsch (1931) |201j by looking for stream functions in the forms 

2 

= x('") + and ^ = x{i~) + 

Based on these Strakchowitsch's and Slezkin's |196j representations, Rosenblatt (1936) [185] 
constructed solutions with 

Both equations admit straightforward integration ones. It was shown that the equations admit 
non-singular solutions given by infinite power series. Therefore, it seems that Slezkin |196| was 
the first who looked for exact stationary fluid flows that are singular at the origin x = 0. 

In 1950, Yaceev [218] performed a detailed construction (in fact, quite similar to Slezkin's 
one; cf. Appendix A) of such exact solutions without torsion and u) = 0, where 

(C.2) n = ^, . = ^, ^ = 0, p=^. 



These assumptions essentially simplify the system (|8.5p that now takes the form 

{-u^ + vv! -v"^ = 2p + v" + cot eii' -2v -2v' -2 cot Ov, 
vv' = -p' + v" + cot eij' - -irg v + 2u', 
u + v' + cot9v = 0. 

The last equation gives the first velocity component 
(C.4) u = -v -cotOv, 

and substituting into the second one, on integration once, yields the pressure, 
(C.5) p = -v' -coteij-^v^ + Co (CoGM). 

Finally, substituting into the first ODE yields the following solutions of (|C.3p [218| : 
(C.6) i, = -y'-,otev, p = -2f)' + ^(^^^, v = -%^, 

where x(^) is given by the hyper-geometric functions 

x(e) = (cos|)-(smf)^+"+^- 
X [ciF(a, /3, 7, cos^ |) + c2F{a + 1 - 7, /3 -M - 7, 2 - 7, cos^ |)] , 
where ci^2 G The relations between constants a, b, c and a, (3, 7 are as follows: 

r a = 72 - (1 + a + /3)7 + i (a + /?)2 - i, 
1 6 = (a + /3 - 1)7 - i (a + /?) + i, c = i [(a - /3)2 - 1]. 



(C.8) 
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The hyper-geometric function F = F[z) in ()C.7p satisfying the ODE 
(C.9) z(l - z)F" + [c - (a + 6 + l)z\F' - abF = 0, 

is regular at the origin and is given by the Kummer power series, 
(CIO) F(a, 6, c, z) = 1 + f z + i ""^""tlc+T^ + ■■■ ' 

where c 7^ and c 7^ — /, with / S N. This converges uniformly in {\z\ < 1} and also on the unit 
circle {\z\ = 1} if a + 6 — c < 0. Note that the reduction to a linear ODE related to (|C.9p was 
already available in Slezkin |196j . It follows from the expression for p in ()C.6P that sufficiently 
regular solutions demand a = b = 0. The solution (jT.ip in these coordinates yields ()7.6p and is 
a particular case of (1C.6P : see |132j and [1331 p. 82]. 

Similar exact singular solutions were obtained by Squire (1951) [200j (see Pai |173t p. 72] for 
extra comments), where (cf. (IC.2P ) 

(C.ll) ^ = £S> v = -J§„ w = 0. 

Then, similarly to ()C.4p . we have from the second equation in (18. ip that 

(C.12) pg = -yyg + l.ue =^ p = -^v^ + lu + ^ (ciGM), 

where the constant of integration is taken in the form ^ for further convenience. The first 
equation in (j8.5p then yields for / = f{C), with ^ = cos9, the following equation: 

(C.13) (/')' + //" = 2/ + [(1 - e)f"]' - 2ci. 

Integrating leads to the first-order quadratic Bernoulli- type equation = 4^/ -|- 2(1 — C^)/' — 

,2,3,) 



2(ci^^ + C2C + C3), so that (a, /3, b depend on ci^: 



(C.14) fiO = a{l + e) + /3(1 + + '-^^^^ [b - If drj] -\ 

The regular case a = (3 = can be interpreted as a jet issuing from a nozzle. 

A further asymptotic extension of the exact solutions ()C.2p was performed two years later by 
Rumer [186j . who looked for solutions on the subspace W3 = Span{i, t-, ^-j. 



(C.15) n = ^ + ^, . = ^ + ^, w = 0, p=^ + ?). 



It turns out that a certain extension of solutions (|C.6p . (IC.7P exists in the case a = b = c = 
in (|C.8p . but these are not exact solutions so that (|C.15p gives the next asymptotic term of 
expansion in terms of K Further, Loicyanskii |143j in 1953 proposed to look for solutions with 
torsion in the cylindrical coordinates {r, ip, z} in terms of power series for u, w ^ small, 

(C.16) {u,v,w): En=i^^ where r? = f , 

which are not converging and are singular at the origin. The formula ()C.16P may supply us with 
some other singular solutions, that, though not admitting explicit representations, can be used 
in blow-up analysis. In general, it seems that the method of formal asymptotic expansions is 
an appropriate way to treat our more general and difficult system (|8.5p and further related ones 
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to be introduced. We refer to Kurdyumov |127j for more recent developments associated with 
such exact solutions of ()l.ip . 
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